Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,169,672 T→C 5.8% intergenic (+61/+21) 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,6720TC 5.8% 304.3 / 7.3 104intergenic (+61/+21)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 T (48/50);  new base C (4/2);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

Only 100 of 106 total aligned reads displayed.

CGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAA  >  NC_000913/2169526‑2169821
                                                                                                                                                  |                                                                                                                                                     
cGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATc                                                                                                                                                    >  1:125548/1‑150 (MQ=255)
cGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCAGTGTAAACACCTTATc                                                                                                                                                    <  2:1105772/150‑1 (MQ=255)
    tattGGTCAAAATTTTCATATTGAAATACTTCCTGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATTATTTGTAAACACCTTATCCGGc                                                                                                                                                <  1:1444662/147‑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)
                             aTACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaat                                                                                                                       <  1:1105075/150‑1 (MQ=255)
                                cTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                                    >  1:918776/1‑150 (MQ=255)
                                cTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGtaattaa                                                                                                                    >  2:1530677/1‑150 (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:1437295/150‑1 (MQ=255)
                                         gcgttgcgttGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCcg                                                                                                           <  2:1445072/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                                                                                                 <  2:937976/150‑1 (MQ=255)
                                                   gTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGc                                                                                                 <  1:799193/150‑1 (MQ=255)
                                                   gTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACATTATCCGGCATACAGAACAATACTACACAGTAATTAATTTTACCCGCGTCTTCTGc                                                                                                 >  2:439998/1‑150 (MQ=255)
                                                       aTTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATg                                                                                             <  1:469530/150‑1 (MQ=255)
                                                          accaccAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCa                                                                                          <  2:420962/150‑1 (MQ=255)
                                                          accaccAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCa                                                                                          >  1:1166857/1‑150 (MQ=255)
                                                           ccaccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGATCAATACTAGACCGTAATTAATTTTACCCGCGTCTTCTGCAATGACAg                                                                                         >  2:1043088/1‑150 (MQ=255)
                                                           ccaccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                                         >  2:1242518/1‑150 (MQ=255)
                                                           ccaccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                                         >  2:1483605/1‑150 (MQ=255)
                                                           ccaccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                                         >  2:45435/1‑150 (MQ=255)
                                                           ccaccaGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                                         >  2:782588/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)
                                                                       ccACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACctc                                                                             <  2:154091/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                                                             >  2:956599/1‑150 (MQ=255)
                                                                                       cAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCg                                                             <  1:1029781/150‑1 (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                                                       <  1:1631412/150‑1 (MQ=255)
                                                                                             aTTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGtaata                                                       <  2:286384/150‑1 (MQ=255)
                                                                                                  ctctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTaataataat                                                  >  1:764894/1‑150 (MQ=255)
                                                                                                  ctctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTaataataat                                                  >  2:1688739/1‑150 (MQ=255)
                                                                                                   tctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATg                                                 >  2:986203/1‑150 (MQ=255)
                                                                                                    ctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGt                                                >  1:1682766/1‑150 (MQ=255)
                                                                                                       gAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCaa                                             <  1:1460363/150‑1 (MQ=255)
                                                                                                          ccAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAttt                                          >  1:1688247/1‑150 (MQ=255)
                                                                                                          ccAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCAACCTCTGACACAGCATCATCGGTAATAATAATGTCAAttt                                          >  2:872708/1‑150 (MQ=255)
                                                                                                           cAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAtttt                                         >  2:871712/1‑150 (MQ=255)
                                                                                                            aGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAttttt                                        >  2:898787/1‑150 (MQ=255)
                                                                                                            aGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAttttt                                        >  1:1017556/1‑150 (MQ=255)
                                                                                                            aGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAttttt                                        >  2:572531/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)
                                                                                                                         aaTGCCGGATGTGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGgtagta                           >  2:705677/1‑150 (MQ=255)
                                                                                                                          aTGCCGGATAAGGTGTTTACACCGCATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCCCCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTAc                          >  2:789330/1‑150 (MQ=255)
                                                                                                                          aTGCCGGATAAGGTGTTTACACCGCATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTAc                          >  2:1462317/1‑150 (MQ=255)
                                                                                                                          aTGCCGGATAAGGTGTTTACACCGCATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTAc                          >  2:331456/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)
                                                                                                                                    cGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTaaaa                >  1:566211/1‑150 (MQ=255)
                                                                                                                                      gtgtAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAAcc              >  1:506184/1‑150 (MQ=255)
                                                                                                                                         tAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACcgcg           >  2:1712998/1‑150 (MQ=255)
                                                                                                                                             cacCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACgg       >  1:1706704/1‑150 (MQ=255)
                                                                                                                                               ccTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCCCCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGtt     >  2:153881/1‑150 (MQ=255)
                                                                                                                                               ccTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAAc               >  1:710732/1‑140 (MQ=255)
                                                                                                                                               ccTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAAc               <  2:710732/140‑1 (MQ=255)
                                                                                                                                                  tATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTaaa  >  2:901844/1‑150 (MQ=255)
                                                                                                                                                  tATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCAGCGTCTTCTGAAATGACAGTTTATACAACTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTATATTGGTAGTACCTGAATAAAACAGCGACGGTTaaa  >  1:1716494/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAA  >  NC_000913/2169526‑2169821

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: