Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,169,665 A→T 5.9% intergenic (+54/+28) 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,6650AT 5.9% 299.0 / 5.4 102intergenic (+54/+28)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 A (44/51);  new base T (4/2);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 103 total aligned reads displayed.

ATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCG  >  NC_000913/2169517‑2169812
                                                                                                                                                    |                                                                                                                                                   
  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)
             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)
                                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)
                                             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:611989/1‑150 (MQ=255)
                                                  gcgttgcgttGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCcg                                                                                                  <  2:640023/150‑1 (MQ=255)
                                                  gcgttgcgttGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCcg                                                                                                  <  2:1437295/150‑1 (MQ=255)
                                                  gcgttgcgttGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCcg                                                                                                  <  2:1445072/150‑1 (MQ=255)
                                                      tgcgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGtc                                                                                              <  1:1675274/150‑1 (MQ=255)
                                                       gcgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGtct                                                                                             <  1:1578817/150‑1 (MQ=255)
                                                       gcgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGtct                                                                                             <  2:1473523/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)
                                                                              gTCCACTGTTGCTTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCAcc                                                                      <  2:1166130/150‑1 (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                                                    <  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                                              <  1:1631412/150‑1 (MQ=255)
                                                                                                      aTTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGtaata                                              <  2:286384/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)
                                                                                                                   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)
                                                                                                                                                    |                                                                                                                                                   
ATGGCGAACCGTTGAGTGGGCAAAATTTTCATATTGAAATACTTCCGGCAGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCG  >  NC_000913/2169517‑2169812

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: