Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,169,657 G→A 5.8% intergenic (+46/+36) yegS → / ← gatR phosphatidylglycerol kinase, metal‑dependent/pseudogene, repressor for gat operon; interrupted by IS3; split galactitol utilization operon repressor, fragment 2; split galactitol utilization operon repressor, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,169,6570GA 5.8% 326.4 / 5.4 104intergenic (+46/+36)yegS/gatRphosphatidylglycerol kinase, metal‑dependent/pseudogene, repressor for gat operon; interrupted by IS3; split galactitol utilization operon repressor, fragment 2; split galactitol utilization operon repressor, interrupted
Reads supporting (aligned to +/- strand):  ref base G (45/52);  new base A (4/2);  total (49/55)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 108 total aligned reads displayed.

TTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATT  >  NC_000913/2169508‑2169803
                                                                                                                                                     |                                                                                                                                                  
ttAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATg                                                                                                                                                    <  2:1687511/150‑1 (MQ=255)
ttAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATg                                                                                                                                                    >  1:1577985/1‑150 (MQ=255)
ttAATCTTGATGGCGAACCGTTGAGTGGGAAAAATTTTCATATTGAAATACTTCCGTCAGCGTTTGGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATg                                                                                                                                                    <  1:1613077/150‑1 (MQ=255)
   aTCTTGATGGCGAACCGTTGAGTGGGCAAAATTATCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCgg                                                                                                                                                 <  2:422086/150‑1 (MQ=255)
     attGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGTATGCGgtg                                                                                                                                               <  1:1690162/149‑1 (MQ=255)
       tGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTa                                                                                                                                             <  2:1266877/150‑1 (MQ=255)
         aTGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTTTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTaaa                                                                                                                                           >  2:1585882/1‑150 (MQ=255)
           ggCGATCTGTTGAGGGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTAGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAaca                                                                                                                                         <  2:238435/150‑1 (MQ=255)
             cGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACAcc                                                                                                                                       <  1:1716577/150‑1 (MQ=255)
                  cGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATc                                                                                                                                  >  1:125548/1‑150 (MQ=255)
                  cGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCAGTGTAAACACCTTATc                                                                                                                                  <  2:1105772/150‑1 (MQ=255)
                       aGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCa                                                                                                                             <  2:585854/150‑1 (MQ=255)
                        gTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTACACACCTTCTCCTGCAt                                                                                                                            >  1:266367/1‑150 (MQ=255)
                        gTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGGTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATAAGGTGTTTACACCGCATCCGGCAt                                                                                                                            <  2:704355/150‑1 (MQ=255)
                         tGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGAATACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATGCTCTCATGAACCAGGGAGAATCCGTAATGCCGGATAAGGTGTTTACACCGCATCCGGCAtt                                                                                                                           <  1:1101502/150‑2 (MQ=255)
                          gggCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATAc                                                                                                                          <  2:1079408/150‑1 (MQ=255)
                           ggCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGATTCAGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACa                                                                                                                         >  2:327830/1‑150 (MQ=255)
                               aaaTTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAAc                                                                                                                     <  1:216800/150‑1 (MQ=255)
                                 aTTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACaa                                                                                                                   >  1:390626/1‑150 (MQ=255)
                                  ttttCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAAt                                                                                                                  <  1:1297793/150‑1 (MQ=255)
                                        tatTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAACGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTaca                                                                                                            >  2:585248/1‑150 (MQ=255)
                                         atTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTacac                                                                                                           >  1:1546956/1‑150 (MQ=255)
                                               aTACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaat                                                                                                     <  1:1105075/150‑1 (MQ=255)
                                                  cTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                  >  1:918776/1‑150 (MQ=255)
                                                  cTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                  >  2:1530677/1‑150 (MQ=255)
                                                    tCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTaattaatt                                                                                                <  1:472448/150‑1 (MQ=255)
                                                      cGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAAtttt                                                                                              <  2:1549800/150‑1 (MQ=255)
                                                       ggCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGTGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTa                                                                                             <  1:27331/150‑1 (MQ=255)
                                                         cAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTAcc                                                                                           >  1:1144899/1‑150 (MQ=255)
                                                          aGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTAccc                                                                                          <  2:918776/150‑1 (MQ=255)
                                                           gcgttgcgttGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCcg                                                                                         <  2:1445072/150‑1 (MQ=255)
                                                           gcgttgcgttGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCcg                                                                                         <  2:640023/150‑1 (MQ=255)
                                                           gcgttgcgttGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCcg                                                                                         <  2:1437295/150‑1 (MQ=255)
                                                           gcgttgcgttGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCcg                                                                                         >  2:611989/1‑150 (MQ=255)
                                                               tgcgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGtc                                                                                     <  1:1675274/150‑1 (MQ=255)
                                                                gcgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGtct                                                                                    <  2:1473523/150‑1 (MQ=255)
                                                                gcgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGtct                                                                                    <  1:1578817/150‑1 (MQ=255)
                                                                  gttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGtcttc                                                                                  >  1:754361/1‑150 (MQ=255)
                                                                   ttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGtcttct                                                                                 >  2:221250/1‑150 (MQ=255)
                                                                   ttgTAGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGtcttct                                                                                 <  2:1318022/150‑1 (MQ=255)
                                                                    tgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTg                                                                                <  1:795052/150‑1 (MQ=255)
                                                                     gTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGc                                                                               <  1:799193/150‑1 (MQ=255)
                                                                     gTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGc                                                                               <  2:937976/150‑1 (MQ=255)
                                                                     gTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACATTATCCGGCATACAGAACAATACTACACAGTAATTAATTTTACCCGCGTCTTCTGc                                                                               >  2:439998/1‑150 (MQ=255)
                                                                         aTTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATg                                                                           <  1:469530/150‑1 (MQ=255)
                                                                            accaccAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCa                                                                        >  1:1166857/1‑150 (MQ=255)
                                                                            accaccAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCa                                                                        <  2:420962/150‑1 (MQ=255)
                                                                             ccaccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGATCAATACTAGACCGTAATTAATTTTACCCGCGTCTTCTGCAATGACAg                                                                       >  2:1043088/1‑150 (MQ=255)
                                                                             ccaccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                       >  2:1242518/1‑150 (MQ=255)
                                                                             ccaccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                       >  2:782588/1‑150 (MQ=255)
                                                                             ccaccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                       >  2:45435/1‑150 (MQ=255)
                                                                             ccaccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                       >  2:1483605/1‑150 (MQ=255)
                                                                              caccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGt                                                                      >  2:1445199/1‑150 (MQ=255)
                                                                               accaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTTTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGtt                                                                     <  1:1585882/150‑1 (MQ=255)
                                                                               accaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGTATGCGGTGTAAACACCTTATACGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGtt                                                                     <  1:1048625/150‑1 (MQ=255)
                                                                               accaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGtt                                                                     <  1:1456789/150‑1 (MQ=255)
                                                                                ccaGATTGTCCACTGTTGCGTTATTCAAACTATTCTCACTTGAACCAGGCAGAATCCGTAATGCCGGATGCGTTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGAAAATAATTTTACCCGCGACTTCTCCAATGCCAGttt                                                                    >  1:280794/1‑150 (MQ=255)
                                                                                    aTTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATct                                                                <  1:439998/150‑2 (MQ=255)
                                                                                       gTCCACTGTTGCTTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCAcc                                                             <  2:1166130/150‑1 (MQ=255)
                                                                                            cTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                        <  1:1145886/150‑1 (MQ=255)
                                                                                            cTGTTGCGTTAATCAAACTATTATCTCTTGAAAAAGGGAGAATCAGTAATTACGGATGCGGTGTAAACACCTTATCCGGCATACAGAAAAATACTAAACCGTAATTAATTTTACCAGAGTCTTATGCAATGACAGTATATCAACCTCTGa                                                        >  1:568690/1‑150 (MQ=255)
                                                                                             tGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGTGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGac                                                       >  2:1547339/1‑150 (MQ=255)
                                                                                               ttGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCACCTCTGAcaca                                                    <  2:1247742/150‑1 (MQ=255)
                                                                                                   gTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGca                                                 <  2:1546956/150‑1 (MQ=255)
                                                                                                        tCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGTCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCatcatc                                            <  1:1118239/150‑1 (MQ=255)
                                                                                                         cAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCg                                           <  1:1029781/150‑1 (MQ=255)
                                                                                                         cAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCg                                           >  2:956599/1‑150 (MQ=255)
                                                                                                         cAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCg                                           >  1:615447/1‑150 (MQ=255)
                                                                                                          aaaCTATTCTCTCTTTAACCAGGAAGAAACCGTAATGCCGGATGAGCTGTAAACACCTTATCCGGCATACAGAACAATACTACACCCTAATTATGTTAACACGCGTCTTCTGCAATGCAAGTTTATCCACCTCTGACAAAGCATCATCgg                                          >  1:115658/1‑150 (MQ=255)
                                                                                                           aaCTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGt                                         <  2:1144899/150‑1 (MQ=255)
                                                                                                             cTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGtaa                                       <  2:1375647/150‑1 (MQ=255)
                                                                                                              tATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTGTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGtaat                                      <  2:629724/150‑1 (MQ=255)
                                                                                                              tATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGtaat                                      <  1:1530677/150‑1 (MQ=255)
                                                                                                               aTTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGtaata                                     <  2:286384/150‑1 (MQ=255)
                                                                                                               aTTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGtaata                                     <  1:1631412/150‑1 (MQ=255)
                                                                                                                    ctctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTaataataat                                >  2:1688739/1‑150 (MQ=255)
                                                                                                                    ctctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTaataataat                                >  1:764894/1‑150 (MQ=255)
                                                                                                                     tctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATg                               >  2:986203/1‑150 (MQ=255)
                                                                                                                      ctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGt                              >  1:1682766/1‑150 (MQ=255)
                                                                                                                         taaCCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCaa                           <  2:1577985/149‑1 (MQ=255)
                                                                                                                         gAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCaa                           <  1:1460363/150‑1 (MQ=255)
                                                                                                                            ccAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCAACCTCTGACACAGCATCATCGGTAATAATAATGTCAAttt                        >  2:872708/1‑150 (MQ=255)
                                                                                                                             cAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAtttt                       >  2:871712/1‑150 (MQ=255)
                                                                                                                              aGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAttttt                      >  2:572531/1‑150 (MQ=255)
                                                                                                                              aGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAttttt                      >  2:898787/1‑150 (MQ=255)
                                                                                                                              aGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAttttt                      >  1:1017556/1‑150 (MQ=255)
                                                                                                                               gggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCTGCATACAGTACAATTCTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAtttttt                     <  2:929759/150‑1 (MQ=255)
                                                                                                                                ggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTc                    <  1:1242518/150‑1 (MQ=255)
                                                                                                                                  agaATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                  >  2:1506983/1‑150 (MQ=255)
                                                                                                                                    aaTCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATAATACACCGTAATTAATTTTACCCGCGTATTCTGCAATTCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGAtt                <  1:1578076/150‑1 (MQ=255)
                                                                                                                                       ccGTAATGCCGGATGCGGTGTAATCCCCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATCCCAGTTTATCCACCTCTCACAAAGCATCATCTGTAATAATAATGTCAATTATTTCGATTGgt             >  2:727411/1‑150 (MQ=255)
                                                                                                                                        cGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACCGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGgta            >  2:465515/1‑150 (MQ=255)
                                                                                                                                        agTAATGCCTGATGCTTTTTAAACACCTTATACTTCATACAGAACAATACTACACCGTAATTAATTTTACCCTCGTCTTCTGCAATGCCATTTTATCCAACTCTGACACATCATCATCTGTAATAATAATTTCAATTTTTTCGATTGgta            <  2:568690/149‑1 (MQ=255)
                                                                                                                                         gTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGgtag           <  1:956599/150‑1 (MQ=255)
                                                                                                                                            aTGCCGGATAAGGTGTTTACACCGCATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCCCCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTAc        >  2:789330/1‑150 (MQ=255)
                                                                                                                                            aTGCCGGATAAGGTGTTTACACCGCATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTAc        >  2:1462317/1‑150 (MQ=255)
                                                                                                                                             tGCCGGATAAGGTGTTTACACCGCATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTAcc       >  1:206594/1‑150 (MQ=255)
                                                                                                                                                 ggATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGAt   >  2:1586978/1‑150 (MQ=255)
                                                                                                                                                 ctaagCGGTGTAAATAGATTATCCGGGATAAAGATCAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATTCCAGTTTATCCACCTCTGAGAATGCATCATCGGTAATAATATTGTCAATTTTTTCGATTGGTAGTACCTGAt   <  2:115658/146‑1 (MQ=255)
                                                                                                                                                  gATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGAtt  >  2:768009/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTAATCTTGATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATT  >  NC_000913/2169508‑2169803

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