Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,846,200 T→G 6.2% N323T (AAC→ACC)  uhpT ← hexose‑6‑phosphate:phosphate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,846,2000TG 6.2% 232.9 / 7.0 81N323T (AAC→ACC) uhpThexose‑6‑phosphate:phosphate antiporter
Reads supporting (aligned to +/- strand):  ref base T (34/42);  new base G (2/3);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

AATCAACAATTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCACTGGTCGATACCAATACGTACC  >  NC_000913/3846055‑3846349
                                                                                                                                                 |                                                                                                                                                     
aaTCAACAATTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGc                                                                                                                                                   >  1:1286446/1‑150 (MQ=255)
 aTCAACAATTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGCCCGTTCGcc                                                                                                                                                  >  2:384922/1‑150 (MQ=255)
 aTCAACAATTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGcc                                                                                                                                                  <  2:559941/150‑1 (MQ=255)
 aTCAACAATTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAACCCAGATAGATATATTCGTTTCTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGcc                                                                                                                                                  >  2:379125/1‑150 (MQ=255)
   caacaaTTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAg                                                                                                                                                >  2:1159347/1‑150 (MQ=255)
     acaaTTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGt                                                                                                                                              <  2:402371/150‑1 (MQ=255)
      caaTTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTc                                                                                                                                             >  1:1187522/1‑150 (MQ=255)
      caaTTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTc                                                                                                                                             >  2:203568/1‑150 (MQ=255)
      caaTTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTc                                                                                                                                             <  2:393291/150‑1 (MQ=255)
      caaTTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGAACGGTCGcggtggc                                                                                                                                             <  2:1258556/150‑7 (MQ=255)
      caaTTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGAACGGTCGccgtggc                                                                                                                                             <  2:1403274/150‑6 (MQ=255)
      caaTTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGAACGGTCGccgtggc                                                                                                                                             <  2:1180975/150‑6 (MQ=255)
         ttGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGGGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCaga                                                                                                                                          <  2:228709/150‑1 (MQ=255)
          tGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGCCCGTTCGCCAGGTCagag                                                                                                                                         >  2:1374173/1‑150 (MQ=255)
            cGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCagagag                                                                                                                                       <  2:1482289/150‑1 (MQ=255)
            cGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCagagag                                                                                                                                       <  1:805210/150‑1 (MQ=255)
              ggCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGagcc                                                                                                                                     >  2:1501527/1‑150 (MQ=255)
              ggCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGagcc                                                                                                                                     <  1:944294/150‑1 (MQ=255)
                    aGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGcccc                                                                                                                               <  2:531013/150‑1 (MQ=255)
                    aGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGcccc                                                                                                                               >  2:1307011/1‑150 (MQ=255)
                    aGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGcccc                                                                                                                               >  2:695932/1‑150 (MQ=255)
                       ccAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATATTCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCcaca                                                                                                                            <  2:540915/150‑1 (MQ=255)
                       ccAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCcaca                                                                                                                            >  1:273975/1‑150 (MQ=255)
                         aGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCAcagc                                                                                                                          >  2:330079/1‑150 (MQ=255)
                               cccAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTa                                                                                                                    >  2:337939/1‑150 (MQ=255)
                               cccAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTa                                                                                                                    >  2:326317/1‑150 (MQ=255)
                                ccAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACATCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTAc                                                                                                                   >  2:1032956/1‑150 (MQ=255)
                                     gCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGAcc                                                                                                              <  1:21413/150‑1 (MQ=255)
                                      cAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCa                                                                                                             >  1:1501417/1‑150 (MQ=255)
                                         agagagAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAgcg                                                                                                          <  1:203568/150‑1 (MQ=255)
                                         agagagAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAgcg                                                                                                          <  2:1424971/150‑1 (MQ=255)
                                         agagagAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGGGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAgcg                                                                                                          <  1:1607480/150‑1 (MQ=255)
                                          gagagaAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAtcgc                                                                                                         <  2:190003/150‑1 (MQ=255)
                                          gagagaAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATCCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAgcgc                                                                                                         <  1:40977/150‑1 (MQ=255)
                                                  ccAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCtt                                                                                                 <  2:591503/150‑1 (MQ=255)
                                                         gataTATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAg                                                                                          <  2:36390/150‑1 (MQ=255)
                                                              atTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTaa                                                                                     >  2:527065/1‑150 (MQ=255)
                                                              atTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTaa                                                                                     >  1:1609394/1‑150 (MQ=255)
                                                              atTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTaa                                                                                     >  2:1490854/1‑150 (MQ=255)
                                                              atTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTaa                                                                                     >  2:1034447/1‑150 (MQ=255)
                                                                tCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGCCCAGCGCACCAGCTTCAAACAGCGTAAAg                                                                                   >  1:1072261/1‑150 (MQ=255)
                                                                    tACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCca                                                                                                                              <  2:692648/103‑1 (MQ=255)
                                                                    tACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCca                                                                                                                              >  1:692648/1‑103 (MQ=255)
                                                                     aCTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTg                                                                              <  1:330079/150‑1 (MQ=255)
                                                                      cTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGa                                                                             <  2:747030/150‑1 (MQ=255)
                                                                       ttgCATGTTGATACACACCGAGCGTGGCTATAATCAGCGCCAGCGCGATGCAGGCCACCATGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGaa                                                                            <  1:765621/148‑1 (MQ=255)
                                                                         gCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACGAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATc                                                                          <  2:574536/150‑1 (MQ=255)
                                                                             gTTGATACACACCGAGCGTGGAGATAATCAGCGCCAGCGCGATGCAGGACACCAGGCCACGGAGCCCGTTCGCCAGGTCAGAGAGCCAGCACCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCTTAAAGCCCTGAATCGCCa                                                                      >  1:1081784/1‑150 (MQ=255)
                                                                                  tACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCt                                                                 <  2:1286446/150‑1 (MQ=255)
                                                                                  tACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCt                                                                 >  2:1300796/1‑150 (MQ=255)
                                                                                    cacacCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGCCCAGCGCACCAGCTTCAAACCGCGTAAAGCCCTGAATCGCCACCGCttt                                                               >  1:286095/1‑150 (MQ=255)
                                                                                      cacCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCCGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTag                                                             >  2:340395/1‑150 (MQ=255)
                                                                                      cacCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTag                                                             >  2:446507/1‑150 (MQ=255)
                                                                                              gTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACGGTTCGCCATGTCTGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGATTCGCCACCGCTTTAGAGAGTTTc                                                     <  2:1250895/150‑1 (MQ=255)
                                                                                              gTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTc                                                     <  1:1490854/150‑1 (MQ=255)
                                                                                              gTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTc                                                     <  1:1034447/150‑1 (MQ=255)
                                                                                                ggCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAg                                                   <  2:917162/150‑1 (MQ=255)
                                                                                                 gCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGt                                                  >  1:471301/1‑150 (MQ=255)
                                                                                                      aaTCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTg                                             <  1:1127738/150‑1 (MQ=255)
                                                                                                      aaTCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTg                                             <  1:911388/150‑1 (MQ=255)
                                                                                                      aaTCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTg                                             <  1:1307011/150‑1 (MQ=255)
                                                                                                       aTCCGCGGCAGCGCGATGCAGGCCACGAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTgg                                            <  2:1044132/150‑1 (MQ=255)
                                                                                                       aTCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTgg                                            >  2:1200038/1‑150 (MQ=255)
                                                                                                         cagcgccagcgcGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGaa                                          >  2:277388/1‑150 (MQ=255)
                                                                                                            cgccagcgcGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAAcgc                                       >  1:1264494/1‑150 (MQ=255)
                                                                                                             gccagcgcGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAAcgcg                                      <  2:273975/150‑1 (MQ=255)
                                                                                                             gccagcgcGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAAcgcg                                      >  1:1413837/1‑150 (MQ=255)
                                                                                                                agcgcGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGtat                                   >  1:931061/1‑150 (MQ=255)
                                                                                                                    cgATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACgg                               <  2:1187522/150‑1 (MQ=255)
                                                                                                                         caggccaccaggccacGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGCCCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGAc                          >  1:1546021/1‑150 (MQ=255)
                                                                                                                          aggccaccaggccacGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGAcc                         <  1:1159347/150‑1 (MQ=255)
                                                                                                                           ggccaccaggccacGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCa                        <  2:1264494/150‑1 (MQ=255)
                                                                                                                                    gccacGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCACTGGTCGAt               <  2:1501417/150‑1 (MQ=255)
                                                                                                                                    gacctggCGAACGGTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCACTGGTCGAt               >  2:1183017/6‑150 (MQ=255)
                                                                                                                                      cctggCGAACGGTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGCCCACTGGTCGATAc             >  2:1123871/4‑150 (MQ=255)
                                                                                                                                        cGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCACTGGTCGATACCa           <  2:1089017/150‑1 (MQ=255)
                                                                                                                                         ggCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCACTGGTCGATACCaa          <  2:471301/150‑1 (MQ=255)
                                                                                                                                          gCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCACTGGTCGATACCAAt         <  1:1300796/150‑1 (MQ=255)
                                                                                                                                             aCCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTTTACGGTTGACCACTGGTCGATACCAATACg      <  1:1032956/150‑1 (MQ=255)
                                                                                                                                              ccGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCATCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCACTGGTCGATACCAATACGt     <  1:528540/150‑1 (MQ=255)
                                                                                                                                              ccGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCACTGGTCGATACCAATACGt     <  1:337939/150‑1 (MQ=255)
                                                                                                                                                 ttCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGTCCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCACTGGTCGATACCAATACGTAcc  <  2:931061/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AATCAACAATTGCGGGCCAAAGACCAGGAAACCCAACGCAAAGAGAGAAGCCAGATAGATATATTCGTTACTGGCATGTTGATACACACCGAGCGTGGCGATAATCAGCGCCAGCGCGATGCAGGCCACCAGGCCACGGCGACCGTTCGCCAGGTCAGAGAGCCAGCCCCACAGCAGCGTACCGACCAGCGCACCAGCTTCAAACAGCGTAAAGCCCTGAATCGCCACCGCTTTAGAGAGTTTCAGTTCCTGGAACGCGTATACGGTTGACCACTGGTCGATACCAATACGTACC  >  NC_000913/3846055‑3846349

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