Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,200,918 T→A 6.0% M39L (ATG→TTG)  xisE ← e14 prophage; putative excisionase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,200,9180TA 6.0% 522.8 / 14.3 162M39L (ATG→TTG) xisEe14 prophage; putative excisionase
Reads supporting (aligned to +/- strand):  ref base T (85/67);  new base A (4/6);  total (89/73)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 162 total aligned reads displayed.

CTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAG  >  NC_000913/1200773‑1201066
                                                                                                                                                 |                                                                                                                                                    
  gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCgg                                                                                                                                                >  2:1393630/1‑150 (MQ=255)
  gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCgg                                                                                                                                                >  1:1076475/1‑150 (MQ=255)
  gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGAAACAAGGCCgg                                                                                                                                                <  1:964633/150‑1 (MQ=255)
  gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGAAACAAGGCCgg                                                                                                                                                <  1:1419021/150‑1 (MQ=255)
   aaaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCggg                                                                                                                                               >  2:2352667/1‑149 (MQ=255)
   aaaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGAAACAAGGCCGGt                                                                                                                                               <  2:107721/150‑1 (MQ=255)
    aaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTg                                                                                                                                              <  1:928243/150‑1 (MQ=255)
      cATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGAAACAAGGCCGGttat                                                                                                                                            <  2:2242748/150‑4 (MQ=255)
      cATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAATAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGa                                                                                                                                            <  2:664250/150‑1 (MQ=255)
       attattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                                                                           <  2:2030640/150‑1 (MQ=255)
        ttattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACa                                                                                                                                          <  1:303223/150‑1 (MQ=255)
         tattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaa                                                                                                                                         <  1:1012358/150‑1 (MQ=255)
             cTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTTACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaata                                                                                                                                     >  1:1999304/1‑150 (MQ=255)
                aCGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaataaat                                                                                                                                  <  1:252275/150‑1 (MQ=255)
                 cgcgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGCCAAAATAAAtt                                                                                                                                 >  1:2327397/1‑150 (MQ=255)
                 cgcgCTCAGCCTGCCATAACTCAAAATCCGTTGCAAAAGAAATGACACATTTATCATGAATACAGCTGTAACCAACACCCCAACTAACTACGAGTCCTCACGCACTCGCCATAACCGGCAATCTAACATTGCAGCTGGACAAAAAAAAtt                                                                                                                                 >  2:1974344/1‑150 (MQ=255)
                       cTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCtt                                                                                                                           <  1:1283672/150‑1 (MQ=255)
                        tGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTa                                                                                                                          <  1:508160/150‑1 (MQ=255)
                         ggCTTCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACGAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAg                                                                                                                         <  1:1563931/150‑1 (MQ=255)
                         ggCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAg                                                                                                                         >  1:502427/1‑150 (MQ=255)
                         ggCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAg                                                                                                                         >  1:1728667/1‑150 (MQ=255)
                           cTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCa                                                                                                                       <  2:1055365/150‑1 (MQ=255)
                                   aCTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgacg                                                                                                               >  1:2137117/1‑150 (MQ=255)
                                    cTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                                                                              >  1:2096856/1‑150 (MQ=255)
                                     tCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCa                                                                                                             <  2:1755386/150‑1 (MQ=255)
                                      cAAAATCCGTTGCAAAAGTAATGCGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCaa                                                                                                            >  2:644779/1‑150 (MQ=255)
                                      cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCaa                                                                                                            >  1:401609/1‑150 (MQ=255)
                                      cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCaa                                                                                                            >  1:2103169/1‑150 (MQ=255)
                                      cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCaa                                                                                                            >  1:2031613/1‑150 (MQ=255)
                                        aaaTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAAtg                                                                                                          <  2:1680524/150‑1 (MQ=255)
                                         aaTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTAGCTACTCCCGAACTAACTAGGCGTCGTCACGCACTCGCCATAAGCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGAATAACGACGGAAtgt                                                                                                         <  1:1974344/150‑1 (MQ=255)
                                                tGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGacaca                                                                                                  >  1:625135/1‑150 (MQ=255)
                                                  cAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACt                                                                                                >  1:2269246/1‑150 (MQ=255)
                                                   aaaaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCACCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACtt                                                                                               >  2:1072078/1‑150 (MQ=255)
                                                     aaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgg                                                                                             <  1:521212/150‑1 (MQ=255)
                                                      aGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgga                                                                                            <  2:1581892/150‑1 (MQ=255)
                                                       gTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGAACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTggag                                                                                           <  2:1599538/150‑1 (MQ=255)
                                                        tAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTggagg                                                                                          <  2:1167598/150‑1 (MQ=255)
                                                          aTGAGTCATTTTTCTTGATTACAGGTGTTACCAACTGCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGAt                                                                                        <  1:1049894/150‑1 (MQ=255)
                                                          aTGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGAt                                                                                        <  1:1835404/150‑1 (MQ=255)
                                                             aGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTAc                                                                                     >  1:1330818/1‑150 (MQ=255)
                                                               tCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACtt                                                                                   >  1:593076/1‑150 (MQ=255)
                                                               tCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCAGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACtt                                                                                   <  1:902332/150‑1 (MQ=255)
                                                                cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTc                                                                                  >  1:1923772/1‑150 (MQ=255)
                                                                cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTc                                                                                  >  1:1032586/1‑150 (MQ=255)
                                                                cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTAGACACACTTGGAGGATTACTTc                                                                                  >  2:88621/1‑150 (MQ=255)
                                                                 aTTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAAATCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCt                                                                                 <  2:404873/150‑1 (MQ=255)
                                                                   ttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGt                                                                               <  1:325891/150‑1 (MQ=255)
                                                                   ttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGt                                                                               <  2:2124231/150‑1 (MQ=255)
                                                                       cTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTAttt                                                                           >  2:587175/1‑150 (MQ=255)
                                                                       cTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTAttt                                                                           >  1:1002253/1‑150 (MQ=255)
                                                                          gATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTc                                                                        <  1:2314269/150‑1 (MQ=255)
                                                                            ttACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcag                                                                      <  1:396201/150‑1 (MQ=255)
                                                                              aCAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcagca                                                                    <  1:2001537/150‑1 (MQ=255)
                                                                                 ggTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc                                                                 >  2:1199166/1‑150 (MQ=255)
                                                                                       aCCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttct                                                           <  1:2389778/150‑1 (MQ=255)
                                                                                       aCCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttct                                                           >  1:1159106/1‑150 (MQ=255)
                                                                                         cAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttcttc                                                         >  2:218109/1‑150 (MQ=255)
                                                                                               cccAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGt                                                   >  2:802299/1‑150 (MQ=255)
                                                                                               cccAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGt                                                   <  2:563897/150‑1 (MQ=255)
                                                                                               cccAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGt                                                   >  2:1242796/1‑150 (MQ=255)
                                                                                                ccAACTAACTCGGCGTCCTCACGCACTCGCCATAACCTGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGtt                                                  >  2:1843083/1‑150 (MQ=255)
                                                                                                 cAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTa                                                 >  2:2319194/1‑150 (MQ=255)
                                                                                                   agtaactCGGCGTCCTCTGGCACTCGCCATAACCGGCCTTGTTTGATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACAGACTTGGAGGATTACTTCTGTATTTTTCAGAAGCCGATTCTTCAAGAGTTAAc                                               <  2:1817328/148‑1 (MQ=255)
                                                                                                   actaactCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAAc                                               >  2:761158/1‑150 (MQ=255)
                                                                                                   actaactCGGCGTCCTAACTCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTTGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAAc                                               <  1:785281/150‑1 (MQ=255)
                                                                                                    ctaactCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttcttc                                                         >  1:1387163/1‑139 (MQ=255)
                                                                                                    ctaactCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACa                                              <  2:2374862/150‑1 (MQ=255)
                                                                                                       actCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACAttt                                           >  2:1504204/1‑150 (MQ=255)
                                                                                                       actCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACAttt                                           >  2:384935/1‑150 (MQ=255)
                                                                                                        ctCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTTGAGTATTACTTCTTTATTTTTCAGCATCCCATTCTTTAAGATTTAAGATTTg                                          <  2:1465148/150‑1 (MQ=255)
                                                                                                          cGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGaa                                        >  2:1348433/1‑150 (MQ=255)
                                                                                                          cGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGaa                                        <  2:133978/150‑1 (MQ=255)
                                                                                                           ggCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAg                                       <  1:1440063/150‑1 (MQ=255)
                                                                                                            gCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGaatc                                      <  2:340176/150‑3 (MQ=255)
                                                                                                               tCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCata                                   >  1:781854/1‑150 (MQ=255)
                                                                                                                ccTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCatat                                  <  2:2171394/150‑1 (MQ=255)
                                                                                                                ccTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCatat                                  >  1:2448586/1‑150 (MQ=255)
                                                                                                                  tCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCGTATGc                                >  2:859517/1‑150 (MQ=255)
                                                                                                                   cACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCg                               >  2:1083538/1‑150 (MQ=255)
                                                                                                                       cACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCa                           <  2:1360159/150‑1 (MQ=255)
                                                                                                                        aCTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTTCTTCTGTATTTTTCAGCAGCCGATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCAc                          >  2:1414114/1‑150 (MQ=255)
                                                                                                                          tCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCt                        <  2:388911/150‑1 (MQ=255)
                                                                                                                           cGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCtt                       >  1:1571529/1‑150 (MQ=255)
                                                                                                                           cGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCtt                       >  2:2155322/1‑150 (MQ=255)
                                                                                                                            gCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCtta                      <  1:218109/150‑1 (MQ=255)
                                                                                                                             ccATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTa                                                 <  2:1788205/122‑1 (MQ=255)
                                                                                                                              cATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCCCGttatt                    >  1:1618154/1‑150 (MQ=255)
                                                                                                                              cATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCttatt                    <  2:1721857/150‑1 (MQ=255)
                                                                                                                                tAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGCCGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGACCATTCTTCAAGAGATAACATTTGAAGCAAATGCGATCACCTTATtac                  >  1:1127622/1‑150 (MQ=255)
                                                                                                                                 caCCGGCCATGAAACAAGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtact                 >  1:350125/2‑150 (MQ=255)
                                                                                                                                 aaCCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtact                 >  1:2483735/1‑150 (MQ=255)
                                                                                                                                  aCCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtacta                <  1:880065/150‑1 (MQ=255)
                                                                                                                                   ccGGCCTTGTATCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtactac               <  1:1390413/150‑1 (MQ=255)
                                                                                                                                     ggCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTacac             <  1:961395/150‑1 (MQ=255)
                                                                                                                                          tGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAAGTACGCATTGTTGACACACTTGGAGGATTACTTCTGTTTTTTTCTGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACt        <  2:1012654/150‑1 (MQ=255)
                                                                                                                                          tGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACt        <  1:1729296/150‑1 (MQ=255)
                                                                                                                                              tCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCtt    >  1:1632155/1‑150 (MQ=255)
                                                                                                                                               cATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTa   >  1:1718402/1‑150 (MQ=255)
                                                                                                                                                aTGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAg  >  2:257192/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAG  >  NC_000913/1200773‑1201066

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