Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,081,057 A→G 10.7% V69A (GTG→GCG) ‡ sbmC ← DNA gyrase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,081,0570AG10.7% 252.0 / 18.4 95V69A (GTG→GCG) ‡sbmCDNA gyrase inhibitor
Reads supporting (aligned to +/- strand):  ref base A (42/42);  new base G (2/8);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 9.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

GCAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCC  >  NC_000913/2080908‑2081197
                                                                                                                                                     |                                                                                                                                            
gCAAGAGACTATTAAAGAACTGATACTATGGTTTAGCAAAATCATCAACTACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCa                                                                                                                                              <  2:791020/150‑1 (MQ=255)
   agagACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCacga                                                                                                                                           >  1:1499083/1‑150 (MQ=255)
   agagACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCacga                                                                                                                                           >  1:1461602/1‑150 (MQ=255)
      gaCTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACgg                                                                                                                                        <  2:355324/150‑1 (MQ=255)
       aCTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGgt                                                                                                                                       <  2:1587555/150‑1 (MQ=255)
       aCTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGACACCGTCGTGACGgt                                                                                                                                       <  1:353613/150‑1 (MQ=255)
        cTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGgtg                                                                                                                                      >  1:1333032/1‑150 (MQ=255)
            tAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGc                                                                                                                                  >  1:882757/1‑150 (MQ=255)
            tAAAGAACTGATACCAGGCTAAAGCAAAATCATCACCGACTACCCGAGCTACCCCCCCCGCATACTGACCACCTGAAAATTCACACAGAATCACGCCCTCACTGTTTTCGGCAAGCGTAAAGTAACCCGCCACCCTCACCACggactccc                                                                                                                                  >  2:1246600/1‑144 (MQ=255)
             aaaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCGTGACGGTGCCggg                                                                                                                                 <  1:1346409/150‑3 (MQ=255)
             aaaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCGTGACGGTGCCggg                                                                                                                                 <  2:722153/150‑3 (MQ=255)
             aaaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGACACCGTCGTGACGGTGCCggg                                                                                                                                 <  1:1354290/150‑3 (MQ=255)
              aaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAg                                                                                                                                <  2:1564097/150‑1 (MQ=255)
              aaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGACACCGTCGTGACGGTGCCgggt                                                                                                                                <  1:729640/150‑4 (MQ=255)
              aaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGACACCGTCGTGACGGTGCCgggt                                                                                                                                <  1:1529249/150‑4 (MQ=255)
               aGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGATGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGc                                                                                                                               >  1:369127/1‑150 (MQ=255)
               aGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGACACCGTCGTGACGGTGCCgggtt                                                                                                                               <  1:969926/150‑5 (MQ=255)
               aGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGACACCGTCGTGACGGTGCCgggtt                                                                                                                               <  2:1423976/150‑5 (MQ=255)
                   cTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGCCCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTaa                                                                                                                           >  2:300029/1‑150 (MQ=255)
                    tGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAAt                                                                                                                          >  1:352072/1‑150 (MQ=255)
                     gATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAAtt                                                                                                                         >  1:1575794/1‑150 (MQ=255)
                       tACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTGACCCGGCACCGTCACGACGGTGTCGCAGCGTAAtttt                                                                                                                       <  1:995075/150‑1 (MQ=255)
                       tACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAAtttt                                                                                                                       >  1:1104692/1‑150 (MQ=255)
                         ccAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTc                                                                                                                     >  1:949142/1‑150 (MQ=255)
                         ccAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTACGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCTGCGTAATTTTTc                                                                                                                     >  2:914768/1‑150 (MQ=255)
                             ggTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcggcg                                                                                                                 >  1:487723/1‑150 (MQ=255)
                             ggTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcggcg                                                                                                                 <  2:564852/150‑1 (MQ=255)
                              gTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGATGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcggcgg                                                                                                                <  2:369127/150‑1 (MQ=255)
                                 tAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGgtg                                                                                                             <  1:1598331/150‑1 (MQ=255)
                                    cAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGCCGGTGTCGCAGCGTAATTTTTCGGCGGGTGttt                                                                                                          >  1:209121/1‑150 (MQ=255)
                                      aaaTCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTTTTTCa                                                                                                        <  1:1123155/150‑1 (MQ=255)
                                        atcatcACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATc                                                                                                      >  1:123980/1‑150 (MQ=255)
                                          catcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTg                                                                                                    <  2:949142/150‑1 (MQ=255)
                                           atcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTgg                                                                                                   >  1:1089404/1‑150 (MQ=255)
                                           atcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTgg                                                                                                   >  1:1271625/1‑150 (MQ=255)
                                           atcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTgg                                                                                                   >  1:1104205/1‑150 (MQ=255)
                                           atcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTgg                                                                                                   <  1:1312015/150‑1 (MQ=255)
                                              aCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGAtt                                                                                                <  2:487723/150‑1 (MQ=255)
                                               ccGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTg                                                                                               >  1:844350/1‑150 (MQ=255)
                                                  aCTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTca                                                                                            <  1:1465711/150‑2 (MQ=255)
                                                   cTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGt                                                                                           <  2:946303/150‑1 (MQ=255)
                                                    tACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTa                                                                                          <  1:1479412/150‑1 (MQ=255)
                                                     acacGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTaa                                                                                         >  2:474443/1‑150 (MQ=255)
                                                     acacGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTAATCTGGATTGTCGTaa                                                                                         >  2:861038/1‑150 (MQ=255)
                                                      cacGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAAt                                                                                        >  1:1341198/1‑150 (MQ=255)
                                                        cGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAg                                                                                      <  2:1104692/150‑1 (MQ=255)
                                                        cGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAg                                                                                      >  1:930666/1‑150 (MQ=255)
                                                        cGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAg                                                                                      >  2:50011/1‑150 (MQ=255)
                                                        cGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAg                                                                                      >  1:748621/1‑150 (MQ=255)
                                                          tgCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCAGTAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCATAATAGAc                                                                                    <  2:1548494/149‑1 (MQ=255)
                                                               ccGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGTATTGTCGTAATAGACAGCaa                                                                               >  2:1130798/1‑150 (MQ=255)
                                                                 gCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGCCGCCCTCACTGTTTTCGGGAAGCGGAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAAcc                                                                             >  1:435430/1‑150 (MQ=255)
                                                                 gCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAAcc                                                                             >  1:1304942/1‑150 (MQ=255)
                                                                 gCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAAcc                                                                             >  1:35387/1‑150 (MQ=255)
                                                                  ccACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAAccc                                                                            <  2:1461602/150‑1 (MQ=255)
                                                                  ccACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAAccc                                                                            >  2:265269/1‑150 (MQ=255)
                                                                     ccGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACt                                                                         <  2:215169/150‑1 (MQ=255)
                                                                     ccGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACt                                                                         >  1:780020/1‑150 (MQ=255)
                                                                       gCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCATAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTcc                                                                       >  2:315109/1‑150 (MQ=255)
                                                                        cATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCt                                                                      >  1:987640/1‑150 (MQ=255)
                                                                            cTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCgg                                                                  <  2:1333032/150‑1 (MQ=255)
                                                                                      gTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATAttt                                                        >  1:103961/1‑150 (MQ=255)
                                                                                        aaTTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATAttttt                                                      <  2:417062/150‑1 (MQ=255)
                                                                                        aaTTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATAttttt                                                      <  2:1499083/150‑1 (MQ=255)
                                                                                         aTTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGTCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTg                                                     <  2:200482/150‑1 (MQ=255)
                                                                                          tttCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGT‑AAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCt                                                   <  1:381306/150‑1 (MQ=255)
                                                                                           ttCTGTCAGAATGCCGCCCTCACTGTTTTGGGGAAGCGTAAAGTAACCCGGCACCGTCCCGACGGTGTCGCAGCGTAATGTTTCGGCGGGGGTTTCATCTGGATTGTCGTAATAGCCAGCACGCCCCTGCTTCGGCACAGTATTTTTggt                                                   >  2:1280496/1‑148 (MQ=255)
                                                                                              tGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATc                                                <  1:549288/150‑1 (MQ=255)
                                                                                                tCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTa                                              <  2:103961/150‑1 (MQ=255)
                                                                                                 cAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTAc                                             <  1:901370/150‑1 (MQ=255)
                                                                                                   gAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTAccc                                           >  1:539543/1‑150 (MQ=255)
                                                                                                    aaTGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCca                                          <  2:844350/150‑1 (MQ=255)
                                                                                                     aTGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCcac                                         <  1:1090296/150‑1 (MQ=255)
                                                                                                           cccTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACatcatc                                   <  1:474443/150‑1 (MQ=255)
                                                                                                            ccTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACAtcatca                                  >  2:734983/1‑150 (MQ=255)
                                                                                                               cACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACt                               >  2:1387817/1‑150 (MQ=255)
                                                                                                               cACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACt                               <  2:1341198/150‑1 (MQ=255)
                                                                                                                  tGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCt                            <  2:94672/150‑1 (MQ=255)
                                                                                                                   gTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTc                           >  1:983178/1‑150 (MQ=255)
                                                                                                                    ttttCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCa                          <  2:930666/150‑1 (MQ=255)
                                                                                                                    ttttCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCa                          <  2:987640/150‑1 (MQ=255)
                                                                                                                    ttttCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCa                          <  1:50011/150‑1 (MQ=255)
                                                                                                                        cGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGc                      <  1:1067317/150‑1 (MQ=255)
                                                                                                                        cGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGc                      <  1:365703/150‑1 (MQ=255)
                                                                                                                         gggAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGcc                     >  1:1077793/1‑150 (MQ=255)
                                                                                                                              gCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCctttct                >  1:378484/1‑150 (MQ=255)
                                                                                                                              gCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGATATCTACCCACATCATCAACTGCTCAAAGCctttct                >  1:422389/1‑150 (MQ=255)
                                                                                                                               cGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGATCAAAGCctttctt               >  2:1221813/1‑150 (MQ=255)
                                                                                                                                 tAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTa             <  2:1089404/150‑1 (MQ=255)
                                                                                                                                  aaaGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTAc            <  1:265269/150‑1 (MQ=255)
                                                                                                                                    aGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCg          <  2:435430/150‑1 (MQ=255)
                                                                                                                                     gTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGt         >  2:1372369/1‑150 (MQ=255)
                                                                                                                                      tAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTc        <  1:335374/150‑1 (MQ=255)
                                                                                                                                       gctgcgacACCGTCGTGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCt       >  1:369021/8‑150 (MQ=255)
                                                                                                                                       aaCCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCt       <  1:98543/150‑1 (MQ=255)
                                                                                                                                        ctgcgacACCGTCGTGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTg      >  1:1400241/7‑150 (MQ=255)
                                                                                                                                          ccGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGtt    >  1:1286296/1‑150 (MQ=255)
                                                                                                                                            ggCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTcc  <  1:1148964/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GCAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCC  >  NC_000913/2080908‑2081197

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