Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,616,668 A→G 17.4% M193V (ATG→GTG)  bepA → OM protein maintenance and assembly metalloprotease and chaperone, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,616,6680AG17.4% 120.4 / 4.7 51M193V (ATG→GTG) bepAOM protein maintenance and assembly metalloprotease and chaperone, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (4/37);  new base G (8/2);  total (12/39)
Fisher's exact test for biased strand distribution p-value = 2.94e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGA  >  NC_000913/2616444‑2616906
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
tttgccTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGGTATGGCGCCCGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCAcc                                                                                                                                                                                                                                    <  1:752437/235‑1 (MQ=255)
 tGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGGTGGCGGGGCTTACCGGTGCACTGGCGGGGGA‑‑GCGGCCGGGGAGGATGAGATTGAccc                                                                                                                                                                                                                                   >  1:629791/1‑239 (MQ=255)
 tGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACCACGCCGGGAAC‑‑CGGGCAAGGGGGGGTCAGGTTCAcc                                                                                                                                                                                                                                    >  1:933600/1‑238 (MQ=255)
       gTTCCGTTATGCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGGCCCCCACCGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCaaa                                                                                                                                                                                                                             <  1:212290/239‑1 (MQ=255)
        ttCCTTTTTTCCGATAACGAAAGTCAACGGGCTTCAGGTTTGGCGCCAGAAATCCCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCTGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCaaaa                                                                                                                                                                                                                            <  2:933600/239‑1 (MQ=255)
        ttCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGGTCTATTTTACTGGCGATGGCCAGTCCGCGGGCGGGGGTGGCGGGCGTGACCGGGGCCATGGGGGGGACCGGCGCCGGGGGAGGAGAGTTTCACCCAACaaa                                                                                                                                                                                                                             >  2:502856/1‑239 (MQ=255)
         tCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGGCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCaaaa                                                                                                                                                                                                                            <  1:1080026/238‑1 (MQ=255)
             ttATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGGTGGCGGGGCTTGCCGGGACACTGGCGGGGAC‑‑CCGGCCGGGGGGGGTCCGGTTTACCCCGCCAAATggaa                                                                                                                                                                                                                       >  1:418900/1‑236 (MQ=255)
             ttATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGGAGGCGGCGCTGACCGGTACACTGGCGGGGAA‑‑ACGGCCGGGGGCGATCAGtttttccccacacaagaaac                                                                                                                                                                                                                       >  1:745486/1‑222 (MQ=255)
                     aTAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGGAC‑‑GCGTCCGGGGGGGATTCGGTTCACCCACCAAAATGGAACGGAAGCg                                                                                                                                                                                                               >  1:288268/1‑239 (MQ=255)
                     aTAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGGACACTGGCGGGGAT‑‑CGGTCCGGGGGAGGTACGTTTTTCCCCCCAAAAAGAACAGGaaagg                                                                                                                                                                                                               >  1:742619/1‑236 (MQ=255)
                       aaCGAATGTCAACTGGCTTCAGTTAGGGCGCACGAAATCTCCCACGTCCCCCAACGTCACCTGGCGCGAGCGATGGAATATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGa                                                                                                                                                                                                             <  1:986265/239‑1 (MQ=255)
                        aCGAAAGTCAAATGGCTTCAGTTAAGGCGCACGAAACCTCCCACGTCACCCAACGGCACCTGGCGCGAGCGGTGGAAGATCAGCAGCGCAGCGCGCCGCTGGCCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGAc                                                                                                                                                                                                            <  1:689163/239‑1 (MQ=255)
                          gAAAGTCAACTGGGTTCAGTTAGGGCGCACGAAATCTCCCCCGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCg                                                                                                                                                                                                          <  2:894990/239‑1 (MQ=255)
                             aGTCTACTTGCTTCAGTTATGGCCCACTAAATCTCCGCCGTCACCCAACGTCCCCGGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGGCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCAt                                                                                                                                                                                                      <  2:453374/239‑1 (MQ=255)
                             aGTCAAGTGGCTTCAGTTAGGGCGCACGAAATCTCCCACGTCCCCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCAt                                                                                                                                                                                                      <  1:593910/239‑1 (MQ=255)
                             aGTCAACTGGCTTTAGTTAGGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCAt                                                                                                                                                                                                      <  1:447067/239‑1 (MQ=255)
                             aGTCAACTGGCTTCAGTTATGGCGCACGAAATCCCCCACGTCCCCCAACGTCCCCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGGCCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCAt                                                                                                                                                                                                      <  1:502856/239‑1 (MQ=255)
                             aGTCAAATGGCTTCAGTTATGGCACACGAAATCTCCCACGGCCCCCAACGTCACCGGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCAt                                                                                                                                                                                                      <  1:1066497/239‑1 (MQ=255)
                                        ttCAGTTTTGGCGCACGAAATCCCCCACCGCCCCCAACGGCACCTGGCGGGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGGGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAgg                                                                                                                                                                                           <  2:1335090/239‑1 (MQ=255)
                                               aTGGCGCCAGAAATCCCCCACGGCCCCCAACCTCGCCGGGCGCGAGGGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCa                                                                                                                                                                                    <  2:745486/239‑1 (MQ=255)
                                                    gcACGAAACTTCCCACGTCCCCCAACGTCACCGGGGGGAGGCGATGGAAGATCAGCAGGGCGGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCt                                                                                                                                                                               <  2:418900/239‑1 (MQ=255)
                                                     cacacaaTCTCCCACGGCCACCAACCTCACCCGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCGGACCTGGGGCGGCGCGTTAGGTTCTATTTGACTGGAGATGGCCAGTCCGCAGGCGGGGATGGCGGAGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTc                                                                                                                                                                              <  1:1082467/234‑1 (MQ=255)
                                                     cACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGGAA‑‑GCGTGAGGGGGTGACCAGATTTTCCCAGCCAAAACAAAAGGGAACGGGCCGCCTTGGTGGTTCTGTGCTGCaaaaccc                                                                                                                                                                              >  1:552431/1‑234 (MQ=255)
                                                           aTCCCCCACGTCCCCCAACGTCCACGGGCGCGAGCGATGGAAGATCAGCAGCGTAGCGCGCCGCTGCCCGGGGTCGGCGCGTTAGGTTCTATTTTACGGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGAtt                                                                                                                                                                        <  1:477732/239‑1 (MQ=255)
                                                                       aCCCAACGTCACCGGGCGCGAGCGAGTGAAGATCAGCAGCGCAGCGCGCCGCTGCCCTGGGTCGGCGCGTTAGGTTCTATTTTACGGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAgg                                                                                                                                                             <  1:583076/238‑1 (MQ=255)
                                                                                           ggcggTGGAAGGTCAGCAGGGCAGCGGCGCGTTGACCGGGGGCGGCGGGTTAGGTTCTATTTGACTGGGGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGaa                                                                                                                                         <  1:579374/237‑1 (MQ=255)
                                                                                                   aaGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACa                                                                                                                                                                                                                      <  1:1345349/154‑1 (MQ=255)
                                                                                                       tCAGCAGCGCAGCGCGGCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATc                                                                                                                            <  1:1159624/239‑1 (MQ=255)
                                                                                                            cgcgcaggcgggcgttGACCTGGGTCGGCCGGTTAGGTTCTATTTTACTGGCTATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTTCACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCCCCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGgc                                                                                                                        <  1:329798/223‑1 (MQ=255)
                                                                                                                                           ctttgttgtttttttCCGGGCGTTGGCCGGTCCGCAGGCGGGGATGGAGGCGGTGAGCGGGACGCTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTa                                                                                        <  2:515588/227‑1 (MQ=255)
                                                                                                                                                                gatggcacatCCGCAGGCGGGGATGGCGGCGCTGACCGTTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCCCCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGAGTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGaa                                                                   <  1:442172/229‑1 (MQ=255)
                                                                                                                                                                gATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGaa                                                                   <  1:1206694/238‑1 (MQ=255)
                                                                                                                                                                gATGGCCAGTCCGCAGGCGGGGATGGCGGCGCCGACCGGTAAACTGGCGGGGAAC‑‑CCGCCGGGGGGGGATCGGTTTACCCCCGCAAAAAAAAAGGAAACGGGCCGCAATTGGTTTTCGTGGCCTCCACGCCCCGGGGTCCGATCGCCGGCGGAGCCAACCCTCCTCCGAAAAAAAATCGGTCAGGAGCGGTGATACTCCCCCCCGCCCGGAAAttttttatatccccccctgcccggaa                                                                   >  2:495670/1‑217 (MQ=255)
                                                                                                                                                                      cAGTCCGCGGGCGGGGATGGCGGCGCTGACCGGTACACGGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAgtcgt                                                             <  1:226762/238‑1 (MQ=255)
                                                                                                                                                                            ccggggggggATGTCGGGGCTGACCGGGACACTGGCGGGAAC‑‑GGGTCAGGGGATGAACAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAg                                                      <  2:288268/233‑1 (MQ=255)
                                                                                                                                                                                       tGGGGGCGCTGACCGGTACACTGGCGGGAAT‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCaa                                            <  1:1203503/238‑1 (MQ=255)
                                                                                                                                                                                       gtatgtgtataagagacagacgggCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATGCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCaa                                            <  1:279116/216‑1 (MQ=255)
                                                                                                                                                                                          cagcgcTGGACGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCg                                         <  1:742705/236‑1 (MQ=255)
                                                                                                                                                                                           cggcgtaaccgggtccatGGCGGGAAC‑‑GCGTCAGGGGGAGAACAGTTTCACCCAGCAAAATGAACAGGAAGGGGAACGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGAAGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCATCGCGCCCGCCGGAAATTTGATTTACTCACCCGTGGCCGGAAAGTCGTCTGGCAGATGCCCCCAACCGTg                                       <  2:1241711/222‑1 (MQ=255)
                                                                                                                                                                                                     ggTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGAt                             <  1:495670/239‑1 (MQ=255)
                                                                                                                                                                                                     ggTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGCGCTAATCAGAt                             <  1:316385/239‑1 (MQ=255)
                                                                                                                                                                                                       tACACTGGCGGTAAC‑‑CCGTCAGGGGAAGATACGTTTTACCCCGCCAAATTAACAGGAAGCGGACCGCATTGGTATTCAGGTACTGCAACGCTCGGGATTCGATCCCCAGGCGACGCCAACCTTACTCGCAAAATTACTCGCTCAGGCGCGCTATTCCCCACGCCCCCCCGGCACTTTTTTTACACTACCCCTGCCCGAAAGTCGCCTGGCAGAGGcccccgcacaggctaaagcaaagg                           >  2:94178/1‑218 (MQ=255)
                                                                                                                                                                                                       tACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTGCTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATg                            <  1:979385/238‑1 (MQ=255)
                                                                                                                                                                                                       gtcacaggggggggAC‑‑GTGTCAGGGGGTGTTCAGTTTCCCCCCACAAAACGAACCGGGAGCGGACCGACTTGGTATTCAGGTGCTGCAACGCTCGGGATTCGACCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCTGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTAGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATg                            <  2:545980/229‑1 (MQ=255)
                                                                                                                                                                                                        ggggaagcggggggacaggtgcAGGGGATGGCCATTTTCCCCCACCAAAAGTAACAGGAAGCGCACCGCATTGGTATTCAGGTGCTGCACCGCTCGGGATTCGATCCGCAGGCCATGCCAACCTTCCTCGAAAAATTACTCGATCATGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCTCCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGAt                             <  2:251746/218‑1 (MQ=255)
                                                                                                                                                                                                           cagggggggAAC‑‑GGGTCGGGGGATGATCAGTTTCCCCCACAAAAATGAACAGGAAGTGGACCGCTTTTGTATCTAGGTGCTCCAACGCTCGGGCGTCGATCCCCAGGCGATGCCACACTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCTCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGcc                        <  2:1363961/233‑1 (MQ=255)
                                                                                                                                                                                                              gCGGGAAC‑‑GCGTCAGGGGGTGATCAGTTTCACCCAGCAAAAAGAACAGGAAGCGGACCGCATTGGTATTCAGGGGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAAACTTCCTCGAAAAATTACTCGATCAGGCGCGTTCCTCCTCCCCCCCCCCGGAAATTTTTTTGAATCACCCGTTGCCGGAAAAGCGTCTGGCCGAGGCCCCCAACCgcggcaaccagcatggcccgac                    >  2:1253528/1‑218 (MQ=255)
                                                                                                                                                                                                                        cggctaaGGGGATTATCAGTTTCACACAACAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGgtggtg             <  1:1079866/233‑1 (MQ=255)
                                                                                                                                                                                                                           gTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATTCGCCCCCTggtggt              >  1:1363961/1‑235 (MQ=255)
                                                                                                                                                                                                                                 ggATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAgtcgt      <  1:130092/237‑1 (MQ=255)
                                                                                                                                                                                                                                   atgatCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGCGCTGCAACTCTCGGGATTCGACACCCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCCCGTTACTCCTCACGCCCGCCGGAAATTTTATTGACTCACCAGTTGCCGGAAAGTCCCCCGGAAGATCCCCGCTAACGAGCTCATATGACGAGCCCgccgggggggcggcccgcgga  >  2:180359/1‑219 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
CTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAC‑‑GCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGCGATGCCAACCTTCCTCGAAAAATTACTCGATCAGGCGCGTTACTCCTCGCGCCCGCCGGAAATTTTATTGACTCACCCGTTGCCGGAAAGTCGTCTGGCAGATGCCCGCAACCGTGCTAATCAGATGCGCCCGATGGTGGTGCAGTCGTCGGA  >  NC_000913/2616444‑2616906

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