Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,081,058 C→T 10.4% V69M (GTG→ATG) ‡ sbmC ← DNA gyrase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,081,0580CT10.4% 270.6 / 14.8 95V69M (GTG→ATG) ‡sbmCDNA gyrase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (43/42);  new base T (2/8);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 9.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCC  >  NC_000913/2080911‑2081197
                                                                                                                                                   |                                                                                                                                           
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:1529249/150‑4 (MQ=255)
           aaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGACACCGTCGTGACGGTGCCgggt                                                                                                                                <  1:729640/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                                                                                                                 <  2:564852/150‑1 (MQ=255)
                          ggTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcggcg                                                                                                                 >  1:487723/1‑150 (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:1104205/1‑150 (MQ=255)
                                        atcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTgg                                                                                                   >  1:1271625/1‑150 (MQ=255)
                                        atcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTgg                                                                                                   >  1:1089404/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                                                                                      >  2:50011/1‑150 (MQ=255)
                                                     cGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAg                                                                                      >  1:748621/1‑150 (MQ=255)
                                                     cGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAg                                                                                      >  1:930666/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                                                                         >  1:780020/1‑150 (MQ=255)
                                                                  ccGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACt                                                                         <  2:215169/150‑1 (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:987640/150‑1 (MQ=255)
                                                                                                                 ttttCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCa                          <  2:930666/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)
                                                                                                                                                   |                                                                                                                                           
AGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCC  >  NC_000913/2080911‑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: ‑