Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,200,913 T→A 5.3% K40N (AAA→AAT) ‡ xisE ← e14 prophage; putative excisionase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,200,9130TA 5.3% 459.5 / 8.4 159K40N (AAA→AAT) ‡xisEe14 prophage; putative excisionase
Reads supporting (aligned to +/- strand):  ref base T (82/65);  new base A (2/7);  total (85/74)
Fisher's exact test for biased strand distribution p-value = 8.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 163 total aligned reads displayed.

AGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGC  >  NC_000913/1200764‑1201062
                                                                                                                                                     |                                                                                                                                                     
aGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCAATAACCGGCCTTGt                                                                                                                                                       >  1:662980/1‑150 (MQ=255)
  ttAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGttt                                                                                                                                                     <  2:1676658/150‑1 (MQ=255)
    aGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCa                                                                                                                                                   <  2:426918/150‑1 (MQ=255)
    aGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCa                                                                                                                                                   >  2:1425873/1‑150 (MQ=255)
    aGGGACTGAAACATTATTCTTACGCGGTCTGGCAGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCa                                                                                                                                                   >  2:561866/1‑150 (MQ=255)
      ggACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACAAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATg                                                                                                                                                 >  2:1503029/1‑150 (MQ=255)
       gACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATgg                                                                                                                                                <  2:1179619/150‑1 (MQ=255)
         cTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGcc                                                                                                                                              >  1:477774/1‑150 (MQ=255)
           gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCgg                                                                                                                                            >  1:366215/1‑150 (MQ=255)
            aaaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGt                                                                                                                                           >  2:635231/1‑150 (MQ=255)
            aaaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGt                                                                                                                                           >  2:1401511/1‑150 (MQ=255)
            aaaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGt                                                                                                                                           >  1:394914/1‑150 (MQ=255)
             aaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTg                                                                                                                                          <  1:1041087/150‑1 (MQ=255)
                 ttattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGAAACAAGGCCGGttatgg                                                                                                                                      <  2:6918/150‑6 (MQ=255)
                  tattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGAAACAAGGCCGGttatggc                                                                                                                                     <  2:5883/150‑7 (MQ=255)
                   attCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGAAACAAGGCCGGttatggcg                                                                                                                                    <  1:1673642/150‑8 (MQ=255)
                       ttACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaataa                                                                                                                                >  1:1637923/1‑150 (MQ=255)
                        tACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaataaa                                                                                                                               >  1:1460330/1‑150 (MQ=255)
                              gTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGc                                                                                                                         >  2:596175/1‑150 (MQ=255)
                                cTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGAACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAAAAAATAAATTCTGcct                                                                                                                       >  2:1285964/1‑148 (MQ=255)
                                      gCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATa                                                                                                                 >  1:793308/1‑150 (MQ=255)
                                         aTTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacg                                                                                                              >  1:1526692/1‑150 (MQ=255)
                                           tACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgac                                                                                                            >  2:525578/1‑150 (MQ=255)
                                            aCTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCACAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgaag                                                                                                           >  2:335412/1‑148 (MQ=255)
                                             cTCAAAATCCGTTGCAAAAGTAATGAGTCATTTATCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                                                                          <  2:864749/150‑1 (MQ=255)
                                              tCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCa                                                                                                         >  2:1676239/1‑150 (MQ=255)
                                              tCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCa                                                                                                         >  2:993674/1‑150 (MQ=255)
                                               cAAAATCCGTTGCAAAAGTAATGCGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCTCGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTATGATTAGCATAACGCCGCat                                                                                                        >  2:296417/1‑149 (MQ=255)
                                               cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCATGGTTTAATGGCCGGTGGACAAAATAAATTCTGCTAACCATAAGGACGCaa                                                                                                        >  2:482161/1‑150 (MQ=255)
                                               cAAAATCCGTTGCAAAAGTAATGAGACATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCCGCGTCCTCACGCACTCGCCATAACCGGCCATGTATCATGCCCGGTCAACAAAAAAAATTCTGCTTAGCATAACGACGCaa                                                                                                        >  2:273164/1‑150 (MQ=255)
                                                   aTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAAtgtg                                                                                                    <  2:520369/150‑1 (MQ=255)
                                                   aTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAAtgtg                                                                                                    <  1:1439770/150‑1 (MQ=255)
                                                    tCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTgg                                                                                                   <  2:1526692/150‑1 (MQ=255)
                                                      cGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGac                                                                                                 >  1:523446/1‑150 (MQ=255)
                                                      cGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGac                                                                                                 >  1:70117/1‑150 (MQ=255)
                                                         tGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGacaca                                                                                              >  1:534260/1‑150 (MQ=255)
                                                           cAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACt                                                                                            >  2:1726037/1‑150 (MQ=255)
                                                           cAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACt                                                                                            >  1:952765/1‑150 (MQ=255)
                                                           cAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACt                                                                                            >  1:1097478/1‑150 (MQ=255)
                                                            aaaaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACtt                                                                                           >  1:302921/1‑150 (MQ=255)
                                                             aaaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTg                                                                                          <  1:1174936/150‑1 (MQ=255)
                                                                gaaaTGAGTCATTTTTCTTGATTACAGTTGTTACCAACTCCCCAACTAACTCGGCGTCATCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACTCTTggag                                                                                       <  1:275287/148‑1 (MQ=255)
                                                                gTAATGAGTCATTTTTCTTGATTACAGGAGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTggag                                                                                       <  2:857790/150‑1 (MQ=255)
                                                                ctactGAGTCATTTTTCTAGATTACGGGTGTTACAAACTCACCATCTTACTCGGCGTCCTCACGCACTCGCCATAACCGGCCGTGCTTCATGGCCGGTGGACAAAATAACTTCTGCTTAGCATAACGACGCAATGTGGACACACTTggag                                                                                       <  2:1066380/146‑1 (MQ=255)
                                                                     gAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTa                                                                                  >  1:250001/1‑150 (MQ=255)
                                                                       gTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACt                                                                                >  2:828823/1‑150 (MQ=255)
                                                                         cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTc                                                                              <  1:187709/150‑1 (MQ=255)
                                                                         cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTc                                                                              <  2:1444467/150‑1 (MQ=255)
                                                                         cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTc                                                                              >  1:696455/1‑150 (MQ=255)
                                                                         cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTGTCGTGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTc                                                                              >  1:978636/1‑150 (MQ=255)
                                                                            ttttCTTGCTTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGTCGCAATGTGGACACACTTGGAGGATTACTTCTGt                                                                           <  1:130220/150‑1 (MQ=255)
                                                                            ttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGt                                                                           <  2:952765/150‑1 (MQ=255)
                                                                             tttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTa                                                                          >  2:1162916/1‑150 (MQ=255)
                                                                               tCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTAtt                                                                        >  2:913241/1‑150 (MQ=255)
                                                                                cTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTAttt                                                                       >  1:1444435/1‑150 (MQ=255)
                                                                                 ttGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTAtttt                                                                      >  1:781400/1‑150 (MQ=255)
                                                                                  tGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTAttttt                                                                     >  2:1305377/1‑150 (MQ=255)
                                                                                     ttACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcag                                                                  <  1:1203695/150‑1 (MQ=255)
                                                                                       aCAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcagca                                                                <  2:1340086/150‑1 (MQ=255)
                                                                                          ggTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAAGCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc                                                             <  2:1637923/150‑1 (MQ=255)
                                                                                          ggTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc                                                             <  2:366215/150‑1 (MQ=255)
                                                                                          ggTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc                                                             >  1:1715796/1‑150 (MQ=255)
                                                                                          ggTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc                                                             >  2:707757/1‑150 (MQ=255)
                                                                                          ggTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc                                                             >  2:1020040/1‑150 (MQ=255)
                                                                                          ggTGTTACCAACTCCCCAACTAACTCGGCGTACTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGcc                                                             >  2:387760/1‑150 (MQ=255)
                                                                                           gtgtTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGccc                                                            <  1:619164/150‑1 (MQ=255)
                                                                                            tgtTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCa                                                           <  2:1598803/150‑1 (MQ=255)
                                                                                               tACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttc                                                        <  1:1380520/150‑1 (MQ=255)
                                                                                               tACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAtta                                                        >  2:815493/1‑149 (MQ=255)
                                                                                               tACCAACTCCCCAACTAACTCGGCGTCCACACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACCCTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttc                                                        >  2:1053135/1‑150 (MQ=255)
                                                                                                  cAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttcttc                                                     >  2:658009/1‑150 (MQ=255)
                                                                                                   aaCTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCa                                                    <  1:1503029/150‑1 (MQ=255)
                                                                                                    aCTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCGGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAGTCTTCaa                                                   >  2:697057/1‑150 (MQ=255)
                                                                                                     cTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAag                                                  <  2:1321887/150‑1 (MQ=255)
                                                                                                     cTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAag                                                  >  2:1337322/1‑150 (MQ=255)
                                                                                                        cccAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGt                                               >  2:1381473/1‑150 (MQ=255)
                                                                                                          cAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTa                                             >  1:998699/1‑150 (MQ=255)
                                                                                                               aactCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACAtt                                        <  1:1544916/150‑1 (MQ=255)
                                                                                                                actCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACAttt                                       <  2:214821/150‑1 (MQ=255)
                                                                                                                 ctCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTg                                      <  2:1444435/150‑1 (MQ=255)
                                                                                                                 ctCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTg                                      >  2:237408/1‑150 (MQ=255)
                                                                                                                    ggCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAg                                   <  2:596910/150‑1 (MQ=255)
                                                                                                                       gTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGCCGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCat                                >  1:186230/1‑150 (MQ=255)
                                                                                                                       gTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCat                                <  1:1250990/150‑1 (MQ=255)
                                                                                                                          cTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATg                             >  2:1627666/1‑150 (MQ=255)
                                                                                                                            cACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCg                           >  1:529055/1‑150 (MQ=255)
                                                                                                                               gCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATc                        <  2:793308/150‑1 (MQ=255)
                                                                                                                                 aCTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCAc                      >  2:947693/1‑150 (MQ=255)
                                                                                                                                  cTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCAcc                     <  1:635231/150‑1 (MQ=255)
                                                                                                                                     gCCATAACCGGCCTTGTTTCATGGACGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAATCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCtta                  >  2:1538738/1‑150 (MQ=255)
                                                                                                                                       cATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCAGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGATTTAACATTTGAAGCATATGCGATCACCttatt                <  2:1260722/150‑1 (MQ=255)
                                                                                                                                         tAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtac              <  2:1097478/150‑1 (MQ=255)
                                                                                                                                          aaCCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtact             <  1:1726037/150‑1 (MQ=255)
                                                                                                                                               gCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACt        <  1:913241/150‑1 (MQ=255)
                                                                                                                                                ccTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCTTATCGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTa       <  1:596175/150‑1 (MQ=255)
                                                                                                                                                ccTTGTTTAATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGGATATGCGATCACCTTATtactaaaata       >  1:457142/1‑145 (MQ=255)
                                                                                                                                                 cTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTaa      >  2:747197/1‑150 (MQ=255)
                                                                                                                                                 cTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTaa      <  1:1667983/150‑1 (MQ=255)
                                                                                                                                                   tGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACt    <  2:1024958/150‑1 (MQ=255)
                                                                                                                                                     tttCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGc  >  1:283854/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGC  >  NC_000913/1200764‑1201062

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