Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,418 T→G 12.0% intergenic (+336/‑36) 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,4180TG12.0% 320.9 / 30.6 124intergenic (+336/‑36)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 (54/55);  new base G (4/11);  total (58/66)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 126 total aligned reads displayed.

TTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGC  >  NC_000913/1426269‑1426566
                                                                                                                                                     |                                                                                                                                                    
ttCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACCCTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAATCATTGCGTCGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACtt                                                                                                                                                      >  2:296245/1‑150 (MQ=255)
  cAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGTCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcg                                                                                                                                                    >  2:1592427/1‑150 (MQ=255)
  cAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcg                                                                                                                                                    >  2:67988/1‑150 (MQ=255)
    ggAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcgcg                                                                                                                                                  <  1:610721/150‑1 (MQ=255)
    ggAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcgcg                                                                                                                                                  >  2:1370569/1‑150 (MQ=255)
     gAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCgg                                                                                                                                                 <  2:1019756/150‑1 (MQ=255)
         aGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTggg                                                                                                                                             <  1:1037848/150‑1 (MQ=255)
         aGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTggg                                                                                                                                             <  1:1478049/150‑1 (MQ=255)
          gtgtgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGt                                                                                                                                            <  1:236814/150‑1 (MQ=255)
           tgtgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtt                                                                                                                                           <  1:1116347/150‑1 (MQ=255)
            gtgtgtTGTTGTTTATTACGTACACTCGGTTACAGGGGTGTGTTGGGGTCTTTTATTAAATAGATGATGCTCCTGGAGAAAAGCATTGCGGGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGttt                                                                                                                                          >  2:1351767/1‑150 (MQ=255)
             tgtgtAGTTGTTATTTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttc                                                                                                                                         <  2:472849/150‑2 (MQ=255)
             tgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttt                                                                                                                                         <  1:778155/150‑1 (MQ=255)
              gtgtAGTTTTTAATTACGTACACTCTGTTACATAGATGAGATTGTGTCTTTAATTAAATCTATGATGCTCCTGGAGAAAAGCATTGCGTGTCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttca                                                                                                                                        <  2:235081/150‑3 (MQ=255)
              gtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttca                                                                                                                                        <  2:1658417/150‑3 (MQ=255)
               tgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcag                                                                                                                                       <  2:1087709/150‑4 (MQ=255)
                ttAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcagg                                                                                                                                      <  1:1611683/149‑5 (MQ=255)
                gtAGTTTTTAATACCGAACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcagg                                                                                                                                      <  1:273250/150‑5 (MQ=255)
                gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttcagg                                                                                                                                      <  1:628068/150‑5 (MQ=255)
                gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttg                                                                                                                                      <  1:529589/150‑1 (MQ=255)
                gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcagg                                                                                                                                      <  1:1717214/150‑5 (MQ=255)
                gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcaggc                                                                                                                                     <  1:720291/150‑6 (MQ=255)
                 tAGTTGTTAATTACGTACACTCTGTTACAGTGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgt                                                                                                                                     <  1:873732/150‑1 (MQ=255)
                 tAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgt                                                                                                                                     >  2:252587/1‑150 (MQ=255)
                 tAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcaggc                                                                                                                                     <  1:1180930/150‑6 (MQ=255)
                 tAGATGTTAATTACGTACACTCTGTTACAGAGATGTGATTGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgt                                                                                                                                     <  1:867812/150‑1 (MQ=255)
                  aGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgtt                                                                                                                                    <  1:143159/150‑1 (MQ=255)
                  aGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgtt                                                                                                                                    <  2:51411/150‑1 (MQ=255)
                      gttAATTACGTACTCTCTGTTACAGAGATGTGATGGTGTCTTTATTTTAATCGATGATGCTCCTTGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAATCCCACTTCGCGTTGGGTTTTTTGTTGTCa                                                                                                                                <  1:1176714/150‑1 (MQ=255)
                       ttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAg                                                                                                                               <  2:1537794/150‑1 (MQ=255)
                        tAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAgg                                                                                                                              <  2:1083570/150‑1 (MQ=255)
                        tAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAgg                                                                                                                              >  2:314166/1‑150 (MQ=255)
                                      ctGTTACAGAGATGTGATGGTGTCTTAAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTAGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACATCGCGGTGGGTTTTTTGTTGTCAGCAGATTTAATAAATg                                                                                                                >  2:1699962/1‑150 (MQ=255)
                                        gTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGc                                                                                                              <  1:67988/150‑1 (MQ=255)
                                            catagaTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCAGGACAAAAGCATTGCGTGGCCTCGTAATCGCTAAATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGCGTTTTTAGTTGTCATGAGTTTTAATAAATGGCatag                                                                                                          >  2:284838/4‑150 (MQ=255)
                                            cAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCagag                                                                                                          >  2:483762/1‑150 (MQ=255)
                                             agagATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACATCGCGGTGGGTTTTTTGTTGTCAGGAGTTTAAATAAATGGCagaga                                                                                                         >  2:332551/1‑149 (MQ=255)
                                              gagaTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCa                                                                                                        >  2:1001660/1‑150 (MQ=255)
                                                gaTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCaaa                                                                                                      >  2:1021679/1‑150 (MQ=255)
                                                  tgtgTTGTTGTCTTTAATTAAATCGATGATGCTACTGGAGAAAAGCATTGCGTGGCCTCGTAATCTCTATATCTTCTATTATGTCGCATGAAACCCACATCGCGGTGGTTTTTTTGTTGGCAGGAGTTTTAATAAATGGCAGAGCAAAcc                                                                                                    <  1:252587/150‑1 (MQ=255)
                                                  tgtgATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAAcc                                                                                                    <  2:922702/150‑1 (MQ=255)
                                                  tgtgATGGTGTATTTAATTAAATCGTTGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGGGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAAcc                                                                                                    <  2:1268774/150‑1 (MQ=255)
                                                      aTGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTcgc                                                                                                <  1:1672333/150‑1 (MQ=255)
                                                         gtgtCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGtc                                                                                             >  1:857701/1‑150 (MQ=255)
                                                            tCTTTAACTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCg                                                                                          >  2:205156/1‑150 (MQ=255)
                                                            gcTTTAATTAAATAGATGATGCTCCTGGAGAAAAGCATTGCGTGGCGTCGTAATCGCTATATCTACTATTATGTCGCCTGACACCCACTTCGCGGTGGGGTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCg                                                                                          <  1:296245/149‑1 (MQ=255)
                                                             cTTTTATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTTTATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGc                                                                                         <  1:65467/150‑1 (MQ=255)
                                                              tttAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCa                                                                                        <  2:644698/150‑1 (MQ=255)
                                                              tttAATTAAAACGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCCTAATCGCTAAATCAACTATTATGTCGCCTGAAACCCAATTCGCGGAGGGTAATTTGTTGTCAGGAGTTTAAATAAATCGCAGAGAAAACCTCGCGTCTCGCa                                                                                        >  2:1368576/1‑150 (MQ=255)
                                                                taattaaATCGATGATGCTCCTGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATATACTATTATGTCGCCTGAAACCCAATTCGCGGAGGGTTTTTTGTTGTAAGGAGTTTTAATAAATGGCAGAGAAAACCTCGAGTCTCACaata                                                                                     >  2:1592175/1‑150 (MQ=255)
                                                                 aactaaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCaata                                                                                     <  1:205156/147‑1 (MQ=255)
                                                                       aTCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAattatt                                                                               <  1:1001660/150‑1 (MQ=255)
                                                                         cGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGa                                                                             >  2:582913/1‑150 (MQ=255)
                                                                          gatgatGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGAt                                                                            >  1:657365/1‑150 (MQ=255)
                                                                             gatgCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGc                                                                         >  2:568045/1‑150 (MQ=255)
                                                                             gatgCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGc                                                                         <  2:804916/150‑1 (MQ=255)
                                                                                gCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCTCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACt                                                                      <  1:647421/150‑1 (MQ=255)
                                                                                    cTGGAGAAAAGCATTGCGTGGCCTCGGAATCGCTTTATCTACTAGTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGGCTGGAGTTTTTATAAATGGCAGAGCACCCCTCGCGTCTCGCAAAAATTATTGATAGCACTggcg                                                                  >  2:332758/1‑148 (MQ=255)
                                                                                     tGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGc                                                                 >  1:1354163/1‑150 (MQ=255)
                                                                                       gagaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCTCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGa                                                               <  1:1547862/150‑1 (MQ=255)
                                                                                        agaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATATTTATTGATAGCACTGGAGCGaa                                                              >  1:385307/1‑150 (MQ=255)
                                                                                         gaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaa                                                             <  2:519815/150‑1 (MQ=255)
                                                                                         gaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaa                                                             <  1:1686208/150‑1 (MQ=255)
                                                                                         gaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaa                                                             >  1:1321372/1‑150 (MQ=255)
                                                                                               cATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTg                                                                              <  2:592632/127‑1 (MQ=255)
                                                                                               cATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTg                                                                              >  1:592632/1‑127 (MQ=255)
                                                                                                aTTGCGTGGTCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATg                                                      <  1:1592427/150‑1 (MQ=255)
                                                                                                   gCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTg                                                   >  1:988447/1‑150 (MQ=255)
                                                                                                    cGTGTCCTCGTAATCGCTATATCTACTATTATGTCTCATGAAACCCACTTCGCGGTTGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGa                                                  <  1:60840/150‑1 (MQ=255)
                                                                                                    cGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGa                                                  >  1:925188/1‑150 (MQ=255)
                                                                                                      tGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACa                                                >  1:1067219/1‑150 (MQ=255)
                                                                                                        gCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAAt                                              >  1:1145106/1‑150 (MQ=255)
                                                                                                        gCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAAt                                              >  1:1567717/1‑150 (MQ=255)
                                                                                                        gCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAAt                                              <  1:535378/150‑1 (MQ=255)
                                                                                                         ccTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATc                                             >  1:317323/1‑150 (MQ=255)
                                                                                                         ccTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATc                                             <  1:1673626/150‑1 (MQ=255)
                                                                                                          cTCGTAATCGCTATATCTACTAATATGTCGCCTGAACCCCACTTCGCGGTGGGTTTTTTGTTGTCAGTAGTTTTAATAAATGGAAGAGCAAACATCGCTTCGCTCAATAATTATTGATAACAATGGAGCGAAAAATAATCCTGCCAATCt                                            >  1:157613/1‑150 (MQ=255)
                                                                                                            cGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGa                                          <  1:777183/150‑1 (MQ=255)
                                                                                                              tAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcc                                        <  1:293945/150‑1 (MQ=255)
                                                                                                               aaTCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcct                                       <  1:1351767/150‑1 (MQ=255)
                                                                                                               aaTCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcct                                       >  2:724789/1‑150 (MQ=255)
                                                                                                                   gCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACctcctc                                   <  2:699464/150‑1 (MQ=255)
                                                                                                                    cTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCa                                  >  1:1320864/1‑150 (MQ=255)
                                                                                                                      atatCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTGTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTAGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTTACATCATCAtt                                >  1:524569/1‑150 (MQ=255)
                                                                                                                         tCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGt                             >  1:1417292/1‑150 (MQ=255)
                                                                                                                         tCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGt                             >  1:482490/1‑150 (MQ=255)
                                                                                                                          ctactaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTc                            >  2:277804/1‑150 (MQ=255)
                                                                                                                            actaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCaa                          >  1:52743/1‑150 (MQ=255)
                                                                                                                             ctaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCaaa                         <  2:1242040/150‑1 (MQ=255)
                                                                                                                              tattatGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCaaaa                        <  2:1567717/150‑1 (MQ=255)
                                                                                                                               attatGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAAt                       <  1:467875/150‑1 (MQ=255)
                                                                                                                                  atGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGAc                    >  1:1590510/1‑150 (MQ=255)
                                                                                                                                   tGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACg                   <  2:1354163/150‑1 (MQ=255)
                                                                                                                                    gTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATATTTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGc                  >  1:1385732/1‑150 (MQ=255)
                                                                                                                                      cGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAATATAATGCTGAAAATCTGACCTCCTCATTAGTCAAAATGACGCAg                >  1:283825/1‑150 (MQ=255)
                                                                                                                                       gCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAgg               >  1:1593228/1‑150 (MQ=255)
                                                                                                                                       gCCTGAAACCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGATGACAATCTGACCTCCTCATTAGTCAAAATGACGCAgga              >  1:1623145/1‑149 (MQ=255)
                                                                                                                                        caaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTACTCATTAGTCAAAATGACGCAGGc              >  2:845169/5‑150 (MQ=255)
                                                                                                                                         aaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCTCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCt             >  2:710948/4‑150 (MQ=255)
                                                                                                                                         aaaaaaCCCACCGCGAAGTGGGGTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGAACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCt             >  1:1186917/4‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGC  >  NC_000913/1426269‑1426566

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: