Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,758,920 C→T 5.2% intergenic (+64/‑65) alpA → / → yfjI CP4‑57 prophage; DNA‑binding transcriptional activator/CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,758,9200CT 5.2% 490.5 / 9.9 154intergenic (+64/‑65)alpA/yfjICP4‑57 prophage; DNA‑binding transcriptional activator/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (68/76);  new base T (5/3);  total (74/80)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.19e-01

Only 100 of 160 total aligned reads displayed.

AGCGAATGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGCAG  >  NC_000913/2758774‑2759068
                                                                                                                                                  |                                                                                                                                                     
aGCGAATGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGATATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCAt                                                                                                                                                     <  2:141349/149‑1 (MQ=255)
   gAATGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTaa                                                                                                                                                  >  1:1420690/1‑149 (MQ=255)
   gAATGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTCACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTaa                                                                                                                                                  <  1:548988/149‑1 (MQ=255)
   gAATGCTCGGAGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTaa                                                                                                                                                  >  1:869314/1‑149 (MQ=255)
     aTGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATa                                                                                                                                                <  1:477540/149‑1 (MQ=255)
        cTCGGAGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCg                                                                                                                                             <  1:699685/149‑1 (MQ=255)
         tCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCgg                                                                                                                                            <  2:15805/149‑1 (MQ=255)
            gCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTCCGCTATTAATGGCGGaaat                                                                                                                                       <  1:990279/149‑3 (MQ=255)
              gTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAg                                                                                                                                      <  1:1017967/149‑1 (MQ=255)
              gTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAg                                                                                                                                      >  2:596064/1‑149 (MQ=255)
               tGATATATGTCGGATGGATTGAGGCCGACATTGATGAGTGGTTATCACAACGCTGTAAACTTATATGAGGGTGTTACATGAATTAATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTATTAGCGGAATAGc                                                                                                                                     >  2:752664/1‑149 (MQ=255)
               tGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGc                                                                                                                                     <  2:973012/149‑1 (MQ=255)
               tGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGc                                                                                                                                     <  2:444826/149‑1 (MQ=255)
                gAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGcc                                                                                                                                    <  1:209853/149‑1 (MQ=255)
                 aaaTCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGGaaatagca                                                                                                                                   <  2:1482577/149‑7 (MQ=255)
                 aaaTCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCa                                                                                                                                   <  2:219750/149‑1 (MQ=255)
                  aaTCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCat                                                                                                                                  <  1:1198428/149‑1 (MQ=255)
                           gatggatTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTaa                                                                                                                         >  2:572038/1‑149 (MQ=255)
                            atggatTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTaaa                                                                                                                        >  1:483159/1‑149 (MQ=255)
                               gatTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACaa                                                                                                                     >  1:1131204/1‑149 (MQ=255)
                               gatTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACaa                                                                                                                     >  1:1219739/1‑149 (MQ=255)
                                   gAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTc                                                                                                                 >  2:819868/1‑149 (MQ=255)
                                   gAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTc                                                                                                                 <  2:1364951/149‑1 (MQ=255)
                                    aGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCt                                                                                                                >  2:141553/1‑149 (MQ=255)
                                      gCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCtta                                                                                                              >  2:128970/1‑149 (MQ=255)
                                      gCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCtta                                                                                                              <  2:1406996/149‑1 (MQ=255)
                                       ccGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTCTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaa                                                                                                             <  2:340110/149‑1 (MQ=255)
                                        cGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaat                                                                                                            >  1:1326428/1‑149 (MQ=255)
                                          agaTTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaattt                                                                                                          <  2:886785/149‑1 (MQ=255)
                                            tttGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGTGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaa                                                                                                        <  2:302928/148‑1 (MQ=255)
                                             ttGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaat                                                                                                       <  2:450191/149‑1 (MQ=255)
                                                 tgaGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTaaa                                                                                                   <  1:1131629/149‑1 (MQ=255)
                                                  gaGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCTGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAAt                                                                                                  >  1:1005661/1‑149 (MQ=255)
                                                  gaGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAAt                                                                                                  >  1:640298/1‑149 (MQ=255)
                                                  gaGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAAt                                                                                                  <  1:1431302/149‑1 (MQ=255)
                                                    gTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATaa                                                                                                >  2:676215/1‑149 (MQ=255)
                                                      ggTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAaga                                                                                              <  2:1403534/149‑1 (MQ=255)
                                                       gTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAagag                                                                                             <  2:1163822/149‑1 (MQ=255)
                                                          aTCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGGTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCa                                                                                          >  1:406338/1‑149 (MQ=255)
                                                           tCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCAt                                                                                         <  2:1593971/149‑1 (MQ=255)
                                                             acaACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTAATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTg                                                                                       >  2:888592/1‑149 (MQ=255)
                                                               aaCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                                                                     <  1:1203131/149‑1 (MQ=255)
                                                                 cGCTGTAATCTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATATGCTAAACAATTTGTTAATTTAATTAATTAAGAGGCATTGCATg                                                                                   >  2:1058534/1‑149 (MQ=255)
                                                                 cGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATg                                                                                   >  2:513943/1‑149 (MQ=255)
                                                                      tAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                                                                              >  1:298696/1‑149 (MQ=255)
                                                                      tAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                                                                              >  1:1501299/1‑149 (MQ=255)
                                                                      tAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGACATTAATAGCGGAATAGACATATTTTCTAAACAATTTCTTAATTTAATTAAATAATAGTAATTTCATGTTCaa                                                                              >  1:969264/1‑149 (MQ=255)
                                                                       aaaCTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                                                                             <  1:265262/149‑1 (MQ=255)
                                                                       aaaCTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                                                                             <  2:741079/149‑1 (MQ=255)
                                                                        aaCTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATg                                                                            <  2:92764/149‑1 (MQ=255)
                                                                          cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgt                                                                          >  1:1374446/1‑149 (MQ=255)
                                                                          cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgt                                                                          >  2:614356/1‑149 (MQ=255)
                                                                          cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgt                                                                          <  2:736660/149‑1 (MQ=255)
                                                                           ttATTTGAGTGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTTTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATTTTCAATGgtc                                                                         <  1:927024/149‑1 (MQ=255)
                                                                           ttATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATACTGATATTTACGCCATTAATAGCTGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGTCATTGCATGTTCAATGgta                                                                         <  2:582562/149‑2 (MQ=255)
                                                                                 gAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCtt                                                                   >  2:233176/1‑149 (MQ=255)
                                                                                  aGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCttt                                                                  <  2:1383125/149‑1 (MQ=255)
                                                                                   gggTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                                                                 >  2:1048314/1‑149 (MQ=255)
                                                                                     gttttACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATACTGCTATTTCAGCCATTAATATCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGTTCGTCCTTTccc                                                               <  2:306333/146‑1 (MQ=255)
                                                                                     gtgtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATAttt                                                                                                                              >  1:726107/1‑86 (MQ=255)
                                                                                     gtgtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATAttt                                                                                                                              <  2:726107/86‑1 (MQ=255)
                                                                                     gtgtTACATGAATTAATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTccc                                                               <  2:301804/149‑1 (MQ=255)
                                                                                       gtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg                                                             >  1:1250333/1‑149 (MQ=255)
                                                                                        ttACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGt                                                            <  1:1372375/149‑1 (MQ=255)
                                                                                           cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGa                                                         >  2:286110/1‑149 (MQ=255)
                                                                                            aTGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGAt                                                        >  1:1355522/1‑149 (MQ=255)
                                                                                            aTGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGAt                                                        >  2:528627/1‑149 (MQ=255)
                                                                                             tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATg                                                       <  1:459739/149‑1 (MQ=255)
                                                                                              taaTTCATACTCAATTGCTGTCATCGGATTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc                                                      <  2:386453/148‑1 (MQ=255)
                                                                                                   cATACTCAATTGCTGTCATCGGCGGGTAAACATCCTGCTATTTCCGCCATTCATAGCGGAATCGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCCTTTc                                                 <  2:233819/149‑1 (MQ=255)
                                                                                                     tACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCt                                               >  2:1464944/1‑149 (MQ=255)
                                                                                                         cAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTaaaa                                           <  1:1258821/149‑1 (MQ=255)
                                                                                                          aaTTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCGTGTAGATGCAGTTCCTaaatt                                          >  2:790224/1‑147 (MQ=255)
                                                                                                              gCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATc                                      <  1:596064/149‑1 (MQ=255)
                                                                                                               cTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCa                                     <  2:1348830/149‑1 (MQ=255)
                                                                                                               cTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCa                                     <  2:1358795/149‑1 (MQ=255)
                                                                                                                   cATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATATCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGaa                                 <  2:1327909/149‑1 (MQ=255)
                                                                                                                        gAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCaa                            >  2:1209565/1‑149 (MQ=255)
                                                                                                                         aGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAAt                           >  1:823911/1‑149 (MQ=255)
                                                                                                                           tgtAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAAttt                         >  1:448038/1‑149 (MQ=255)
                                                                                                                             tAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTAt                       >  2:381031/1‑149 (MQ=255)
                                                                                                                               aaCTTCCTTCTATTTCCGCCATTAATAGCGTAATAGCCATTTTTGCTAATAAATATCTTAATTTAATAAAATTATAGGCATTGAATGTTCATTGGTCGTCCTTTCGCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGa                     <  1:752664/149‑1 (MQ=255)
                                                                                                                                aCATCCTGCTATTTCCGCCATTAATAGCGGATTAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGaa                    >  1:135729/1‑149 (MQ=255)
                                                                                                                                aCATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGTCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGaa                    <  1:184743/149‑1 (MQ=255)
                                                                                                                                aCATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGaa                    <  1:572038/149‑1 (MQ=255)
                                                                                                                                 cATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAg                   <  2:1501299/149‑1 (MQ=255)
                                                                                                                                  aTCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGt                  <  1:233176/149‑1 (MQ=255)
                                                                                                                                     atggctattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGa               >  1:144340/8‑149 (MQ=255)
                                                                                                                                      tggctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGa               >  2:391461/7‑149 (MQ=255)
                                                                                                                                       tcTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAAc             <  1:546366/148‑1 (MQ=255)
                                                                                                                                        gctattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACa            >  2:862167/5‑149 (MQ=255)
                                                                                                                                         ctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACa            >  2:154507/4‑149 (MQ=255)
                                                                                                                                          tattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGc          >  1:1494040/3‑149 (MQ=255)
                                                                                                                                                gCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGc    <  2:489098/149‑1 (MQ=255)
                                                                                                                                                gCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGAATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGc    >  2:478988/1‑149 (MQ=255)
                                                                                                                                                 ccATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGTTCGTCCTTTACCTGTAGATGCATTTCCTAAAATTATCAGGAATGAAATTTATGAAGTGCAACAGCATACGCa   >  1:814454/1‑149 (MQ=255)
                                                                                                                                                 ccATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAAtt                          <  1:959759/126‑1 (MQ=255)
                                                                                                                                                 ccATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAAtt                          >  2:959759/1‑126 (MQ=255)
                                                                                                                                                 ccATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGCa   <  2:188836/149‑1 (MQ=255)
                                                                                                                                                  cATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGCAg  <  1:528627/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCGAATGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGCAG  >  NC_000913/2758774‑2759068

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