Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 394,473 G→T 5.0% 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 5.0% 381.7 / 5.9 134intergenic (+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 (61/65);  new base T (6/2);  total (67/67)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

GAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATG  >  NC_000913/394325‑394620
                                                                                                                                                    |                                                                                                                                                   
gAACTTCAGCGCCTATACCGATGCTAACTACGTCGGTGTTGGTGACGTATATCAAGACTGGTACGCGAACGTGTGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGa                                                                                                                                                    <  2:845567/150‑1 (MQ=255)
    ttCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTg                                                                                                                                                <  2:1858399/150‑1 (MQ=255)
    ttCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTg                                                                                                                                                >  1:1356577/1‑150 (MQ=255)
      cAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGAc                                                                                                                                              <  2:1653990/150‑1 (MQ=255)
        gcgcCTATACCGATGCTAACTACATCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACgg                                                                                                                                            <  1:1052486/150‑1 (MQ=255)
         cgcCTATATCGATGCTATCTACCTCGCTGGTGGTGACTTAGATAAAGTCTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACggg                                                                                                                                           <  2:190503/150‑1 (MQ=255)
         cgcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACggg                                                                                                                                           <  2:1179300/150‑1 (MQ=255)
          gcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTTGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGc                                                                                                                                          <  2:1284153/150‑1 (MQ=255)
          gcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGc                                                                                                                                          >  1:491231/1‑150 (MQ=255)
          gcCTATACCGAAGCTAACTACCTCGGTGGTTGTGACGTAGATCAAGACTGGTCAGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGc                                                                                                                                          <  1:76358/150‑1 (MQ=255)
           ccTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCa                                                                                                                                         >  2:1293774/1‑150 (MQ=255)
           ccTAAACCGATGCTAACTACCTCGGTGGTGTTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCAAAATAATGCCCGTCAAATCCTTGACGGGCa                                                                                                                                         <  2:754621/150‑1 (MQ=255)
             tataCCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCAtt                                                                                                                                       <  2:1881259/150‑1 (MQ=255)
               taCCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTAc                                                                                                                                     >  2:1375884/1‑150 (MQ=255)
               taCCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTAc                                                                                                                                     >  1:2063182/1‑150 (MQ=255)
               taCCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTGCACTCCGAAGAAATCCTGGTAATTTAATCTAAATAATGCCCGTCAAGGAGTTGCCGGGCAttta                                                                                                                                     >  1:1006922/1‑148 (MQ=255)
               taCCGATGCAAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTAc                                                                                                                                     <  1:1290523/150‑1 (MQ=255)
                 ccGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGTTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTg                                                                                                                                   <  1:1348370/150‑1 (MQ=255)
                 ccGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTg                                                                                                                                   >  2:1537871/1‑150 (MQ=255)
                   gATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCa                                                                                                                                 >  1:280605/1‑150 (MQ=255)
                      gCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCTAAg                                                                                                                              >  2:868943/1‑150 (MQ=255)
                      gCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAg                                                                                                                              <  2:543853/150‑1 (MQ=255)
                       cTAACTACCTCGGTGGTGGTGACGTAGATCAATACTGGTCAGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATTCCCGTCAAGGATTTGAAGGTCATTACTTCAAAgg                                                                                                                             >  2:201351/1‑150 (MQ=255)
                       cTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAgg                                                                                                                             <  2:914534/150‑1 (MQ=255)
                       cTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAgg                                                                                                                             >  1:538345/1‑150 (MQ=255)
                            tACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGAcgcg                                                                                                                        <  2:111578/150‑1 (MQ=255)
                            tACCTCGGAGGTGGGGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGAcgcg                                                                                                                        <  1:1690182/150‑1 (MQ=255)
                               cTCGGTGGTGGTGACGTAGATAAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGATGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCaa                                                                                                                     >  2:1467821/1‑150 (MQ=255)
                                  ggtggtggtGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATg                                                                                                                  <  1:408231/150‑1 (MQ=255)
                                  ggtgggggtGACGTAGATCATGAATGGTCCGCGAACGTTGGTGTTAAATATACCTGTTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATg                                                                                                                  <  1:1794654/150‑1 (MQ=255)
                                    tggtggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGtt                                                                                                                <  2:1108891/150‑1 (MQ=255)
                                     ggtggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTa                                                                                                               >  1:809740/1‑150 (MQ=255)
                                      gtggagACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTAt                                                                                                              <  2:356319/150‑1 (MQ=255)
                                        ggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCt                                                                                                            >  2:639848/1‑150 (MQ=255)
                                        ggtgACGTAGATCAAGACTGGTACGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGGCGCGCAAATGTTATcg                                                                                                            >  1:798466/1‑149 (MQ=255)
                                         gtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTg                                                                                                           >  2:219549/1‑150 (MQ=255)
                                         gtgACGTAGATCAAGACTGATCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTg                                                                                                           <  2:2037370/150‑1 (MQ=255)
                                          tgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGt                                                                                                          >  2:1906581/1‑150 (MQ=255)
                                             cGTAGATCAAGACTGGTCCGCGATCGTGGGTGTTAATTATACCTGGTAATATTCTTCACTCAGAAGAAATACTGGTAATTTAATCTAAATAATGCCAGTCAGGTATTTGACGGGAATTACTTCAAAGGACGGGCTAATGTTATCTTTAGt                                                                                                       >  1:331113/1‑150 (MQ=255)
                                             cGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGt                                                                                                       <  2:827280/150‑1 (MQ=255)
                                                 gATCAAGACTGGTCCGCGAACGTGGTTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTaaa                                                                                                   <  1:1265626/150‑1 (MQ=255)
                                                 gATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTaaa                                                                                                   <  2:488393/150‑1 (MQ=255)
                                                  aTCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAAc                                                                                                  >  1:1065448/1‑150 (MQ=255)
                                                  aTCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAAc                                                                                                  >  1:1808190/1‑150 (MQ=255)
                                                    cAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAAGCTAAATAATGCCCGTCATGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACct                                                                                                >  2:826877/1‑150 (MQ=255)
                                                      aGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTATTCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACctct                                                                                              >  2:805343/1‑150 (MQ=255)
                                                       gtgtGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGTCGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCtt                                                                                             <  1:1464303/147‑1 (MQ=255)
                                                         cTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCa                                                                                           >  1:616555/1‑150 (MQ=255)
                                                         cTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCa                                                                                           <  1:416184/150‑1 (MQ=255)
                                                           ggTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAgg                                                                                         >  1:978277/1‑150 (MQ=255)
                                                                 cgAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTg                                                                                   >  2:786394/1‑150 (MQ=255)
                                                                    aCGTGGTTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGTCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTa                                                                                <  1:1382062/150‑1 (MQ=255)
                                                                     cGTGGGTTTTAAAAATACCTGGTAATATTCTTCACTCCTAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGTTTTGACGGGAATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTaa                                                                               <  2:1762392/150‑1 (MQ=255)
                                                                       tGGTTGTTAAATATAGCTGGAATTTTTCTTCATTCCTATGAATTACTGGAATTTTAATCTAATTATTGCCCGTCATGGATTTGAGGGGGATTACTGCAAAGGACGCGCTAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGc                                                                             <  1:1596741/150‑1 (MQ=255)
                                                                       ggggTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGc                                                                             <  2:1964844/149‑1 (MQ=255)
                                                                         ggTGTTAAATATACCTGTTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTc                                                                           <  2:775367/150‑1 (MQ=255)
                                                                         ggTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTc                                                                           <  1:1697618/150‑1 (MQ=255)
                                                                          gtgtTAAATATACCTGGTAATATTGTTCTCTCCGAAGAAATTCTGTTAATTTAATCTATATAATGCCCGTCATGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTTAACCTCTTCAGGAATTTGGTAAGCTCg                                                                          <  2:1564516/150‑1 (MQ=255)
                                                                            gtTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGTTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAACCCACTAAAGGAATTTGGTAAGCTCGAt                                                                        >  2:1623477/1‑150 (MQ=255)
                                                                                 atatACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGTACGCGCAAATGTTCTCTGTAGTTACACCTCTTCAGGAATTTGGTCAGCTCGATAAATg                                                                   >  1:336783/1‑150 (MQ=255)
                                                                                 aaataCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGACCTTGAAGTATTTGACGGGCATTACTGAAAAGGACTGGAATATTTTATCTGTAGTTACAAGTCTTCAGGAATTTGGTAAGCTCGATAAAtt                                                                   >  1:341193/3‑149 (MQ=255)
                                                                                     aacTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTATTCTAAATAATGCCCGTCTAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTg                                                               <  1:721404/148‑1 (MQ=255)
                                                                                     aCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTg                                                               >  1:1715478/1‑150 (MQ=255)
                                                                                      ccTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGt                                                              <  2:434424/150‑1 (MQ=255)
                                                                                       cTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTc                                                             >  2:1505235/1‑150 (MQ=255)
                                                                                        tGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAATCCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTcc                                                            <  1:266581/150‑1 (MQ=255)
                                                                                          gTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAg                                                          <  2:1644895/150‑1 (MQ=255)
                                                                                              tatTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATAc                                                      <  1:1911639/150‑1 (MQ=255)
                                                                                                ttcttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACgg                                                    >  1:2001046/1‑150 (MQ=255)
                                                                                                ttcttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACgg                                                    <  2:289340/150‑1 (MQ=255)
                                                                                                 tcttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGt                                                   >  2:1929068/1‑150 (MQ=255)
                                                                                                  cttcACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTa                                                  >  1:841130/1‑150 (MQ=255)
                                                                                                     cACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTAcgc                                               <  1:955481/150‑1 (MQ=255)
                                                                                                      aCTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTAAAACCTCTTCAGTAATTTGGTAAGCTCGATAAATGTTTGTCCACATACGGTAcgcg                                              >  2:1756128/1‑150 (MQ=255)
                                                                                                       cTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACgcgc                                             >  2:365002/1‑150 (MQ=255)
                                                                                                          cGAAGAAATACTGGTAATTTAATCTAAATAATGCCCTTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTc                                          <  2:1578828/150‑1 (MQ=255)
                                                                                                              gaaATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGt                                      >  2:1319929/1‑150 (MQ=255)
                                                                                                                aaTACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGttt                                    >  1:785275/1‑150 (MQ=255)
                                                                                                                 aTACTGGTAATTTAATCTAAATAATTCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTta                                   >  2:594451/1‑150 (MQ=255)
                                                                                                                  tACTGGTAATTAAATCTAAATAATGCCCTTCAATGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTtaa                                  <  2:855629/150‑1 (MQ=255)
                                                                                                                          aTTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGaaa                          <  1:1573564/150‑1 (MQ=255)
                                                                                                                              aaTCTAAATAATGCCCGTCAAGGATTTGTCGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCa                      >  1:934034/1‑150 (MQ=255)
                                                                                                                               aTCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCaa                     >  1:1634779/1‑150 (MQ=255)
                                                                                                                                 cTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGt                   >  2:387795/1‑150 (MQ=255)
                                                                                                                                   cagtaatGCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTta                 >  1:2074846/4‑150 (MQ=255)
                                                                                                                                   aaaTAATGCCCGTCAAGTATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTta                 <  1:1318037/150‑1 (MQ=255)
                                                                                                                                      taatGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCTATTGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATATATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTatat              >  1:447467/1‑150 (MQ=255)
                                                                                                                                      taatGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTatat              <  1:2067763/150‑1 (MQ=255)
                                                                                                                                          gCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATATTGAAAAGCAAGTTATATTTTc          >  2:1382007/1‑150 (MQ=255)
                                                                                                                                          gCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCTCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTc          >  2:395702/1‑150 (MQ=255)
                                                                                                                                           ccggtaaaGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc         <  2:1065448/144‑1 (MQ=255)
                                                                                                                                           cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACTCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc         <  1:2103635/150‑1 (MQ=255)
                                                                                                                                           cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc         >  1:1969134/1‑150 (MQ=255)
                                                                                                                                           cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTca         >  1:1977034/1‑149 (MQ=255)
                                                                                                                                             cGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCTCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAg       >  2:1060224/1‑150 (MQ=255)
                                                                                                                                             cGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAg       <  1:1375884/150‑1 (MQ=255)
                                                                                                                                             cGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAg       <  2:1359258/150‑1 (MQ=255)
                                                                                                                                                cAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGa    >  2:936631/1‑150 (MQ=255)
                                                                                                                                                  aGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATg  <  1:1293774/150‑1 (MQ=255)
                                                                                                                                                  aGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATg  <  2:1415138/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCTAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATG  >  NC_000913/394325‑394620

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: