Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,200,910 T→A 6.4% Q41H (CAA→CAT xisE ← e14 prophage; putative excisionase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,200,9100TA 6.4% 416.0 / 13.7 151Q41H (CAA→CATxisEe14 prophage; putative excisionase
Reads supporting (aligned to +/- strand):  ref base T (71/69);  new base A (6/4);  total (78/73)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

Only 100 of 151 total aligned reads displayed.

AAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAA  >  NC_000913/1200763‑1201058
                                                                                                                                                   |                                                                                                                                                    
aaGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTg                                                                                                                                                    >  2:1780511/1‑150 (MQ=255)
   ttAGTGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGgtt                                                                                                                                                 <  1:90739/150‑1 (MQ=255)
    tAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTATTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGttta                                                                                                                                                >  1:148148/1‑149 (MQ=255)
       ggACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATg                                                                                                                                             <  1:1044597/150‑1 (MQ=255)
        gACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATgg                                                                                                                                            <  1:627473/150‑1 (MQ=255)
          cTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCCCGCACTCGCCATAACCGGCCTTGTTTCATGGcc                                                                                                                                          >  2:460393/1‑150 (MQ=255)
             aaaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGt                                                                                                                                       <  2:323632/150‑1 (MQ=255)
             aaaCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGAAACAAGGCCGGt                                                                                                                                       <  2:1248138/150‑1 (MQ=255)
               aCATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTgg                                                                                                                                     <  2:97958/150‑1 (MQ=255)
                cATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATCGCCGGTggc                                                                                                                                    >  1:296975/1‑149 (MQ=255)
                cATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTAGCAAAAGAAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGa                                                                                                                                    <  1:1610162/150‑1 (MQ=255)
                 attattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGAc                                                                                                                                   <  1:1685413/150‑1 (MQ=255)
                  ttattCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACa                                                                                                                                  >  2:230507/1‑150 (MQ=255)
                    attCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACaaa                                                                                                                                >  1:57567/1‑150 (MQ=255)
                      tCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTATCTTGATTACAGGTGTGACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaat                                                                                                                              <  1:1341723/150‑1 (MQ=255)
                      tCTTACGCGGTCTGGCGGCCATTACTCAACATCCGGTGCAAAAGTAATGAGTCAGTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCTCGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaat                                                                                                                              <  1:1349230/150‑1 (MQ=255)
                       cTTACTCGTTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaata                                                                                                                             <  2:238718/150‑1 (MQ=255)
                        ttACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAaaataa                                                                                                                            >  2:1022022/1‑150 (MQ=255)
                           cgcgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAAtt                                                                                                                         <  1:141546/150‑1 (MQ=255)
                            gcgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTc                                                                                                                        <  1:543351/150‑1 (MQ=255)
                            gcgGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTc                                                                                                                        >  2:203244/1‑150 (MQ=255)
                                tCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCCCGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCt                                                                                                                    >  2:945801/1‑150 (MQ=255)
                                   ggCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAg                                                                                                                 <  2:1369811/150‑1 (MQ=255)
                                    gCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGc                                                                                                                <  2:608542/150‑1 (MQ=255)
                                     cTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCa                                                                                                               <  1:1040431/150‑1 (MQ=255)
                                        ccATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATaa                                                                                                            >  2:249019/1‑150 (MQ=255)
                                         cATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAac                                                                                                           >  2:1078645/1‑150 (MQ=255)
                                            tACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgac                                                                                                        <  1:1177660/150‑1 (MQ=255)
                                            tACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAacgac                                                                                                        >  1:519567/1‑150 (MQ=255)
                                             aCTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGCACAAAATAAATTCTGCTTATCATAacgacg                                                                                                       >  2:311952/1‑150 (MQ=255)
                                              cTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGc                                                                                                      >  1:856637/1‑150 (MQ=255)
                                                cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCTTAACGACGCaa                                                                                                    >  2:1306572/1‑150 (MQ=255)
                                                cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACTACGCaa                                                                                                    >  1:1097624/1‑150 (MQ=255)
                                                cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCaa                                                                                                    >  2:1152140/1‑150 (MQ=255)
                                                cAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCaa                                                                                                    >  2:718750/1‑150 (MQ=255)
                                                 aaaaTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAAt                                                                                                   >  1:1125734/1‑150 (MQ=255)
                                                  aaaTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAAtg                                                                                                  >  2:949700/1‑150 (MQ=255)
                                                    aTCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCCGCCTTGTTTCATGGCCGGAGGACAAAATAAAATCTGCTTAGCATAACGACGCAAtgtg                                                                                                >  2:1344928/1‑150 (MQ=255)
                                                     tCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTgg                                                                                               <  2:675350/150‑1 (MQ=255)
                                                     tCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTgg                                                                                               <  2:184538/150‑1 (MQ=255)
                                                     tCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTgg                                                                                               <  1:1390927/150‑1 (MQ=255)
                                                          tGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGacaca                                                                                          <  2:1780679/150‑1 (MQ=255)
                                                           gCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCAGAACGCCGCAATGTGGacaccc                                                                                         >  1:931883/1‑148 (MQ=255)
                                                            cAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACt                                                                                        >  1:1692320/1‑150 (MQ=255)
                                                            cAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGAAGCAATGTGGACACACt                                                                                        >  1:1731229/1‑150 (MQ=255)
                                                             aaaaGTAATGAGTCATTTTTCTTGATTATAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGGCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGCCGCAATGTGGACACACtt                                                                                       >  1:64631/1‑150 (MQ=255)
                                                             aaaaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACtt                                                                                       <  1:810300/150‑1 (MQ=255)
                                                             aaaaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACtt                                                                                       >  1:1467457/1‑150 (MQ=255)
                                                               aaGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTgg                                                                                     <  1:392746/150‑1 (MQ=255)
                                                                     tGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGAtt                                                                               <  2:662975/150‑1 (MQ=255)
                                                                      gAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACAAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTa                                                                              <  2:1737140/150‑1 (MQ=255)
                                                                            tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTg                                                                        <  2:383606/150‑1 (MQ=255)
                                                                            tttttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTAGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTg                                                                        <  1:1331881/150‑1 (MQ=255)
                                                                              tttCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTa                                                                      <  1:1745799/150‑1 (MQ=255)
                                                                                     aTTACAGGTGTTACCAACTCCCCAACTACCTCGGCGTCCTCACGCACTCGCCATAACCGGCATTGGTTCATGGCAGTTGGCGAAAATAAATTCTGCTTAGAAGACCTCCGCATTGTGGCCCCAGGTGGAGGATGACTTCTGTATTTTTcc                                                               >  1:299480/1‑149 (MQ=255)
                                                                                      ttACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcag                                                              <  1:412768/150‑1 (MQ=255)
                                                                                       tACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcagc                                                             <  2:1322671/150‑1 (MQ=255)
                                                                                       tACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcagc                                                             <  2:1392131/150‑1 (MQ=255)
                                                                                       tACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcagc                                                             >  1:447690/1‑150 (MQ=255)
                                                                                       tACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTGTCATGGCCGGTGGACATAATAAATGCTGCTTAGCATAACGCCGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTcagc                                                             >  1:328378/1‑150 (MQ=255)
                                                                                         cAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTATCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGAAGCACTGGGAGGATTACTTCTCTATTTTTcagcat                                                           >  1:316600/1‑149 (MQ=255)
                                                                                                 aCCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttct                                                   >  2:554928/1‑150 (MQ=255)
                                                                                                  ccAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttctt                                                  >  1:1259933/1‑150 (MQ=255)
                                                                                                   cAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCAttcttc                                                 >  2:1527675/1‑150 (MQ=255)
                                                                                                     aCTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCaa                                               <  1:982894/150‑1 (MQ=255)
                                                                                                      cTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAag                                              <  2:803109/150‑1 (MQ=255)
                                                                                                        ccccAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAagag                                            >  1:1530213/1‑150 (MQ=255)
                                                                                                        ccccAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCATGTGTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTTTTCAagag                                            >  1:1277766/1‑150 (MQ=255)
                                                                                                           cAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTa                                         >  2:751138/1‑150 (MQ=255)
                                                                                                           cAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTa                                         >  1:979175/1‑150 (MQ=255)
                                                                                                            aactaactCGGCGTCCTCACGCTCTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTc                                                                <  2:766673/126‑1 (MQ=255)
                                                                                                            aactaactCGGCGTCCTCACGCGCTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTc                                                                >  1:766673/1‑126 (MQ=255)
                                                                                                             actaactCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAAc                                       <  1:949700/150‑1 (MQ=255)
                                                                                                                aactCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATATCGACGCAATGTGGCGACACTTGGAGGAGTACTTCTGTATTTTTCAGCAGCCCGTTCTTCAAGAGTTAACAtt                                    >  2:901004/1‑150 (MQ=255)
                                                                                                                aactCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACAtt                                    >  2:197092/1‑150 (MQ=255)
                                                                                                                  ctCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTg                                  <  1:1176142/150‑1 (MQ=255)
                                                                                                                  ctCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTg                                  <  1:1782932/150‑1 (MQ=255)
                                                                                                                   tCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGa                                 >  2:332963/1‑150 (MQ=255)
                                                                                                                     ggCGTCCTCACGCACTCGACATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAg                               <  2:1185337/150‑1 (MQ=255)
                                                                                                                      gCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACCCTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGc                              >  2:1445725/1‑150 (MQ=255)
                                                                                                                      gCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGc                              >  2:802976/1‑150 (MQ=255)
                                                                                                                           cTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATg                         >  2:1069453/1‑150 (MQ=255)
                                                                                                                                ggtcTCGCCAAAACCGGCCTTGTTTTATGGCCGGTGTACAAAATAAATTCTGTTTAGCATAACGACGCAATGTGGACACACTTGGTTGAATACTTCTTTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATc                    <  2:182770/147‑1 (MQ=255)
                                                                                                                                   cTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGCCGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCCGCCCATTCTTCAAGAGTGAACATTTGAAGCAGATGCGAGCAcc                 >  2:304108/1‑150 (MQ=255)
                                                                                                                                   cTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCAcc                 >  1:1549898/1‑150 (MQ=255)
                                                                                                                                   cTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCAcc                 <  2:1467457/150‑1 (MQ=255)
                                                                                                                                   cTCGCCATAACATGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCAcc                 <  1:1101943/150‑1 (MQ=255)
                                                                                                                                    tCGCCATAACCGGACTTGTTTCATGGCCGGTGGACAAAAAAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCt                <  1:311952/150‑1 (MQ=255)
                                                                                                                                       ccATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGATTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCttat             <  2:708981/150‑1 (MQ=255)
                                                                                                                                        gtccaCCGGCCATGAAACAAGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCttatt            >  2:1680933/5‑150 (MQ=255)
                                                                                                                                        cATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCttatt            <  1:736335/150‑1 (MQ=255)
                                                                                                                                         tccaCCGGCCATGAAACAAGGCCGGTGGCCAAAATAAATTCTGCTTAGCATAACGACGCAATGTGAACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCCTATGCGCTCACCttatta           >  1:164673/4‑150 (MQ=255)
                                                                                                                                           caCCGGCCATGAAACAAGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtact         >  2:495888/2‑150 (MQ=255)
                                                                                                                                           caCCGGCCATGAAACAAGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGAGTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAAGATTTGAAGCATATGCGATCACCTTATtact         >  2:302141/2‑150 (MQ=255)
                                                                                                                                            aCCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTATCATAACGAAGCAATGTGGACACACTTGGAGGATTACTTCTGTATATTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtacta        <  2:1561022/150‑1 (MQ=255)
                                                                                                                                             ccGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATTACCTTATtactac       >  2:8913/1‑150 (MQ=255)
                                                                                                                                             ccGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATtactac       <  2:57567/150‑1 (MQ=255)
                                                                                                                                              cGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTaca      >  1:1446627/1‑150 (MQ=255)
                                                                                                                                               ggCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGAGATCACCTTATTACTacaa     >  1:1146676/1‑149 (MQ=255)
                                                                                                                                                  cTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCATAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTaa  <  1:203244/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTCATGGCCGGTGGACAAAATAAATTCTGCTTAGCATAACGACGCAATGTGGACACACTTGGAGGATTACTTCTGTATTTTTCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAA  >  NC_000913/1200763‑1201058

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: