Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,200,918 T→A 5.8% 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 5.8% 390.8 / 9.9 121M39L (ATG→TTG) xisEe14 prophage; putative excisionase
Reads supporting (aligned to +/- strand):  ref base T (62/52);  new base A (3/4);  total (65/56)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 125 total aligned reads displayed.

GGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGT  >  NC_000913/1200769‑1201067
                                                                                                                                                     |                                                                                                                                                     
gggACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCAt                                                                                                                                                       >  1:606249/1‑150 (MQ=255)
      gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCgg                                                                                                                                                 >  2:350185/1‑150 (MQ=255)
      gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCgg                                                                                                                                                 >  1:1488288/1‑150 (MQ=255)
      gAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCACCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCATTGTTTCATGGCCgg                                                                                                                                                 >  2:1018205/1‑150 (MQ=255)
       aaaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGAAACAAGGCCGGt                                                                                                                                                <  2:105255/150‑1 (MQ=255)
        aaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTAGTTTCATGGCCGGTg                                                                                                                                               >  2:647011/1‑150 (MQ=255)
         aCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGGAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTgg                                                                                                                                              <  2:541250/150‑1 (MQ=255)
         aCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTgg                                                                                                                                              <  1:526399/150‑1 (MQ=255)
          cATTATTCTTACGCGGTCTGGCTGCCATTACTCAAACTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGGCGGTGGa                                                                                                                                             >  2:593704/1‑150 (MQ=255)
           attattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                                                                            <  2:482051/150‑1 (MQ=255)
           attattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                                                                            <  1:176748/150‑1 (MQ=255)
           attattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                                                                            <  2:1488288/150‑1 (MQ=255)
           attattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGAAACAAGGCCGGttatg                                                                                                                                            <  2:885814/150‑5 (MQ=255)
           attattCTTACGCGGGCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                                                                            <  1:673874/150‑1 (MQ=255)
              attCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                                                                                         <  2:1038552/150‑1 (MQ=255)
               ttCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaaa                                                                                                                                        >  2:307521/1‑150 (MQ=255)
                tCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaat                                                                                                                                       <  2:918933/150‑1 (MQ=255)
                  ttACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTAGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaataa                                                                                                                                     >  1:143895/1‑150 (MQ=255)
                     cgcgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                                                                                  >  2:550627/1‑150 (MQ=255)
                         gTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATCAGTCATTTTTCTTGATTACATGTGTTACCAACTACCCAACTAACTCTGCCTACTCACGCACTAGCCATAACCGGCCTTGTTTCATGCCCGGTCGACAAAATAAATTATGc                                                                                                                              >  2:1409189/1‑150 (MQ=255)
                              gCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGc                                                                                                                         <  1:1099521/150‑1 (MQ=255)
                              gCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGc                                                                                                                         <  1:944835/150‑1 (MQ=255)
                                 gCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATa                                                                                                                      >  1:64585/1‑150 (MQ=255)
                                  ccATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATaa                                                                                                                     >  2:1197774/1‑150 (MQ=255)
                                   cATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAac                                                                                                                    >  2:788642/1‑150 (MQ=255)
                                        cTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                                                                               >  1:1521598/1‑150 (MQ=255)
                                          cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                                                                               <  2:582332/148‑1 (MQ=255)
                                          cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                                                                               >  1:582332/1‑148 (MQ=255)
                                          cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCaa                                                                                                             >  2:795132/1‑150 (MQ=255)
                                          cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCaa                                                                                                             >  1:505252/1‑150 (MQ=255)
                                               tCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTgg                                                                                                        <  1:187924/150‑1 (MQ=255)
                                               tCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTgg                                                                                                        <  1:1190289/150‑1 (MQ=255)
                                               tCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTgg                                                                                                        >  1:443255/1‑150 (MQ=255)
                                                 cGTTGCGAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGac                                                                                                      >  2:112394/1‑150 (MQ=255)
                                                   ttGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGacac                                                                                                    <  2:709149/150‑1 (MQ=255)
                                                    tGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGTACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGacaca                                                                                                   >  1:434812/1‑150 (MQ=255)
                                                     gCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGacacac                                                                                                  >  1:1417323/1‑150 (MQ=255)
                                                     gCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGacacac                                                                                                  >  1:1438217/1‑150 (MQ=255)
                                                      cAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACt                                                                                                 >  1:677562/1‑150 (MQ=255)
                                                        aaaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTg                                                                                               >  1:456824/1‑150 (MQ=255)
                                                        aaaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACAAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACAGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTATGCTTAGCATAACGACGCAATGTGGACACACTTg                                                                                               >  1:570369/1‑150 (MQ=255)
                                                         aaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgg                                                                                              >  1:138457/1‑150 (MQ=255)
                                                             aaTGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTggagga                                                                                          >  2:1505437/1‑150 (MQ=255)
                                                               tGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGAtt                                                                                        >  1:1490458/1‑150 (MQ=255)
                                                                 aGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTAc                                                                                      >  1:345000/1‑150 (MQ=255)
                                                                  gTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACt                                                                                     >  2:721801/1‑150 (MQ=255)
                                                                    cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTc                                                                                   <  1:1354876/150‑1 (MQ=255)
                                                                    cATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTc                                                                                   >  2:347254/1‑150 (MQ=255)
                                                                      tttttCTTGATTACTGGTGTTACGAACTCCCCAACTAACTCGTCGTCCTCACGGACTCGCCATTACCGGCCTTGTTTTATGGAGGGTGGACAAAATAAATTCTGCTTAGCATATGGATGCTATGTGGTCACAGTTGGATGATTAATTCTg                                                                                 <  2:324005/150‑1 (MQ=255)
                                                                             tGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTAttttt                                                                          >  2:752973/1‑150 (MQ=255)
                                                                              taTTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCTCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGTACACACTTGGAGGATTACTTCTGTATTTTTc                                                                         <  1:189125/149‑1 (MQ=255)
                                                                               aTTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTca                                                                        <  1:1158906/150‑1 (MQ=255)
                                                                                ttACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcag                                                                       >  1:578686/1‑150 (MQ=255)
                                                                                 tACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcagc                                                                      <  1:1468198/150‑1 (MQ=255)
                                                                                 tACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcagc                                                                      <  1:112394/150‑1 (MQ=255)
                                                                                          tACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttc                                                             >  1:1038913/1‑150 (MQ=255)
                                                                                          tACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGCACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttc                                                             >  1:973560/1‑150 (MQ=255)
                                                                                          tACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAAAGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttc                                                             >  1:1212569/1‑150 (MQ=255)
                                                                                            ccAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttctt                                                           >  2:250214/1‑150 (MQ=255)
                                                                                              aaCTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCa                                                         <  1:75503/150‑1 (MQ=255)
                                                                                              aaCTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCa                                                         <  1:190548/150‑1 (MQ=255)
                                                                                                cTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAag                                                       >  1:797922/1‑150 (MQ=255)
                                                                                                   cccAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGGCGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGt                                                    <  1:1407503/150‑1 (MQ=255)
                                                                                                   cccAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGAAATGTCGACACACTTTGAGGATTACTTCTGTATTTTTCAGCAGCACATTCTTCAAGAGt                                                    <  2:213043/150‑1 (MQ=255)
                                                                                                    ccAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGtt                                                   >  2:1339498/1‑150 (MQ=255)
                                                                                                     cAACTAACTTGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTa                                                  >  2:865460/1‑150 (MQ=255)
                                                                                                     cAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCa                                                         >  2:1209413/1‑143 (MQ=255)
                                                                                                     cAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTa                                                  <  2:138457/150‑1 (MQ=255)
                                                                                                       actaactCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAagag                                                     <  1:780438/145‑1 (MQ=255)
                                                                                                       actaactCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAagag                                                     >  2:780438/1‑145 (MQ=255)
                                                                                                       actaactCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAAc                                                <  2:1034718/150‑1 (MQ=255)
                                                                                                            ctCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACGCACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTg                                           >  1:698820/1‑150 (MQ=255)
                                                                                                             taggCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACTCAATGTGGACACAGTTGGATTATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGa                                          <  2:143895/148‑1 (MQ=255)
                                                                                                              cGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGaa                                         <  2:64585/150‑1 (MQ=255)
                                                                                                               ggCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAg                                        >  1:1041583/1‑150 (MQ=255)
                                                                                                                   tCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCata                                    <  1:488462/150‑1 (MQ=255)
                                                                                                                     cTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATg                                  >  2:48635/1‑150 (MQ=255)
                                                                                                                      tCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGc                                 >  2:778134/1‑150 (MQ=255)
                                                                                                                       cACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCg                                >  2:1398421/1‑150 (MQ=255)
                                                                                                                       cACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCg                                <  1:865460/150‑1 (MQ=255)
                                                                                                                          gCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATc                             <  2:345000/150‑1 (MQ=255)
                                                                                                                           cACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCa                            >  1:1030400/1‑150 (MQ=255)
                                                                                                                           cACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCa                            >  1:1429223/1‑150 (MQ=255)
                                                                                                                              tCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCt                         <  1:688645/150‑1 (MQ=255)
                                                                                                                               cGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCtt                        >  1:931410/1‑150 (MQ=255)
                                                                                                                                 ccATAAACGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATTCGATCACCttat                      >  2:758405/1‑150 (MQ=255)
                                                                                                                                  gtccaCCGGCCATGAAACAAGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCttatt                     >  1:487929/5‑150 (MQ=255)
                                                                                                                                   aTAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCttatta                    <  1:1018205/150‑1 (MQ=255)
                                                                                                                                    tAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtac                   <  1:810179/150‑1 (MQ=255)
                                                                                                                                    tAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtac                   <  1:307521/150‑1 (MQ=255)
                                                                                                                                    tAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtac                   <  1:250214/150‑1 (MQ=255)
                                                                                                                                    ccaCCGGCCATGAAACAAGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtac                   >  2:956068/3‑150 (MQ=255)
                                                                                                                                         ggCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTacac              <  1:350185/150‑1 (MQ=255)
                                                                                                                                           ccTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTa            >  2:158593/1‑150 (MQ=255)
                                                                                                                                            cTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTaa           >  2:1478515/1‑150 (MQ=255)
                                                                                                                                              tGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCATCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACt         <  2:114252/150‑1 (MQ=255)
                                                                                                                                                tttCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGc       >  1:733120/1‑150 (MQ=255)
                                                                                                                                                    aTGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAg   <  1:66685/150‑1 (MQ=255)
                                                                                                                                                     tGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGt  >  1:823766/1‑150 (MQ=255)
                                                                                                                                                     tGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGt  >  1:443133/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGT  >  NC_000913/1200769‑1201067

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: