Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,968,379 T→A 5.6% N19I (AAT→ATT)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,968,3790TA 5.6% 361.3 / 4.5 119N19I (AAT→ATT) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base T (58/54);  new base A (5/2);  total (63/56)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

CGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTCGCG  >  NC_000913/2968230‑2968525
                                                                                                                                                     |                                                                                                                                                  
cgAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCAt                                                                                                                                                    >  2:945/1‑150 (MQ=255)
    cGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAg                                                                                                                                                >  2:39769/1‑150 (MQ=255)
      ttAAAGTGCGACCGCGGGGTGTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAAGAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATAATTAACGCCTCATTCAGGc                                                                                                                                              >  2:505108/1‑150 (MQ=255)
        aCAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCgt                                                                                                                                            <  2:720909/150‑1 (MQ=255)
          aGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCgtgg                                                                                                                                          >  1:1647436/1‑150 (MQ=255)
             gCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGAAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTGAATGAGGCgttaaat                                                                                                                                       <  2:817538/150‑6 (MQ=255)
              cGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTgcg                                                                                                                                      >  1:443019/1‑150 (MQ=255)
                    cgTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGGACAGCCCTCGGTGTCCAGGGCACGACAGATGTCGGTAACCAGCAGTTTTAACGCGTCATGCGGGCGTGGTGCGGTGGCt                                                                                                                                >  2:61701/1‑150 (MQ=255)
                        gTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCt                                                                                                                            <  2:513995/150‑1 (MQ=255)
                        gTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCt                                                                                                                            <  1:165747/150‑1 (MQ=255)
                           ttttCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCtttt                                                                                                                         <  2:124288/150‑1 (MQ=255)
                             ttCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCg                                                                                                                       >  2:167218/1‑150 (MQ=255)
                             ttCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCg                                                                                                                       <  2:1480888/150‑1 (MQ=255)
                                aGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACt                                                                                                                    <  1:908956/150‑1 (MQ=255)
                                   ccccGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTAtt                                                                                                                 >  1:1500602/1‑150 (MQ=255)
                                   ccccGGGTCGCCATCAGGTAGTAACAACGTCGATCATAATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTAtt                                                                                                                 >  1:1032077/1‑150 (MQ=255)
                                    cccGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTAttt                                                                                                                >  1:313435/1‑150 (MQ=255)
                                        ggTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTcgcg                                                                                                            <  2:13420/150‑1 (MQ=255)
                                         gTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCgcgc                                                                                                           >  1:1299921/1‑150 (MQ=255)
                                         gTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCgcgc                                                                                                           <  2:1291388/150‑1 (MQ=255)
                                                  aGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTg                                                                                                  >  1:1068336/1‑150 (MQ=255)
                                                   ggTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCTCGCAGGCGAGTGa                                                                                                 >  2:1507300/1‑150 (MQ=255)
                                                   ggTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGa                                                                                                 >  1:817714/1‑150 (MQ=255)
                                                     tagtaACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGc                                                                                               <  1:1429250/150‑1 (MQ=255)
                                                     tagtaACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGc                                                                                               >  1:1464332/1‑150 (MQ=255)
                                                     tagtaACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGc                                                                                               >  1:1304763/1‑150 (MQ=255)
                                                     tagtaACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTATTCGACTATTTCGCGCAGGCGAGTGAGc                                                                                               >  1:49423/1‑150 (MQ=255)
                                                            aacGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGTGAGTGAGCATAATTg                                                                                        >  1:1604553/1‑150 (MQ=255)
                                                            aacGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTg                                                                                        <  2:1068336/150‑1 (MQ=255)
                                                             acGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGTCGAGTGAGCATAATTgg                                                                                       <  2:1014413/150‑1 (MQ=255)
                                                             acGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTgg                                                                                       <  2:177409/150‑1 (MQ=255)
                                                                 cGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTg                                                                                   <  1:167218/150‑1 (MQ=255)
                                                                 cGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTg                                                                                   >  2:1549889/1‑150 (MQ=255)
                                                                 cGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTg                                                                                   >  2:35846/1‑150 (MQ=255)
                                                                  gATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGa                                                                                  >  2:851862/1‑150 (MQ=255)
                                                                   aTCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGAc                                                                                 >  1:722817/1‑150 (MQ=255)
                                                                       tGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTaa                                                                             >  1:1065704/1‑150 (MQ=255)
                                                                         aTCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAAcc                                                                           <  2:1032077/150‑1 (MQ=255)
                                                                         aTCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAAcc                                                                           <  1:1015861/150‑1 (MQ=255)
                                                                         aTCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAAcc                                                                           >  2:751854/1‑150 (MQ=255)
                                                                           cGGCCAGGTAGACCGAACAGCCCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATATCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACcta                                                                         >  1:1255235/1‑149 (MQ=255)
                                                                            ggCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGAGAGTGAGCATAATTGGCGTGACTTAACctct                                                                        >  2:562161/1‑150 (MQ=255)
                                                                               cAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCtttt                                                                     <  2:1500602/150‑1 (MQ=255)
                                                                                       ccGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGAt                                                             <  1:297817/150‑1 (MQ=255)
                                                                                            cAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAg                                                          >  2:1204449/1‑148 (MQ=255)
                                                                                            cAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAg                                                          >  2:1180002/1‑148 (MQ=255)
                                                                                            cAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAg                                                          <  1:1204449/148‑1 (MQ=255)
                                                                                            cAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAg                                                          <  1:1180002/148‑1 (MQ=255)
                                                                                            cAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAg                                                          <  1:1173443/148‑1 (MQ=255)
                                                                                              gACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGa                                                      >  1:734590/1‑150 (MQ=255)
                                                                                                       aTCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTtg                                             <  1:562161/150‑1 (MQ=255)
                                                                                                       aTCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTtg                                             >  1:1494621/1‑150 (MQ=255)
                                                                                                       aTCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTAGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTtg                                             >  1:1486686/1‑150 (MQ=255)
                                                                                                         ccATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTtgtg                                           <  2:956197/150‑1 (MQ=255)
                                                                                                         ccATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTtgtg                                           <  2:1157666/150‑1 (MQ=255)
                                                                                                         ccATCGCAAGACAGATGTCGGTAACCAGAATATATAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTATCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTtgtg                                           >  1:885123/1‑150 (MQ=255)
                                                                                                          cATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTgg                                          <  1:922839/150‑1 (MQ=255)
                                                                                                           aTCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGt                                         >  1:1456235/1‑150 (MQ=255)
                                                                                                             cGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGttt                                       >  1:1028362/1‑150 (MQ=255)
                                                                                                              gCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGtttt                                      <  1:1549889/150‑1 (MQ=255)
                                                                                                              gCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGtttt                                      <  1:61701/150‑1 (MQ=255)
                                                                                                              gCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGtttt                                      >  1:953568/1‑150 (MQ=255)
                                                                                                               cAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGAt                                                     >  1:990768/1‑134 (MQ=255)
                                                                                                               cAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTg                                     >  1:1208930/1‑150 (MQ=255)
                                                                                                               cAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTg                                     <  1:204031/150‑1 (MQ=255)
                                                                                                                 aGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATACGCAGATGCGTTTTGTGGTTTTGGc                                   >  2:1057446/1‑150 (MQ=255)
                                                                                                                  gACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTATCCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCg                                  >  2:725304/1‑150 (MQ=255)
                                                                                                                  gACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCg                                  <  2:49423/150‑1 (MQ=255)
                                                                                                                  gACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCAAAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCg                                  <  2:1464332/150‑1 (MQ=255)
                                                                                                                    cAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTa                                >  2:528451/1‑150 (MQ=255)
                                                                                                                    cAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTa                                <  1:851862/150‑1 (MQ=255)
                                                                                                                      gATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTAtt                              >  2:485045/1‑150 (MQ=255)
                                                                                                                          tCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTcc                          <  2:1647436/150‑1 (MQ=255)
                                                                                                                           cGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCt                         >  1:571371/1‑150 (MQ=255)
                                                                                                                           cGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCt                         <  2:952925/150‑1 (MQ=255)
                                                                                                                            ggTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTg                        >  2:21310/1‑150 (MQ=255)
                                                                                                                            ggTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTg                        >  1:753979/1‑150 (MQ=255)
                                                                                                                            ggTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTg                        <  1:974073/150‑1 (MQ=255)
                                                                                                                            ggTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTg                        <  2:1594073/150‑1 (MQ=255)
                                                                                                                               aaCCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAg                     <  2:1576343/150‑1 (MQ=255)
                                                                                                                                    gAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGGTTTGTGGTTTTGGAGTATTTTCCTGCAGTGACa                >  2:535271/1‑150 (MQ=255)
                                                                                                                                      gcgcaccaCGCCTGAATGAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATc              >  2:1714054/8‑150 (MQ=255)
                                                                                                                                       tatTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATca             <  2:551678/150‑1 (MQ=255)
                                                                                                                                        atTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATcac            <  1:1454041/150‑1 (MQ=255)
                                                                                                                                         caccaCGCCTGAATGAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTACCCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATcccc           >  1:576356/5‑150 (MQ=255)
                                                                                                                                         caccaCGCCTGAATGAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATcacc           >  1:556820/5‑150 (MQ=255)
                                                                                                                                         caccaCGCCTGAATGAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGAAGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATcacc           >  1:972916/5‑150 (MQ=255)
                                                                                                                                          ttAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATcacca          >  1:925106/1‑150 (MQ=255)
                                                                                                                                           tatcGCCTCATTCTGGCGTGGTGTGCTGGCGACCTTTTCGACTTTTTCGAGCAGGCGAGTGAGCATAATTGGCGTGACTTAACATCTTTTACGTCGATAAGCAGTTGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATcaccac         <  2:1042615/147‑1 (MQ=255)
                                                                                                                                           tAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATcaccac         <  2:722817/150‑1 (MQ=255)
                                                                                                                                           tAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATcaccac         <  2:1456235/150‑1 (MQ=255)
                                                                                                                                               gCCTTATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTc     >  2:278426/1‑150 (MQ=255)
                                                                                                                                               gCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTc     >  1:590711/1‑150 (MQ=255)
                                                                                                                                                ccTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTATCCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTcg    <  1:725304/150‑1 (MQ=255)
                                                                                                                                                ccTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTcg    <  2:1237136/150‑1 (MQ=255)
                                                                                                                                                catCATTCAGGCGGGGTGCGCTGGCTAGCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATACGCAGATGCGTTTTGAGGGTTTGGCGTATTTTCCTGCAGTGACATCACCACACTcg    <  2:313435/148‑1 (MQ=255)
                                                                                                                                                  tCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTcgcg  >  2:475622/1‑150 (MQ=255)
                                                                                                                                                  tCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTcgcg  >  2:787539/1‑150 (MQ=255)
                                                                                                                                                  tCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTcgcg  >  2:1623746/1‑150 (MQ=255)
                                                                                                                                                  tCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGAAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTcgcg  >  2:1421553/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTGCGCTGGCTACCTTTTCGACTATTTCGCGCAGGCGAGTGAGCATAATTGGCGTGACTTAACCTCTTTTACGTCGATAAGCAGATGCGTTTTGTGGTTTTGGCGTATTTTCCTGCAGTGACATCACCACACTCGCG  >  NC_000913/2968230‑2968525

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