Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,491,690 A→G 6.0% intergenic (‑70/+35) yhfG ← / ← ppiA DUF2559 domain‑containing protein YhfG/peptidyl‑prolyl cis‑trans isomerase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,491,6900AG 6.0% 200.9 / 2.8 66intergenic (‑70/+35)yhfG/ppiADUF2559 domain‑containing protein YhfG/peptidyl‑prolyl cis‑trans isomerase A
Reads supporting (aligned to +/- strand):  ref base A (21/41);  new base G (2/2);  total (23/43)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGACGGCGACTGGCCTGAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTC  >  NC_000913/3491541‑3491839
                                                                                                                                                     |                                                                                                                                                     
aGACGGCGACTGGCCTGAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCa                                                                                                                                                       <  2:1408441/150‑1 (MQ=255)
         cTGGCCTGAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCAGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCgg                                                                                                                                              >  2:379008/1‑150 (MQ=255)
          tGGCCTGAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCggg                                                                                                                                             <  2:1300586/150‑1 (MQ=255)
             ccTGAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCaa                                                                                                                                          >  2:1490721/1‑150 (MQ=255)
                    attgcgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGcc                                                                                                                                   >  1:1631811/1‑150 (MQ=255)
                     ttacgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGccc                                                                                                                                  <  1:978665/150‑1 (MQ=255)
                     ttacgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGccc                                                                                                                                  <  1:508164/150‑1 (MQ=255)
                        cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCcgcg                                                                                                                               <  1:757443/150‑1 (MQ=255)
                          attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcgct                                                                                                                             <  1:1385575/150‑2 (MQ=255)
                          attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcgct                                                                                                                             <  2:1229005/150‑2 (MQ=255)
                          attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcgct                                                                                                                             <  1:748857/150‑2 (MQ=255)
                             acgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCgaga                                                                                                                          <  2:938752/150‑1 (MQ=255)
                              cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGTAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCCCGCCAGaa                                                                                                                         >  2:730159/1‑150 (MQ=255)
                              cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGCTTGCCCGCGCgagca                                                                                                                         >  1:1597643/1‑148 (MQ=255)
                              cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGaa                                                                                                                         >  2:1166200/1‑150 (MQ=255)
                               gCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGaaa                                                                                                                        <  2:1418299/150‑1 (MQ=255)
                                 tGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCgag                                                                                                                           >  1:38437/1‑145 (MQ=255)
                                 tGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCgag                                                                                                                           <  2:38437/145‑1 (MQ=255)
                                 ggcaGCTCCCCGCGCCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATc                                                                                                                      <  2:394224/147‑1 (MQ=255)
                                   aaGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                    <  2:1373759/150‑1 (MQ=255)
                                      ctctCAGAGACGGGACTTTTGCTTATCGGTGAGGTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTAc                                                                                                                 <  1:1034161/150‑1 (MQ=255)
                                          cAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCa                                                                                                             <  2:385035/150‑1 (MQ=255)
                                          cAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCa                                                                                                             <  1:824950/150‑1 (MQ=255)
                                            cagaCGGGACTTTTGCTTTTCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCTCAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAgg                                                                                                           <  1:1271332/149‑1 (MQ=255)
                                              gaCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCAGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGAc                                                                                                         <  1:379008/150‑1 (MQ=255)
                                               aCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACt                                                                                                        <  1:1899842/150‑1 (MQ=255)
                                                 gggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACttt                                                                                                      <  1:588718/150‑1 (MQ=255)
                                                  ggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACtttt                                                                                                     >  2:132035/1‑149 (MQ=255)
                                                     cTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCg                                                                                                  <  1:1611012/150‑1 (MQ=255)
                                                     cTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCg                                                                                                  >  1:671199/1‑150 (MQ=255)
                                                      ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCgg                                                                                                 <  2:1780966/150‑1 (MQ=255)
                                                        ttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCTGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGaa                                                                                               <  1:904866/150‑1 (MQ=255)
                                                         tGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGaaa                                                                                              <  1:235891/150‑1 (MQ=255)
                                                          gCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAg                                                                                             <  1:1065449/150‑1 (MQ=255)
                                                           cgatctggTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAgg                                                                                            <  1:1039979/144‑1 (MQ=255)
                                                             tATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGAt                                                                                          >  1:1478602/1‑150 (MQ=255)
                                                               tCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATaa                                                                                        <  2:1388404/150‑1 (MQ=255)
                                                                cGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATATTGCCCGCGCGAGATATGATTAAGGGAGGAGTTTAGGGGAAATGTTaaa                                                                                       >  1:590290/1‑149 (MQ=255)
                                                                cGGTGAGTTTCTTCACTTGATGCCACCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAAc                                                                                       <  1:1980634/150‑1 (MQ=255)
                                                                 ggTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCCGGCAATAATTCCCCCGCGAGAAAAAATTAAGGCAGGACGTTAGCGGAAAGGagatag                                                                                      >  2:1039979/1‑145 (MQ=255)
                                                                    gAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCGGAAGCGCGGGCAAGATTTCCCGCGCGAGAAATAATTAAGGCAGGACTTTTGCGGAAAGGATAACTaac                                                                                   >  2:413385/1‑148 (MQ=255)
                                                                    gAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCGGAAGCGCGGGCAAGATTGCCCCCGCGAGAAATAAATACGGCAGGACTTTAGAGGAAAGGATAACAAcc                                                                                   >  2:406092/1‑150 (MQ=255)
                                                                      gTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGGAAGAATGCCCGCGCGAGAAATCATTACGGCAGGACATTAGCGGAAAGGATAAATAcccg                                                                                 >  1:342844/1‑148 (MQ=255)
                                                                          cttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCgg                                                                                 >  1:161846/1‑146 (MQ=255)
                                                                          cttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCgg                                                                                 <  2:161846/146‑1 (MQ=255)
                                                                           ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTg                                                                            <  2:1478602/150‑1 (MQ=255)
                                                                                 tGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGca                                                                      <  1:1779709/150‑1 (MQ=255)
                                                                                  gATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGcac                                                                     <  2:1803360/150‑1 (MQ=255)
                                                                                    tGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACAt                                                                   >  2:1356027/1‑150 (MQ=255)
                                                                                              aaCGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACgg                                                         <  2:671199/150‑1 (MQ=255)
                                                                                                    tCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGGAATCATTACGGCAGGGCGTTTGCGGGAAGGATAACTACCGGTTTTGACGGCCCATTCTGGTACGGGCCAAc                                                   >  1:1441129/1‑150 (MQ=255)
                                                                                                        tttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCa                                               <  1:1490721/150‑1 (MQ=255)
                                                                                                          tGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCGGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                             <  1:406092/150‑1 (MQ=255)
                                                                                                          tGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCGGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                             <  1:413385/150‑1 (MQ=255)
                                                                                                          tGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                             <  2:752394/150‑1 (MQ=255)
                                                                                                                      aaGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCCAGAAATAATTACGGCAGGAAATTAGCGGAAAGGATAACAACCCGTTTTGACGGCACATTATGGTACGGACCAACGTCAAGAGTGGGCAccgg                                 >  2:131777/1‑148 (MQ=255)
                                                                                                                       aGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAATAAGGCACGACTTTAGCGGAAAGGATAACAACCCGTTTTGACGGCACAATATGGTACCGACCAACCTCAAGACTAGGCACCggg                                >  2:299601/1‑150 (MQ=255)
                                                                                                                           gCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGaaa                            <  2:528494/150‑1 (MQ=255)
                                                                                                                              aaTCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                         <  1:1923497/150‑1 (MQ=255)
                                                                                                                                     gTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGc                  <  2:670356/150‑1 (MQ=255)
                                                                                                                                         cGCCGGGAGAGCAGAAGCGCGGGCAAGATTTCCCGCGCGAGAAATAATTACGGCAGGACTTTAGGGGAAAGGATAACAACCGGTTTTTACGGCACATTATGGTACGGACCAACGTCAAGAGTCGGGACCTGGGAAAACTTATAGGCaaag              >  2:382816/1‑148 (MQ=255)
                                                                                                                                               gagaGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATg        <  1:299601/150‑1 (MQ=255)
                                                                                                                                                 tagCAGAAGCGCTGGAAAGATTTCCCGCGAGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGcc      <  2:286499/149‑1 (MQ=255)
                                                                                                                                                 tagCAGAAGCGCGGGCAAGATTGCCCGCGCTAGAAATCATTACGGCAGTACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGTAAATCTTATCGGCAACGTCCATGcc      <  2:279088/149‑1 (MQ=255)
                                                                                                                                                 gagCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGcc      <  2:61532/150‑1 (MQ=255)
                                                                                                                                                  agCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATAAATACGGCAGGCCTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGGCCAAAGTCATTAGTAGGCACCGGGGAAATCTTATCGGGAACGTCCATGccg     >  2:1939887/1‑149 (MQ=255)
                                                                                                                                                    cAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCttt   >  1:864872/1‑150 (MQ=255)
                                                                                                                                                     aGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc  >  2:1703215/1‑150 (MQ=255)
                                                                                                                                                     aGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc  <  2:1441129/150‑1 (MQ=255)
                                                                                                                                                     aGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTAAGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACCTCATGAGTCGGCACCGGGGAAATCTTATAGGCAACGTCCATGCCTTTc  >  1:199138/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGACGGCGACTGGCCTGAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTC  >  NC_000913/3491541‑3491839

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