Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,393 T→A 5.9% intergenic (+58/+23) pheP → / ← ybdG phenylalanine transporter/mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913603,3930TA 5.9% 1149.5 / 46.1 340intergenic (+58/+23)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base T (166/152);  new base A (9/11);  total (176/164)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

Only 100 of 352 total aligned reads displayed.

CACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCATA  >  NC_000913/603245‑603541
                                                                                                                                                    |                                                                                                                                                    
 aagatggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAACACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAtt                                                                                                                                                     <  2:260341/147‑1 (MQ=255)
    atggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa                                                                                                                                                  <  2:785220/149‑1 (MQ=255)
     tggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAGCCGCTCTCATCCATTCGAt                                                                                                                                                 <  2:2114764/149‑1 (MQ=255)
      ggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                                                                                >  2:1479398/1‑149 (MQ=255)
       gatgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                                                                               <  3:1971792/149‑1 (MQ=255)
       gatgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                                                                               <  2:1331824/149‑1 (MQ=255)
       gatgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                                                                               >  3:765352/1‑149 (MQ=255)
         tgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgaga                                                                                                                                             <  4:476836/149‑1 (MQ=255)
           atatGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGc                                                                                                                                           <  4:1717923/149‑1 (MQ=255)
             atGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                                                                         <  2:888628/149‑1 (MQ=255)
                   ttGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttttt                                                                                                                                   >  2:1785165/1‑149 (MQ=255)
                      tCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTaa                                                                                                                                <  1:1785165/149‑1 (MQ=255)
                         gCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttt                                                                                                                                      <  3:2211221/140‑1 (MQ=255)
                         gCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTa                                                                                                                             >  4:426530/1‑149 (MQ=255)
                         gCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTa                                                                                                                             <  4:763322/149‑1 (MQ=255)
                          cGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTAc                                                                                                                            <  2:1057993/149‑1 (MQ=255)
                           gATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACt                                                                                                                           >  4:264617/1‑149 (MQ=255)
                            aTCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTg                                                                                                                          >  4:1054720/1‑149 (MQ=255)
                             tCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                                                                 >  3:202257/1‑125 (MQ=255)
                             tCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTAc                                                                                                                            >  4:601915/1‑146 (MQ=255)
                                  ctgcCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                                                                    <  4:1378186/149‑1 (MQ=255)
                                   tgcCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                                                                   <  4:783063/149‑1 (MQ=255)
                                   tgcCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                                                                   <  1:1982573/149‑1 (MQ=255)
                                    gcCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                                                                  >  4:1821211/1‑149 (MQ=255)
                                          gtgGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                                                                            >  2:62499/1‑149 (MQ=255)
                                             gATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgaga                                                                                                                                             >  3:726892/1‑113 (MQ=255)
                                             gATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgaga                                                                                                                                             <  4:726892/113‑1 (MQ=255)
                                             gATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCa                                                                                                         >  4:906849/1‑149 (MQ=255)
                                               ttGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCaga                                                                                                       >  3:665023/1‑149 (MQ=255)
                                                 gTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagaga                                                                                                     <  1:1599975/149‑1 (MQ=255)
                                                  tATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                                                                 <  3:1393358/104‑1 (MQ=255)
                                                    ttCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCg                                                                                                  >  4:551071/1‑149 (MQ=255)
                                                     tCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                                                                   <  4:2312776/131‑1 (MQ=255)
                                                      ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGaa                                                                                                <  3:1486245/149‑1 (MQ=255)
                                                         gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAatat                                                                                             <  1:61208/149‑1 (MQ=255)
                                                            tATGGCATTTAAAACGCTGCGTCAGAAATAAGCCATTCACGCTACATCTGACAAAACGATGTCAACTATCCGAAACCGCGCACATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAt                                                                                          <  2:371995/149‑1 (MQ=255)
                                                               ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                                                                       <  4:2034267/149‑1 (MQ=255)
                                                               ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                                                                       <  2:985041/149‑1 (MQ=255)
                                                                 cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAATGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCg                                                                                     <  2:305829/149‑1 (MQ=255)
                                                                 cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCg                                                                                     >  1:2165113/1‑149 (MQ=255)
                                                                  aTTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAtt                                                                                                                                                     >  1:2290150/1‑84 (MQ=255)
                                                                       aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGc                                                                               <  1:1479398/149‑1 (MQ=255)
                                                                                 tCGGAAATAAGGCATTCACGCTCCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                                                                           >  1:1529163/1‑79 (MQ=37)
                                                                                 tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                                                                        <  1:1252270/146‑1 (MQ=255)
                                                                                 tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                                                                     <  3:1894625/149‑1 (MQ=255)
                                                                                 tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                                                                     >  4:2448135/1‑149 (MQ=255)
                                                                                 tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                                                                     <  1:1932100/149‑1 (MQ=255)
                                                                                 tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                                                                     <  1:1731585/149‑1 (MQ=255)
                                                                                    gAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTc                                                                  <  3:373975/149‑1 (MQ=255)
                                                                                     aaaTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCg                                                                 >  1:1006976/1‑149 (MQ=255)
                                                                                         aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa                                                             <  2:1006976/149‑1 (MQ=255)
                                                                                           ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                                                                  >  3:559693/1‑126 (MQ=255)
                                                                                              aTTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTAGTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaaa                                                        >  4:1469992/1‑149 (MQ=255)
                                                                                              aTTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaaa                                                        >  4:2212244/1‑149 (MQ=255)
                                                                                               ttCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc                                                       >  3:1306975/1‑149 (MQ=255)
                                                                                                  aCGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                                                                     >  1:1414419/1‑132 (MQ=255)
                                                                                                  aCGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActct                                                    >  3:5397/1‑149 (MQ=255)
                                                                                                  aCGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActct                                                    >  2:119039/1‑149 (MQ=255)
                                                                                                    gCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTc                                                  >  3:225614/1‑149 (MQ=255)
                                                                                                     cTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg                                                 <  2:232453/149‑1 (MQ=255)
                                                                                                     cTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg                                                 <  4:199822/149‑1 (MQ=255)
                                                                                                        cATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACa                                              >  2:1869815/1‑149 (MQ=255)
                                                                                                         aTCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACaa                                             <  3:912705/149‑1 (MQ=255)
                                                                                                          tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                                                                  >  3:2420382/1‑79 (MQ=255)
                                                                                                            cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATc                                           >  2:743382/1‑148 (MQ=255)
                                                                                                            cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCg                                          >  2:1525368/1‑149 (MQ=255)
                                                                                                             gACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATACACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGc                                         <  1:590430/149‑1 (MQ=255)
                                                                                                                aaaaCGATGTCAACCATCCGAAACCGCTCTCATCTATTCGATGAGAGCGGTTTTTTAATTACTGCTTAAATGCACCCGCCCGAGAGCGAATATCATTGCCGGTTGGGGGCTGGTGAAGTCGCAGAACAAACTCTTCGACCATAGCaaaa                                     >  3:2333517/1‑148 (MQ=255)
                                                                                                                aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCaaa                                      <  2:1267019/149‑1 (MQ=255)
                                                                                                                aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCaaa                                      <  3:1640451/149‑1 (MQ=255)
                                                                                                                  aaCGATTTCAACCATCCGAAAGCGCTCTCATCCATTCGTTTAGTGCGGTTTTTTTAATTACTGCTTAGATGCACCCGCCAGAGAGCGAATGCTATTGACGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTACGACAATCGCAAata                                    >  2:2288647/1‑149 (MQ=255)
                                                                                                                     gATGTCAAACATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATAtgt                                 <  2:2056128/149‑1 (MQ=255)
                                                                                                                       tGTCAACCATCCGAAACCGCTCTCATCCATGCGAGGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGa                               <  1:442460/149‑1 (MQ=255)
                                                                                                                        gTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGAt                              <  1:472648/149‑1 (MQ=255)
                                                                                                                         tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATc                             >  4:2261710/1‑149 (MQ=255)
                                                                                                                           aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGa                           >  1:1128782/1‑149 (MQ=255)
                                                                                                                           aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGa                           >  2:443653/1‑149 (MQ=255)
                                                                                                                               taaaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCaaa                                      >  2:1614983/6‑134 (MQ=38)
                                                                                                                               aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAatat                       <  1:268933/149‑1 (MQ=255)
                                                                                                                                   gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAat                         >  1:907771/1‑143 (MQ=255)
                                                                                                                                   gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCa                   >  3:1222563/1‑149 (MQ=255)
                                                                                                                                        cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCaaa                                      >  3:1694721/1‑125 (MQ=255)
                                                                                                                                         gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGa             >  2:1839091/1‑149 (MQ=255)
                                                                                                                                          ctctCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGaa            >  1:828784/1‑149 (MQ=255)
                                                                                                                                           tctcATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGa                                                <  2:2386317/112‑1 (MQ=255)
                                                                                                                                           tctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATc                             >  1:367496/1‑131 (MQ=255)
                                                                                                                                             tcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGc         >  1:714724/1‑149 (MQ=255)
                                                                                                                                              cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCt        <  3:437858/149‑1 (MQ=255)
                                                                                                                                               aTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaaa                                                        >  4:1306066/1‑100 (MQ=255)
                                                                                                                                               aTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGc                 <  1:2474180/139‑1 (MQ=255)
                                                                                                                                               aTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCtt       >  4:1624253/1‑149 (MQ=255)
                                                                                                                                               aTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCtt       <  3:845529/149‑1 (MQ=255)
                                                                                                                                               aTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCtt       <  4:687323/149‑1 (MQ=255)
                                                                                                                                                tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCttt      >  3:1542675/1‑149 (MQ=255)
                                                                                                                                                 ccATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTc     <  1:2366318/149‑1 (MQ=255)
                                                                                                                                                 ccATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATACACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTc     >  4:2141472/1‑149 (MQ=255)
                                                                                                                                                   aTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCat   >  3:1827414/1‑149 (MQ=255)
                                                                                                                                                    ttCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTc                                                                  >  4:1132392/1‑85 (MQ=255)
                                                                                                                                                    ttCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCata  <  3:1054720/149‑1 (MQ=255)
                                                                                                                                                    ttCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCata  >  1:1835524/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCATA  >  NC_000913/603245‑603541

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