Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,359,314 G→A Q240* (CAG→TAG)  cadB ← lysine:cadaverine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,359,3140GA92.5% 155.2 / 2.7 54Q240* (CAG→TAG) cadBlysine:cadaverine antiporter
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (15/34);  total (18/36)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

TTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGTCATC  >  NC_000913/4359070‑4359555
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
ttttccgaaaTACCGTTGCGTTCGACTTCACCATAAACTTTCGGGAGGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCg                                                                                                                                                                                                                                              <  1:314294/242‑1 (MQ=255)
            ccGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAAc                                                                                                                                                                                                                                 >  1:420313/1‑251 (MQ=255)
            ccGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAAc                                                                                                                                                                                                                                 >  1:411812/1‑251 (MQ=255)
                       gtcTTCTCCATAAACGTCCGGGAAGTTACCGTCGTTACCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGGCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATATTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTa                                                                                                                                                                                                                      <  1:393451/249‑1 (MQ=255)
                           cccccaaaaaCTTTCGGAAAGTTACCGTCGTTCGCGGCCCGTACACCTGCCTGGCCTACAAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAAcc                                                                                                                                                                                                                  <  2:410306/244‑1 (MQ=255)
                              ccATAACCTTTCGGGAAGTTACCGTCGTTAGCGGCCCGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGt                                                                                                                                                                                                               <  1:102273/251‑1 (MQ=255)
                                         cGGGAAGTTCCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATg                                                                                                                                                                                                     <  1:12348/250‑1 (MQ=255)
                                         cGGGAAGTTACCGTCGTTAGCGGCACGTACACCGGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATg                                                                                                                                                                                                     <  1:385268/250‑1 (MQ=255)
                                          ggtaggTTACCGGCGTTAGCGGCGCGTACACCTGCCGGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCGGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGt                                                                                                                                                                                                   <  1:101537/246‑1 (MQ=255)
                                          gggAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTAGCCTACCAACATCGTCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCAGTGAATGCAGAAACCAGCGGCGAAGCCCAGATACCGCGGATAGTCAAAGAACAGATTGCGAACGACACACCGAGAGCCGCCATTACAGAAGACGGATACATACACGACAGCACCTCAGACCCAGAGTGGTAAACAATACCTGCTAAACCAGTACCCAACAtcgt                                                                                                                                                                                                   >  2:257412/1‑248 (MQ=255)
                                               gTTACCGCCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGGGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCa                                                                                                                                                                                              <  1:337267/251‑1 (MQ=255)
                                                   ccGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACATGAGTCGCAGCGATGTAAACAATACCTGCTAAACCCGTACCCAGCATGGTTGCCAGCgg                                                                                                                                                                                          >  2:181173/1‑251 (MQ=255)
                                                              ggCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGtt                                                                                                                                                                               <  1:396352/251‑1 (MQ=255)
                                                                           gcctgGCCTACCAACATCATCCAGGAGCCCAGGAAAGTCAGGCACGCAAAGGCGGTGAAGGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTAGAAGCCCTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTa                                                                                                                                                                    <  2:236598/249‑1 (MQ=255)
                                                                           gcatgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTaa                                                                                                                                                                   <  1:88409/247‑1 (MQ=255)
                                                                            cctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAAcc                                                                                                                                                                 <  1:181173/251‑1 (MQ=255)
                                                                                               tccgatccccAGAGAAGGCGGGCACGCAAAGGCGGGGAAGGCAGAAACCAGCGGCGCAGCCCAGTTACCTAGGATAGTTGAAGCACTGATTGAAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCt                                                                                                                                               <  2:180221/243‑1 (MQ=255)
                                                                                               ccAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCt                                                                                                                                               <  1:444904/250‑1 (MQ=255)
                                                                                                     gCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGAtt                                                                                                                                        <  1:400788/251‑1 (MQ=255)
                                                                                                       ccAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCTAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCt                                                                                                                                               >  1:180221/1‑242 (MQ=255)
                                                                                                           agaggTCAGGCACGCAAAGTCGGTGAGTGCAGAACCCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAacac                                                                                                                                  <  1:536901/247‑1 (MQ=255)
                                                                                                               gTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCAc                                                                                                                              <  1:240872/251‑1 (MQ=255)
                                                                                                                    gCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAgg                                                                                                                         >  1:313605/1‑251 (MQ=255)
                                                                                                                       cGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTCCCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGccc                                                                                                                      <  1:215786/251‑1 (MQ=255)
                                                                                                                          aaaGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCcac                                                                                                                    >  1:148945/1‑250 (MQ=255)
                                                                                                                            ggggtgggagtgctgaacccaccgacgcAGCCCAGTTACCTAGGATATTTGAGTAACTGATGACAAAGGGAGCCCGGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAgg                                                                                                                 <  2:485135/226‑1 (MQ=255)
                                                                                                                             ggCGGTGAATGCAGAAACCAGCGGCGCAACCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAACCGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGgc                                                                                                                <  1:543118/251‑1 (MQ=255)
                                                                                                                                   aaaTGCAGAAACCAGCGGCCCCGCCCAGTTTCCGAGGATAGTTGAAGCCCTTATTGCAACCGGAGCCCCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATGCCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGgcagagca                                                                                                          <  2:76612/250‑1 (MQ=255)
                                                                                                                                       gCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGaa                                                                                                       <  1:262557/250‑1 (MQ=255)
                                                                                                                                        cAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATa                                                                                                     <  1:251437/251‑1 (MQ=255)
                                                                                                                                         aGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCCCCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATAc                                                                                                    <  1:220043/251‑1 (MQ=255)
                                                                                                                                          aaacccAGGGGCGCGGCCCTGTTACCGCGGTAAGTGGAAGCACTGATTGCAAACGGAACCCCGGAAGCCGCCATTACAGAAGACGGATACATCCCGTAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACt                                                                                                   <  2:374442/247‑1 (MQ=255)
                                                                                                                                           aaaCCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACAACCACGAAGGCCCACAGACCGAGCAGAATaatt                                                                                                  >  2:455765/1‑248 (MQ=255)
                                                                                                                                             aCCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGACCACAGGCAGAGCAGAATACtttt                                                                                                <  1:275844/251‑1 (MQ=255)
                                                                                                                                                 gcggcgCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATg                                                                                            >  1:252360/1‑251 (MQ=255)
                                                                                                                                                  cggcgCAGCCCAGTTACCGAGGATAGTTGAAGCCCAGATTGAAAACGGAACCCCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGa                                                                                           <  2:499727/251‑1 (MQ=255)
                                                                                                                                                       cAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAATCCGCCATTACAGAAGACGGATACATACCGGAAAGAACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCGAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGAAGAGCAGAATACTTTTAAGGATCGc                                                                                       >  1:286123/1‑250 (MQ=255)
                                                                                                                                                        aGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGCAGCCGCCATTACAGAAGACGGATACATACCGGCAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGAACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCAt                                                                                     >  2:202097/1‑251 (MQ=255)
                                                                                                                                                                        agaGTTGAGGCCCTGATGGCAAACGGGGCACCAGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATc                                                                     <  2:260896/249‑1 (MQ=255)
                                                                                                                                                                                agaacggATTGCAAACGGAGGTCCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTAt                                                             <  1:514497/245‑1 (MQ=255)
                                                                                                                                                                                  cACTGATTGCAAACGGAACCCAGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATTACCATCAGTGGTATCCGCAGTATTc                                                           <  1:145904/251‑1 (MQ=255)
                                                                                                                                                                                      gATTGCAAACGGAGCACCGGAACCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGt                                                       <  1:273785/251‑1 (MQ=255)
                                                                                                                                                                                                     aCCGGAAGCCGCCATTACAGAAGACGGATCCATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCCGCTAAACCAGTACCCAGCATGGTTGCCAGCGGCACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCAAAAGTTg                                        >  2:264787/1‑251 (MQ=255)
                                                                                                                                                                                                         gAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGc                                    >  1:203944/1‑251 (MQ=255)
                                                                                                                                                                                                          aaGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCca                                                                                                                     <  2:151903/169‑1 (MQ=255)
                                                                                                                                                                                                          aaGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCca                                                                                                                     >  1:151903/1‑169 (MQ=255)
                                                                                                                                                                                                             ccgccATTACAGAAGGCGGATACATACCGGACAGCACCTAAGCCGCAGCGATGTAAACAATTCCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACGCCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCa                                <  2:130873/251‑1 (MQ=255)
                                                                                                                                                                                                                          agaCGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATgccagc                   <  1:435183/251‑1 (MQ=255)
                                                                                                                                                                                                                                atacattcCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCCGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTTGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCcaacaa             <  2:411812/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tacCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCGGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGc         >  2:411093/1‑250 (MQ=255)
                                                                                                                                                                                                                                     tacCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGc         >  1:372843/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ggAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGCATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGtc     >  1:40495/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gAAAGCACCTAAGTCGCAGCGATGTAAACAATCCCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGtcat   <  2:420313/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGATGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCTACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAAAAGCAGtcatc  >  2:5028/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
TTTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGTCATC  >  NC_000913/4359070‑4359555

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑