Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 394,473 G→T 8.3% intergenic (+55/‑33) yaiT → / → yaiV pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein/putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913394,4730GT 8.3% 296.8 / 14.0 108intergenic (+55/‑33)yaiT/yaiVpseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein/putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (59/40);  new base T (5/4);  total (64/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 113 total aligned reads displayed.

CTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTC  >  NC_000913/394328‑394622
                                                                                                                                                 |                                                                                                                                                     
cTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACTTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGAAATATTATTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGAttt                                                                                                                                                   >  1:1450195/1‑150 (MQ=255)
cTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGAttt                                                                                                                                                   >  1:1117172/1‑150 (MQ=255)
cTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGAttt                                                                                                                                                   >  1:1132525/1‑150 (MQ=255)
    aGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACg                                                                                                                                               <  2:1099940/150‑1 (MQ=255)
     gcgcCTATACCGATGCTAACTACCTCTTTGGTGGTGACGTAGATCAAGACTGGTCCGCGTACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTACGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTTACgg                                                                                                                                              <  1:560660/150‑1 (MQ=255)
      cgcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACggg                                                                                                                                             <  1:671484/150‑1 (MQ=255)
      cgcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAATCCTTGACggg                                                                                                                                             <  1:1578966/150‑1 (MQ=255)
       gcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTTTTAAATATTCATGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGc                                                                                                                                            <  2:1223566/150‑1 (MQ=255)
       gcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAATCCTTGACGGGc                                                                                                                                            <  2:405779/150‑1 (MQ=255)
       gcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAATCCTTGACGGGc                                                                                                                                            <  1:1410396/150‑1 (MQ=255)
        ccTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCa                                                                                                                                           <  2:730520/150‑1 (MQ=255)
        ccTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCa                                                                                                                                           <  1:52130/150‑1 (MQ=255)
          tataCCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCAtt                                                                                                                                         >  2:493702/1‑150 (MQ=255)
             aCCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACt                                                                                                                                      <  1:1246701/150‑1 (MQ=255)
             aCCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACt                                                                                                                                      >  1:1449603/1‑150 (MQ=255)
              ccGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTg                                                                                                                                     >  1:1358184/1‑150 (MQ=255)
                gATGCTATCTCCCTCGGTGGGTGGGACGTAGATCAAGACTGGTCCGCGAACGGGGGTGTTAAATATACCTGGTAATATTCTTCACTCAGAATAAATACTGTTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCa                                                                                                                                   <  2:1128774/150‑1 (MQ=255)
                gATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCa                                                                                                                                   <  1:396363/150‑1 (MQ=255)
                 aTGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCaa                                                                                                                                  >  1:213074/1‑150 (MQ=255)
                    cTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAATGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGAAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAgg                                                                                                                               >  1:1270471/1‑150 (MQ=255)
                    cTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAgg                                                                                                                               >  1:693183/1‑150 (MQ=255)
                     tAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGa                                                                                                                              >  2:535443/1‑150 (MQ=255)
                      aaCTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGAc                                                                                                                             >  2:674927/1‑150 (MQ=255)
                       aCTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGAcg                                                                                                                            >  1:1677690/1‑150 (MQ=255)
                        cTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGAcgc                                                                                                                           <  1:1777568/150‑1 (MQ=255)
                          aCCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACgcgc                                                                                                                         >  1:476155/1‑150 (MQ=255)
                                   gtggttACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTAt                                                                                                                <  2:1388952/150‑1 (MQ=255)
                                       tgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATTCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGt                                                                                                            >  2:1654380/1‑150 (MQ=255)
                                        gACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAg                                                                                                          >  1:1510484/1‑150 (MQ=255)
                                          catAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                         >  1:800961/3‑150 (MQ=255)
                                          cGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                         >  2:946800/1‑150 (MQ=255)
                                          cGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                         <  2:1193890/150‑1 (MQ=255)
                                            tAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCATAAGTTATCTGTAGTTa                                                                                                       >  2:1688551/1‑150 (MQ=255)
                                            tAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTa                                                                                                       >  1:1861191/1‑150 (MQ=255)
                                              taTCAAGACTGGGCAGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTaaa                                                                                                     <  1:1635069/149‑1 (MQ=255)
                                               aTCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAAc                                                                                                    <  2:1862023/150‑1 (MQ=255)
                                                 cAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACct                                                                                                  >  2:572564/1‑150 (MQ=255)
                                                  aaGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACctc                                                                                                 >  1:447267/1‑150 (MQ=255)
                                                   aGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACctct                                                                                                <  1:459147/150‑1 (MQ=255)
                                                        ggTCCGCGAACGTGGGTGGTATATATCCCTTGTAATATTCTTCACTCCGAAGAAATACTGTTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAgg                                                                                           <  2:172817/150‑1 (MQ=255)
                                                          tCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGaa                                                                                         >  2:503478/1‑150 (MQ=255)
                                                           ccGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCAATCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCAGTCAAGGATTTGACGGGCATTACTGCAAAGGCCGCGCAAATGTTATCTGTAGTTAAACCTATTCAGGAAt                                                                                        >  2:879920/1‑150 (MQ=255)
                                                            cgcgAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAAtt                                                                                       >  2:794505/1‑150 (MQ=255)
                                                               gAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTgg                                                                                    <  2:63293/150‑1 (MQ=255)
                                                                    tGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGc                                                                               <  1:420841/150‑1 (MQ=255)
                                                                     gggTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCt                                                                              >  1:879767/1‑150 (MQ=255)
                                                                      ggTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTc                                                                             <  1:555805/150‑1 (MQ=255)
                                                                      ggTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTc                                                                             <  2:1042291/150‑1 (MQ=255)
                                                                       gtgtTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCg                                                                            >  1:1249155/1‑150 (MQ=255)
                                                                       gtgtTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCg                                                                            <  2:800961/150‑1 (MQ=255)
                                                                          ttAAATATACCTGGTAATATTATTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGGCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATa                                                                         <  1:1375463/150‑1 (MQ=255)
                                                                               tataCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGt                                                                    >  1:619900/1‑150 (MQ=255)
                                                                                 taCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGttt                                                                  >  1:1598113/1‑150 (MQ=255)
                                                                                  aCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATCTGGTAAGCTCGATAAATGTTTg                                                                 <  1:228194/150‑1 (MQ=255)
                                                                                    cTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTc                                                               >  2:1470980/1‑150 (MQ=255)
                                                                                     tGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATATAAATAATGCCCTTCAATGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTcc                                                              <  1:1613993/150‑1 (MQ=255)
                                                                                      ggTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTcaa                                                             >  1:7541/1‑148 (MQ=255)
                                                                                            atTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACg                                                       >  1:1729009/1‑150 (MQ=255)
                                                                                             ttcttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACgg                                                      >  1:1289344/1‑150 (MQ=255)
                                                                                               cttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGTACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGATCGATAAATGTTTGTCCAGATACGGTa                                                    >  1:1175076/1‑150 (MQ=255)
                                                                                               catcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTa                                                    <  2:306480/148‑1 (MQ=255)
                                                                                                ttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTAc                                                   >  2:1537101/1‑150 (MQ=255)
                                                                                                 tcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTAcg                                                  >  2:1256533/1‑150 (MQ=255)
                                                                                                       cGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTc                                            <  1:1251101/150‑1 (MQ=255)
                                                                                                        gaagaaATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCTAATGTTATGTGTAGTTAATCCTCTTCAGGAATTTGGTAAGCTCGATATATGTTTGTCGAGATACGGTAcgagggtcg                                           >  2:232476/1‑143 (MQ=255)
                                                                                                         aagaaATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGa                                          >  2:1861064/1‑150 (MQ=255)
                                                                                                              aTACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTta                                     <  1:1704978/150‑1 (MQ=255)
                                                                                                              aTACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTta                                     <  2:1677690/150‑1 (MQ=255)
                                                                                                               tACTGGTAATTTAATCTAAATAATGCCCGTCATGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTtaa                                    >  1:1465988/1‑150 (MQ=255)
                                                                                                               tACTGGTAATTTAATCTAAATAATGCCAGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTtaa                                    <  1:106653/150‑1 (MQ=255)
                                                                                                                  tGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTtaataa                                 >  2:709405/1‑150 (MQ=255)
                                                                                                                   ggTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTaataat                                >  2:112777/1‑150 (MQ=255)
                                                                                                                    gTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATg                               >  2:322929/1‑150 (MQ=255)
                                                                                                                     tAATTTAATCTAAATAATGCCCTTAAAGGATTTGACGGGAATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGa                              <  2:1365034/150‑1 (MQ=255)
                                                                                                                      aaTTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGaa                             <  1:950892/150‑1 (MQ=255)
                                                                                                                        tttAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGaaaa                           <  1:1460095/150‑1 (MQ=255)
                                                                                                                         ttAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAg                          <  1:1871916/150‑1 (MQ=255)
                                                                                                                         ttAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAg                          <  1:1864752/150‑1 (MQ=255)
                                                                                                                           aaTCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCa                        >  2:644112/1‑150 (MQ=255)
                                                                                                                              cTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGt                     >  2:1136429/1‑150 (MQ=255)
                                                                                                                               tAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGtt                    >  1:1793874/1‑150 (MQ=255)
                                                                                                                               tAAATAATGCCAGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGtt                    <  2:1541234/150‑1 (MQ=255)
                                                                                                                                cagtaatGCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTta                   >  2:1867693/4‑150 (MQ=255)
                                                                                                                                 aataatGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTtat                  >  2:306778/1‑150 (MQ=255)
                                                                                                                                 aataatGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTtat                  >  2:1727770/1‑150 (MQ=255)
                                                                                                                                    aatGCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAtt               >  1:1312051/1‑150 (MQ=255)
                                                                                                                                     atGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAttt              >  2:1320740/1‑150 (MQ=255)
                                                                                                                                      tGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAtttt             >  1:1334052/1‑150 (MQ=255)
                                                                                                                                      tGCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAtttt             >  1:1796750/1‑150 (MQ=255)
                                                                                                                                       gCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTc            >  1:1014327/1‑150 (MQ=255)
                                                                                                                                       gCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTc            >  1:79710/1‑150 (MQ=255)
                                                                                                                                       gCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAACCAAGTTATATTTTc            >  2:784578/1‑150 (MQ=255)
                                                                                                                                        cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc           >  2:558702/1‑150 (MQ=255)
                                                                                                                                            tCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTg       <  1:794505/150‑1 (MQ=255)
                                                                                                                                               tggATTTGACTGGCATTACTGCAAATGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGATTTTTGTAAGATCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATg    <  1:1421661/149‑1 (MQ=255)
                                                                                                                                               aGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATg    <  1:946800/150‑1 (MQ=255)
                                                                                                                                                ggATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTCCGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGt   >  1:1146821/1‑150 (MQ=255)
                                                                                                                                                ggATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAAAGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAATAGCAAGTTATATTTTCAAGTGATGt   >  1:1152793/1‑150 (MQ=255)
                                                                                                                                                 gATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGTAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTc  <  1:1178864/150‑1 (MQ=255)
                                                                                                                                                 gATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTc  <  2:1358184/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTC  >  NC_000913/394328‑394622

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: