Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,417 T→C 5.1% intergenic (+335/‑37) ydaY → / → tmpR pseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,4170TC 5.1% 390.9 / 7.8 118intergenic (+335/‑37)ydaY/tmpRpseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base T (64/48);  new base C (3/3);  total (67/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

Only 100 of 118 total aligned reads displayed.

GAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAA  >  NC_000913/1426274‑1426562
                                                                                                                                               |                                                                                                                                                 
gAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcgcg                                                                                                                                              >  3:807139/1‑149 (MQ=255)
  gCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGt                                                                                                                                            <  4:253149/149‑1 (MQ=255)
     gtgtgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTggg                                                                                                                                         >  1:282290/1‑149 (MQ=255)
      tgtgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGt                                                                                                                                        <  1:433437/149‑1 (MQ=255)
           gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttttt                                                                                                                                   >  4:573258/1‑149 (MQ=255)
              gttgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgtt                                                                                                                                >  4:299449/1‑149 (MQ=255)
               ttgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgttg                                                                                                                               >  1:1151054/1‑149 (MQ=255)
               ttgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgttg                                                                                                                               >  4:1120774/1‑149 (MQ=255)
               ttgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttca                                                                                                                                    <  4:172652/144‑3 (MQ=255)
               ttgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttca                                                                                                                                    >  3:172652/1‑142 (MQ=255)
                tgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTtgttgt                                                                                                                              <  1:152905/149‑1 (MQ=255)
                   tAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGAAGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAg                                                                                                                           <  1:850836/149‑1 (MQ=255)
                         cGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGttt                                                                                                                     <  3:745886/149‑1 (MQ=255)
                          gTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGtttt                                                                                                                    >  1:666545/1‑149 (MQ=255)
                          gTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGtttt                                                                                                                    >  4:4562/1‑149 (MQ=255)
                              acTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTtaata                                                                                                                >  1:830999/1‑149 (MQ=255)
                                tctGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATaaa                                                                                                              >  1:1047424/1‑149 (MQ=255)
                                 ctGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAAt                                                                                                             <  2:1151054/149‑1 (MQ=255)
                                 ctGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAAt                                                                                                             <  3:882785/149‑1 (MQ=255)
                                     tACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCa                                                                                                         <  1:514642/149‑1 (MQ=255)
                                     tACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCa                                                                                                         <  3:554607/149‑1 (MQ=255)
                                       cAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCaga                                                                                                       <  1:433174/149‑1 (MQ=255)
                                        agagATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCagag                                                                                                      >  2:221981/1‑149 (MQ=255)
                                         gagaTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAAt                                                                                                             >  4:1067718/1‑141 (MQ=255)
                                         gagaTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAAt                                                                                                             <  3:1067718/141‑1 (MQ=255)
                                             tgtgATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAAc                                                                                                 >  2:749330/1‑149 (MQ=255)
                                                  tGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTcgc                                                                                            >  4:1124064/1‑149 (MQ=255)
                                                    gtgtCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGt                                                                                          >  2:371973/1‑149 (MQ=255)
                                                         tttAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGc                                                                                     >  3:285016/1‑149 (MQ=255)
                                                         tttAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGc                                                                                     >  2:618120/1‑149 (MQ=255)
                                                          ttaattaaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCa                                                                                                                            <  1:855181/109‑1 (MQ=255)
                                                          ttaattaaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCa                                                                                                                            >  2:855181/1‑109 (MQ=255)
                                                              ttaaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCaataa                                                                                >  3:89865/1‑149 (MQ=255)
                                                               taaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCaataat                                                                               <  1:1095945/149‑1 (MQ=255)
                                                                  aTCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAattat                                                                            >  1:16828/1‑149 (MQ=255)
                                                                  aTCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAattat                                                                            <  1:371973/149‑1 (MQ=255)
                                                                    cGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTATTAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTg                                                                          <  1:41613/149‑1 (MQ=255)
                                                                    cGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTg                                                                          <  3:573258/149‑1 (MQ=255)
                                                                     gatgatGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgtt                                                                                                                                >  2:328695/1‑94 (MQ=255)
                                                                     gatgatGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgtt                                                                                                                                <  1:328695/94‑1 (MQ=255)
                                                                      atgatgCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGAt                                                                        <  1:221981/149‑1 (MQ=255)
                                                                       tgatgCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATa                                                                       <  4:89865/149‑1 (MQ=255)
                                                                        gatgCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAg                                                                      >  1:244415/1‑149 (MQ=255)
                                                                        gatgCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAg                                                                      <  3:1120774/149‑1 (MQ=255)
                                                                        gatgCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAg                                                                      >  4:1033710/1‑149 (MQ=255)
                                                                          tgCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCa                                                                    >  2:521489/1‑149 (MQ=255)
                                                                             tCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGc                                                                                                     <  4:558655/113‑1 (MQ=255)
                                                                             tCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGc                                                                                                     >  3:558655/1‑113 (MQ=255)
                                                                             tCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAattatt                                                                           >  2:198731/1‑139 (MQ=255)
                                                                             tCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAg                                                                      >  2:813063/1‑144 (MQ=255)
                                                                             tCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAg                                                                      <  1:813063/144‑1 (MQ=255)
                                                                               cTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTtaa                                                                                                                  <  3:1111867/98‑1 (MQ=255)
                                                                               cTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTtaa                                                                                                                  >  4:1111867/1‑98 (MQ=255)
                                                                                  gagaAAAGCATTGCGTGGCCTCGTAATCGCTCTATCTACTATTATGTCGCCTGCAACCCACTTCGCGGTGAGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCg                                                            >  3:677815/1‑149 (MQ=255)
                                                                                     aaaaGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATaaa                                                                                                              <  4:1143652/96‑1 (MQ=255)
                                                                                     aaaaGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATaaa                                                                                                              >  3:1143652/1‑96 (MQ=255)
                                                                                       aaGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCaga                                                                                                       >  4:406864/1‑101 (MQ=255)
                                                                                       aaGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCaga                                                                                                       <  3:406864/101‑1 (MQ=255)
                                                                                        aGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAaat                                                      <  1:749330/149‑1 (MQ=255)
                                                                                         gCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATgg                                                                                                           >  1:698302/1‑95 (MQ=255)
                                                                                         gCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAaata                                                     >  2:592499/1‑149 (MQ=255)
                                                                                           aTTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAaataat                                                   <  4:811444/149‑1 (MQ=255)
                                                                                              gCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCGAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCt                                                <  3:299449/149‑1 (MQ=255)
                                                                                                gTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGa                                              <  4:285016/149‑1 (MQ=255)
                                                                                                   gCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACaa                                           >  2:488309/1‑149 (MQ=255)
                                                                                                    ccTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAAt                                          <  4:359506/149‑1 (MQ=255)
                                                                                                      tCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCt                                        <  2:666545/149‑1 (MQ=255)
                                                                                                        gTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGa                                      <  3:4562/149‑1 (MQ=255)
                                                                                                          aaTCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAaataat                                                   >  2:818721/1‑134 (MQ=255)
                                                                                                            tCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcctc                                  <  4:849729/149‑1 (MQ=255)
                                                                                                              gCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGc                                                             >  1:300547/1‑120 (MQ=255)
                                                                                                              gCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGc                                                             <  2:300547/120‑1 (MQ=255)
                                                                                                               cTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGa                                                                                                                         >  4:63036/1‑59 (MQ=255)
                                                                                                                 atatCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCAt                             <  3:1033087/149‑1 (MQ=255)
                                                                                                                 atatCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCAt                             <  4:639596/149‑1 (MQ=255)
                                                                                                                 atatCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCAt                             <  1:9084/149‑1 (MQ=255)
                                                                                                                  tatCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCAtt                            >  3:218554/1‑149 (MQ=255)
                                                                                                                    tCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAg                          >  2:131924/1‑149 (MQ=255)
                                                                                                                    tCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTATTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGTGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAg                          >  3:745491/1‑149 (MQ=255)
                                                                                                                     ctactaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGt                         >  1:571344/1‑149 (MQ=255)
                                                                                                                      tactaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTc                        >  3:1004756/1‑149 (MQ=255)
                                                                                                                           ttatGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaaaa                                                       >  2:337026/1‑113 (MQ=255)
                                                                                                                           ttatGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaaaa                                                       <  1:337026/113‑1 (MQ=255)
                                                                                                                           ttatGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcc                                    >  4:980376/1‑132 (MQ=255)
                                                                                                                           ttatGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcc                                    <  3:980376/132‑1 (MQ=255)
                                                                                                                             atGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGa                 <  2:1047424/149‑1 (MQ=255)
                                                                                                                                tCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGc              >  2:817795/1‑149 (MQ=255)
                                                                                                                                tCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGc              >  3:784177/1‑149 (MQ=255)
                                                                                                                                    aaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGc          >  3:859529/4‑149 (MQ=255)
                                                                                                                                     tGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCt         >  2:765134/1‑149 (MQ=255)
                                                                                                                                     tGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCt         >  1:1147893/1‑149 (MQ=255)
                                                                                                                                     aaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCt         >  3:721678/3‑149 (MQ=255)
                                                                                                                                         aCCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTgggg     >  1:788182/1‑149 (MQ=255)
                                                                                                                                          cccACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGa    >  3:332086/1‑149 (MQ=255)
                                                                                                                                           ccACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGaa   >  2:418215/1‑149 (MQ=255)
                                                                                                                                           ccACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGaa   >  2:251920/1‑149 (MQ=255)
                                                                                                                                           ccACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGaa   >  4:930759/1‑149 (MQ=255)
                                                                                                                                           ccACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGaa   >  1:851244/1‑149 (MQ=255)
                                                                                                                                           ccACTTCGCGGGGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGaa   >  3:918632/1‑149 (MQ=255)
                                                                                                                                            cACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGaaa  >  2:354519/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAA  >  NC_000913/1426274‑1426562

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