Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,725,453 T→A 5.3% L98F (TTA→TTT)  kgtP ← alpha‑ketoglutarate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,725,4530TA 5.3% 242.7 / 2.9 75L98F (TTA→TTT) kgtPalpha‑ketoglutarate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (32/39);  new base A (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

TAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGG  >  NC_000913/2725304‑2725601
                                                                                                                                                     |                                                                                                                                                    
taGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGAt                                                                                                                                                      <  1:725547/150‑1 (MQ=255)
 aGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATa                                                                                                                                                     <  1:502844/150‑1 (MQ=255)
   gTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAAc                                                                                                                                                   >  2:964337/1‑150 (MQ=255)
    tGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACa                                                                                                                                                  >  1:416658/1‑150 (MQ=255)
       cgcTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCa                                                                                                                                               >  2:360254/1‑150 (MQ=255)
          tGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTGCATATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCg                                                                                                                                            <  1:389055/150‑1 (MQ=255)
             tGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAtt                                                                                                                                         >  2:1336160/1‑150 (MQ=255)
               ccATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAtttt                                                                                                                                       <  1:1024497/150‑1 (MQ=255)
                cATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAttttt                                                                                                                                      >  2:274891/1‑150 (MQ=255)
                cATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAttttt                                                                                                                                      >  2:320088/1‑150 (MQ=255)
                 atatTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAtttttt                                                                                                                                     <  2:852771/150‑1 (MQ=255)
                    ttCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCg                                                                                                                                  >  1:558728/1‑150 (MQ=255)
                     tctcCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGa                                                                                                                                 >  1:1064016/1‑150 (MQ=255)
                       tcCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGAcc                                                                                                                               >  1:106723/1‑150 (MQ=255)
                         cGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCAt                                                                                                                             >  1:907613/1‑150 (MQ=255)
                          gCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATg                                                                                                                            >  1:863473/1‑150 (MQ=255)
                           ccAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGt                                                                                                                           >  1:533505/1‑150 (MQ=255)
                            cAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGAt                                                                                                                                                      <  1:1353449/122‑1 (MQ=255)
                            cAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGAt                                                                                                                                                      >  2:1353449/1‑122 (MQ=255)
                            cAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGtt                                                                                                                          <  1:273176/150‑1 (MQ=255)
                            cAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGtt                                                                                                                          <  1:278255/150‑1 (MQ=255)
                                   tAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGc                                                                                                                   <  1:862104/150‑1 (MQ=255)
                                   tAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGc                                                                                                                   <  2:1197892/150‑1 (MQ=255)
                                    aaTCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCt                                                                                                                  <  2:863473/150‑1 (MQ=255)
                                      tCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTAt                                                                                                                <  2:1271090/150‑1 (MQ=255)
                                        ccTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGc                                                                                                              >  2:615783/1‑150 (MQ=255)
                                         cagAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGCCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCg                                                                                                             >  1:80534/3‑150 (MQ=255)
                                         cTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCg                                                                                                             >  1:1208921/1‑150 (MQ=255)
                                              aTAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCAGGGAGGCAGGCGAAAACCACCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATCTTTATCGGATATCCGCCCaa                                                                                                        >  1:113964/1‑150 (MQ=255)
                                               tAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCaaa                                                                                                       >  2:1251666/1‑150 (MQ=255)
                                               tAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGAAAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCaaa                                                                                                       >  2:540809/1‑150 (MQ=255)
                                               tAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGAAAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCaaa                                                                                                       >  2:535966/1‑150 (MQ=255)
                                                aaaCTAGCGAGAAGAAATAATGCCGTAGACCAAGTACCTATAGTTTAATAACCTGTGAGGCATGCGATAACCATAGATCCGAAAAACATCATACACACCGATAACAGCATCTATTTTTTGCTACCATTTTTATCTGCTATTCGGCCAAAt                                                                                                      <  2:398197/150‑1 (MQ=255)
                                                   cgagcgagAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGc                                                                                                   <  2:1102819/150‑1 (MQ=255)
                                                        gagaAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAAcc                                                                                              >  2:328150/1‑150 (MQ=255)
                                                         agaAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCg                                                                                             >  1:1292063/1‑150 (MQ=255)
                                                         agaAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACACCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCg                                                                                             <  1:964337/150‑1 (MQ=255)
                                                               aataatGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTAtt                                                                                       <  2:599128/150‑1 (MQ=255)
                                                                   atGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGgc                                                                                   >  1:1167312/1‑150 (MQ=255)
                                                                    tGCCGGAGCCCACGTACCTATAGCTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGgcg                                                                                  <  1:1054952/150‑1 (MQ=255)
                                                                    tGCCGGAGCCCACGTACCTATAGCTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGgcg                                                                                  <  1:1060009/150‑1 (MQ=255)
                                                                      ccGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCa                                                                                >  1:1217073/1‑150 (MQ=255)
                                                                           gCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAgg                                                                           <  1:1336160/150‑1 (MQ=255)
                                                                           gCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAgg                                                                           <  1:154142/150‑1 (MQ=255)
                                                                            cccACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGa                                                                          <  2:911971/150‑1 (MQ=255)
                                                                                     cTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCgcag                                                                 <  2:1080270/150‑1 (MQ=255)
                                                                                      tataGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCgcagc                                                                >  1:1087314/1‑150 (MQ=255)
                                                                                       ataGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGCCCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCgcagca                                                               >  1:952826/1‑150 (MQ=255)
                                                                                            ttCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAac                                                          >  1:1422185/1‑150 (MQ=255)
                                                                                             tCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAaca                                                         <  2:1064016/150‑1 (MQ=255)
                                                                                             tCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAaca                                                         >  2:271213/1‑150 (MQ=255)
                                                                                              cATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAacac                                                        >  2:612855/1‑150 (MQ=255)
                                                                                               aTAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACAcc                                                       <  1:615783/150‑1 (MQ=255)
                                                                                                    cTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACctgctg                                                  >  2:557459/1‑150 (MQ=255)
                                                                                                        gAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTg                                              <  1:940787/150‑1 (MQ=255)
                                                                                                         aggcaggcGATAACCAGCGATCCGAAACACATCATACACACCGAAAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTgt                                             <  1:535966/150‑1 (MQ=255)
                                                                                                         aggcaggcGATAACCAGCGATCCGAAACACATCATACACACCGAAAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTgt                                             <  1:540809/150‑1 (MQ=255)
                                                                                                          ggcaggcGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTgta                                            <  1:274891/150‑1 (MQ=255)
                                                                                                           gcaggcGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTgtag                                           <  2:558728/150‑1 (MQ=255)
                                                                                                            caggcGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTgtagt                                          >  1:429713/1‑150 (MQ=255)
                                                                                                            caggcGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTgtagt                                          <  2:1167312/150‑1 (MQ=255)
                                                                                                            caggcGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTgtagt                                          <  2:1327421/150‑1 (MQ=255)
                                                                                                              ggcGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGtagtag                                        <  2:254396/150‑1 (MQ=255)
                                                                                                              ggcGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGtagtag                                        <  1:316531/150‑1 (MQ=255)
                                                                                                               gcGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTagtagt                                       <  1:297255/150‑1 (MQ=255)
                                                                                                                 gATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTg                                     <  1:665627/150‑1 (MQ=255)
                                                                                                                 gATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTg                                     <  1:802531/150‑1 (MQ=255)
                                                                                                                 gATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTg                                     >  1:1027068/1‑150 (MQ=255)
                                                                                                                   tAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAg                                   <  2:801921/150‑1 (MQ=255)
                                                                                                                    aaCCAGCGATGCGAAACACATCATACACAGCGATTACAGGATCGATTTTTTGCGAGGATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATGCGGCAGCAATAACAACTGGTGTTTGTAGTTGTTGAgt                                  <  2:113964/150‑1 (MQ=255)
                                                                                                                      ccAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGAATCAGGAATCCCGCAGCAAAAACACCTGCTTTTTGTATTAGTTGAgtcg                                >  2:1419641/1‑150 (MQ=255)
                                                                                                                                 aaaCACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCg                     <  2:590948/150‑1 (MQ=255)
                                                                                                                                  aaCACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCgg                    <  2:416658/150‑1 (MQ=255)
                                                                                                                                     acaTCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCggaag                 >  2:1179203/1‑150 (MQ=255)
                                                                                                                                          atACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCggaagggaag            <  2:907613/150‑1 (MQ=255)
                                                                                                                                             cacacCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGaagaag         <  2:1292063/150‑1 (MQ=255)
                                                                                                                                               cacCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGAt       >  2:640527/1‑150 (MQ=255)
                                                                                                                                                    aTAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTggg  <  2:106723/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGG  >  NC_000913/2725304‑2725601

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: ‑