Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 570,284 C→T I461I (ATC→ATT ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913570,2840CT98.6% 494.7 / ‑2.3 141I461I (ATC→ATTybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (59/80);  total (60/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 141 total aligned reads displayed.

AAGGGCTTCGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATCATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAATATATGGTTGATGGTGAGATT  >  NC_000913/570149‑570419
                                                                                                                                        |                                                                                                                                       
aaGGGCTTCGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTAtaa                                                                                                                                     >  2:1648094/1‑140 (MQ=255)
  gggCTTCGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTAtaatt                                                                                                                                   >  1:989800/1‑140 (MQ=255)
   ggCTTCGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTAtaatta                                                                                                                                  >  2:268135/1‑140 (MQ=255)
   ggCTTCGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTAtaatta                                                                                                                                  >  1:444495/1‑140 (MQ=255)
    gCTTCGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTTAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTAtaattaa                                                                                                                                 <  1:1750830/140‑1 (MQ=255)
     cTTCGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTAtaattaat                                                                                                                                <  2:1290600/140‑1 (MQ=255)
     cTTCGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTAtaattaat                                                                                                                                <  2:1755998/140‑1 (MQ=255)
      ttCGGTATTCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTAtaattaat                                                                                                                                <  1:1179595/140‑1 (MQ=255)
      ttCGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATa                                                                                                                               >  1:412754/1‑140 (MQ=255)
        cGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACg                                                                                                                             >  2:1241339/1‑140 (MQ=255)
         ggTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTATACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACgg                                                                                                                            <  2:177710/140‑1 (MQ=255)
         ggTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACgg                                                                                                                            <  2:975817/140‑1 (MQ=255)
         ggTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACgg                                                                                                                            <  2:676693/140‑1 (MQ=255)
         ggTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACgg                                                                                                                            <  1:1822181/140‑1 (MQ=255)
              tCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGataa                                                                                                                        >  2:432398/1‑140 (MQ=255)
              tCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGataa                                                                                                                        >  1:1512244/1‑140 (MQ=255)
               cATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGataat                                                                                                                       >  1:1461584/1‑140 (MQ=255)
               cATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGataat                                                                                                                       <  2:919864/140‑1 (MQ=255)
               cATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGataat                                                                                                                       <  1:1309163/140‑1 (MQ=255)
                 ttGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGataat                                                                                                                       >  1:726576/1‑138 (MQ=255)
                 ttGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGataat                                                                                                                       <  2:726576/138‑1 (MQ=255)
                    aaaaCTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAAc                                                                                                                  <  2:1444399/140‑1 (MQ=255)
                      aaCTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCt                                                                                                                <  1:268265/140‑1 (MQ=255)
                      aaCTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCt                                                                                                                <  1:33758/140‑1 (MQ=255)
                      aaCTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCt                                                                                                                <  1:1241339/140‑1 (MQ=255)
                        cTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCtgt                                                                                                              <  2:580507/140‑1 (MQ=255)
                          ttGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGa                                                                                                            >  2:1765222/1‑140 (MQ=255)
                            gAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGata                                                                                                          <  1:119018/140‑1 (MQ=255)
                                  tAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTAtt                                                                                                    >  1:1624409/1‑140 (MQ=255)
                                    aaGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTAtttt                                                                                                  <  2:729637/140‑1 (MQ=255)
                                        tgtAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGa                                                                                              <  1:1011418/140‑1 (MQ=255)
                                         gtAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGaa                                                                                             <  2:463894/140‑1 (MQ=255)
                                          tAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAAt                                                                                            >  2:337741/1‑140 (MQ=255)
                                           aaaTTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGatat                                                                                                         <  2:935407/126‑1 (MQ=255)
                                           aaaTTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGatat                                                                                                         >  1:935407/1‑126 (MQ=255)
                                            aaTTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATgg                                                                                          >  1:1331647/1‑140 (MQ=255)
                                               tcattcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCAt                                                                                       <  1:193875/140‑1 (MQ=255)
                                               tcattcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCAt                                                                                       <  1:1696902/140‑1 (MQ=255)
                                               tcattcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCAt                                                                                       <  1:947574/140‑1 (MQ=255)
                                                 attcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTg                                                                                     <  1:1505584/140‑1 (MQ=255)
                                                 attcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTg                                                                                     <  2:517796/140‑1 (MQ=255)
                                                 attcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTg                                                                                     <  2:824048/140‑1 (MQ=255)
                                                 attcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTg                                                                                     <  1:1037199/140‑1 (MQ=255)
                                                 attcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAAACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGt                                                                                    >  1:1251552/1‑140 (MQ=255)
                                                  ttcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTg                                                                                     <  1:617482/139‑1 (MQ=255)
                                                  ttcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGt                                                                                    >  1:30618/1‑140 (MQ=255)
                                                  ttcaGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGt                                                                                    >  1:1506882/1‑140 (MQ=255)
                                                    caGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGttt                                                                                  >  1:137698/1‑140 (MQ=255)
                                                      gTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGttttt                                                                                <  2:874605/140‑1 (MQ=255)
                                                      gTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGttttt                                                                                <  2:1691239/140‑1 (MQ=255)
                                                       tGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTa                                                                               >  1:728913/1‑140 (MQ=255)
                                                           ttttATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAaca                                                                           <  2:1512244/140‑1 (MQ=255)
                                                             ttATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACt                                                                         <  2:1840592/140‑1 (MQ=255)
                                                             ttATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACt                                                                         <  1:1784961/140‑1 (MQ=255)
                                                                 gCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATcc                                                                     >  1:1651605/1‑140 (MQ=255)
                                                                 gCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATcc                                                                     <  2:103135/140‑1 (MQ=255)
                                                                 gCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATcc                                                                     >  1:1244195/1‑140 (MQ=255)
                                                                     ttGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTa                                                                 >  2:1386708/1‑140 (MQ=255)
                                                                     ttGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTa                                                                 >  2:1389466/1‑140 (MQ=255)
                                                                       gACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAAt                                                               <  1:568492/140‑1 (MQ=255)
                                                                        aCTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATg                                                              <  2:1604217/140‑1 (MQ=255)
                                                                        aCTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATg                                                              <  1:827114/140‑1 (MQ=255)
                                                                        aCTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATg                                                              <  1:1367646/140‑1 (MQ=255)
                                                                         cTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGa                                                             >  2:612564/1‑140 (MQ=255)
                                                                         cTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGa                                                             <  2:655253/140‑1 (MQ=255)
                                                                             aCCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGt                                                         >  1:16020/1‑140 (MQ=255)
                                                                                  tAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAatat                                                    >  2:1043403/1‑140 (MQ=255)
                                                                                  tAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAatat                                                    >  1:631216/1‑140 (MQ=255)
                                                                                    aaaaTGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATcc                                                  >  1:208309/1‑140 (MQ=255)
                                                                                     aaaTGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCg                                                 <  2:1201611/140‑1 (MQ=255)
                                                                                     aaaTGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCg                                                 <  2:819824/140‑1 (MQ=255)
                                                                                          gACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCcag                                            <  2:1190283/140‑1 (MQ=255)
                                                                                          gACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCcag                                            <  2:444495/140‑1 (MQ=255)
                                                                                          gACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGACAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCcag                                            <  1:1005429/140‑1 (MQ=255)
                                                                                                  ttATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGAt                                    <  1:1716041/140‑1 (MQ=255)
                                                                                                  ttATGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGAt                                    <  2:535104/140‑1 (MQ=255)
                                                                                                    aTGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATaa                                  <  2:631216/140‑1 (MQ=255)
                                                                                                    aTGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATaa                                  >  1:419236/1‑140 (MQ=255)
                                                                                                     tGCCGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAg                                 >  1:1459909/1‑140 (MQ=255)
                                                                                                       ccGTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTg                               <  1:1389466/140‑1 (MQ=255)
                                                                                                         gTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGct                             <  2:254795/140‑1 (MQ=255)
                                                                                                         gTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGct                             <  2:30618/140‑1 (MQ=255)
                                                                                                         gTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGct                             <  1:80965/140‑1 (MQ=255)
                                                                                                         gTACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGct                             <  2:1506882/140‑1 (MQ=255)
                                                                                                          tACACTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGctc                            >  2:114299/1‑140 (MQ=255)
                                                                                                           tcacTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGctct                           <  2:1430056/139‑1 (MQ=255)
                                                                                                            cacTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCtctc                          <  1:1768341/140‑1 (MQ=255)
                                                                                                            cacTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCtctc                          <  2:1244195/140‑1 (MQ=255)
                                                                                                              cTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCaa                        <  2:412754/140‑1 (MQ=255)
                                                                                                              cTTGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCaa                        <  2:400519/140‑1 (MQ=255)
                                                                                                                tGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAat                      >  2:1119434/1‑140 (MQ=255)
                                                                                                                tGCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGCGATATCTATTTTATGAATGGCATTGTTTTTAAATACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAat                      <  1:432398/140‑1 (MQ=255)
                                                                                                                 gCCTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAata                     >  2:1800979/1‑140 (MQ=255)
                                                                                                                  ccTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAatat                    >  2:789153/1‑140 (MQ=255)
                                                                                                                   cTATAAAACATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAatata                   <  2:329367/140‑1 (MQ=255)
                                                                                                                          aCATTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAATATATGGTTGa            <  1:1648094/140‑1 (MQ=255)
                                                                                                                            aTTCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAATATATGGTTGATg          >  2:1746519/1‑140 (MQ=255)
                                                                                                                              tCGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAATATATGGTTGATGGt        >  2:1782299/1‑140 (MQ=255)
                                                                                                                               cGAAAAAATTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAATATATGGTTGATGGTg       <  2:1852806/140‑1 (MQ=255)
                                                                                                                                    aaaTTATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAATATATGGTTGATGGTGAGatt  <  2:1145632/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AAGGGCTTCGGTATCATTGAAAACTTTGAAGAGTAAAGGTGTAAATTCATTCAGTGATTTTTATGCTATTGACTTAACCAGTAAAAATGGACGAGAGTTATGCCGTACACTTGCCTATAAAACATTCGAAAAAATCATAATTAATACGGATAATAAAACCTGTGATATCTATTTTATGAATGGCATTGTTTTTAAACACTATCCTTTAATGAAAGTAATATCCGCCCAGCAGGCGATAAGTGCTCTCAAATATATGGTTGATGGTGAGATT  >  NC_000913/570149‑570419

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