Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,322,605 C→T pseudogene (210/372 nt) phnE ← defective phosphonate ABC transporter permease;transport; Transport of small molecules: Anions; membrane channel protein component of Pn transporter; membrane component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,322,6050CT97.0% 1348.9 / 17.8 403pseudogene (210/372 nt)phnEdefective phosphonate ABC transporter permease;transport; Transport of small molecules: Anions; membrane channel protein component of Pn transporter; membrane component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base T (181/208);  total (189/214)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 403 total aligned reads displayed.

TACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCC  >  NC_000913/4322457‑4322752
                                                                                                                                                    |                                                                                                                                                   
   gCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGa                                                                                                                                                <  8:813399/151‑1 (MQ=255)
    cAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCg                                                                                                                                                 >  4:301056/1‑149 (MQ=255)
    cAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCg                                                                                                                                                 >  5:747225/1‑149 (MQ=255)
    cAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAg                                                                                                                                               <  6:861295/151‑1 (MQ=255)
    cAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAAGCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAg                                                                                                                                               <  1:637305/151‑1 (MQ=255)
           tGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGt                                                                                                                                        >  2:680549/1‑151 (MQ=255)
             agagagGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTg                                                                                                                                      >  7:364063/1‑151 (MQ=255)
                  ggAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTc                                                                                                                                 <  3:872382/151‑1 (MQ=255)
                   gAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCg                                                                                                                                >  2:346341/1‑151 (MQ=255)
                     aaTCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAACCCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGaa                                                                                                                              <  3:68452/151‑1 (MQ=255)
                     aaTCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGa                                                                                                                               >  6:171816/1‑150 (MQ=255)
                     aaTCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGa                                                                                                                               >  1:246351/1‑150 (MQ=255)
                        ccAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCGTGCCGACGACGGTTGCAGAGCGGACGTTGGATTCGAAGCg                                                                                                                           >  8:1661690/1‑151 (MQ=255)
                        ccAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGc                                                                                                                            <  8:741773/150‑1 (MQ=255)
                           gcagGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATa                                                                                                                        <  5:998975/151‑1 (MQ=255)
                             agGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGAt                                                                                                                         >  6:1102076/1‑148 (MQ=255)
                             agGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATaga                                                                                                                      <  4:1198357/151‑1 (MQ=255)
                              ggCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATagag                                                                                                                     <  7:1134221/151‑1 (MQ=255)
                                  gagcGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGa                                                                                                                  <  6:698536/147‑1 (MQ=255)
                                  gACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAggg                                                                                                                   >  3:20099/1‑149 (MQ=255)
                                  gACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAggg                                                                                                                   >  4:706147/1‑149 (MQ=255)
                                   aCCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGt                                                                                                                >  2:1542304/1‑151 (MQ=255)
                                      gTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGCCCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAgg                                                                                                             >  2:1073565/1‑151 (MQ=255)
                                       tGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAgg                                                                                                             >  5:591527/1‑150 (MQ=255)
                                          cGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAgg                                                                                                             >  7:720607/1‑147 (MQ=255)
                                             tgatAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGCGAGTAGGAGATc                                                                                                        >  2:1373606/1‑149 (MQ=255)
                                              gatAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATcagc                                                                                                     >  7:1625156/1‑151 (MQ=255)
                                              gatAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATcagc                                                                                                     >  5:650045/1‑151 (MQ=255)
                                               atAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCCCAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATcagca                                                                                                    <  3:1016639/151‑1 (MQ=255)
                                                 aaGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGt                                                                                                  <  1:1000358/151‑1 (MQ=255)
                                                  aGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTg                                                                                                 >  8:1618226/1‑151 (MQ=255)
                                                      ccATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCAt                                                                                             >  1:895765/1‑151 (MQ=255)
                                                        aTCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATc                                                                                            >  8:600909/1‑150 (MQ=255)
                                                          cAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCAcc                                                                                         >  3:1470359/1‑151 (MQ=255)
                                                             ggCGCAGGTTTGATGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCtgt                                                                                      >  1:399243/1‑151 (MQ=255)
                                                              gcgcAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCtgtg                                                                                     <  6:619157/151‑1 (MQ=255)
                                                               cgcAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGGTCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTgg                                                                                    <  7:487706/151‑1 (MQ=255)
                                                                  aGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGgcag                                                                                 <  1:1504024/151‑1 (MQ=255)
                                                                  aGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGgcag                                                                                 <  4:558782/151‑1 (MQ=255)
                                                                   ggTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGgcagc                                                                                <  3:513152/151‑1 (MQ=255)
                                                                        gttgGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGcc                                                                           >  3:97844/1‑151 (MQ=255)
                                                                         ttgGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGCCGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCg                                                                          >  8:41881/1‑151 (MQ=255)
                                                                         ttgGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCg                                                                          >  4:1110525/1‑151 (MQ=255)
                                                                          tgGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCg                                                                          <  4:441178/150‑1 (MQ=255)
                                                                             aaCTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGt                                                                         >  8:535064/1‑148 (MQ=255)
                                                                              aCTGGAAACCGCGAAGCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTagag                                                                     <  7:1248065/151‑1 (MQ=255)
                                                                                 ggAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAgg                                                                    <  8:329888/149‑1 (MQ=255)
                                                                                 ggAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGAt                                                                  <  4:865706/151‑1 (MQ=255)
                                                                                     aCCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATctct                                                              >  7:1039406/1‑151 (MQ=255)
                                                                                     aCCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATctct                                                              >  3:300820/1‑151 (MQ=255)
                                                                                       cgcgAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTc                                                            <  1:1601152/151‑1 (MQ=255)
                                                                                       cgcgAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTc                                                            <  4:819931/151‑1 (MQ=255)
                                                                                        gcgAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCg                                                           <  6:1408554/151‑1 (MQ=255)
                                                                                         cgAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGa                                                          <  7:280527/151‑1 (MQ=255)
                                                                                          gAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCg                                                           >  6:942908/1‑149 (MQ=255)
                                                                                          gAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCg                                                           >  3:1313265/1‑149 (MQ=255)
                                                                                          gAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCg                                                           >  5:182717/1‑149 (MQ=255)
                                                                                              cGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCttg                                                     <  5:1046693/151‑1 (MQ=255)
                                                                                              cGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCtt                                                      <  2:782164/150‑1 (MQ=255)
                                                                                               gCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCttgt                                                    >  2:123588/1‑151 (MQ=255)
                                                                                                cTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCttgtt                                                   >  1:1331040/1‑151 (MQ=255)
                                                                                                   cccACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGc                                                <  7:629997/151‑1 (MQ=255)
                                                                                                   cccACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGc                                                >  7:1384531/1‑151 (MQ=255)
                                                                                                        aGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCgg                                           <  4:1385572/151‑1 (MQ=255)
                                                                                                         gggTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGt                                          <  1:1668420/151‑1 (MQ=255)
                                                                                                          ggTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTg                                         <  6:1211784/151‑1 (MQ=255)
                                                                                                            ggACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGgc                                       <  2:1596558/150‑1 (MQ=255)
                                                                                                             gACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGgc                                       <  5:467434/150‑1 (MQ=255)
                                                                                                              aCGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGGTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGgcgc                                     <  3:1076323/151‑1 (MQ=255)
                                                                                                              aCGCCGATCCCGCCTGCGCAGACAATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATAAGCAGTGGCATCACCTGTGGCAGCACGACGTAGAGTATCTCTTCGAGCTTGTTGGCACACGTGgcgc                                     >  6:15724/1‑151 (MQ=255)
                                                                                                               cGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGGATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGcgcg                                    <  2:232404/151‑1 (MQ=255)
                                                                                                               cGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGcgcg                                    <  3:592467/151‑1 (MQ=255)
                                                                                                                 ccGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGaa                                  <  6:471378/151‑1 (MQ=255)
                                                                                                                  cGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAAt                                 <  8:304817/151‑1 (MQ=255)
                                                                                                                   gATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATg                                >  5:1518755/1‑151 (MQ=255)
                                                                                                                     tCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTCGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGcc                              >  8:646744/1‑151 (MQ=255)
                                                                                                                        cGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTc                           >  7:806432/1‑151 (MQ=255)
                                                                                                                         gCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTcc                          >  8:898016/1‑151 (MQ=255)
                                                                                                                         gCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTcc                          <  2:44237/151‑1 (MQ=255)
                                                                                                                         gCCTGCGCCGACCATGCCGACGACGGTTGACGAGAGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGAGCGAATGCGTTcc                          >  7:1329850/1‑151 (MQ=255)
                                                                                                                           cTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCAc                        >  4:103551/1‑151 (MQ=255)
                                                                                                                            tGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCAcc                       >  5:1018844/1‑151 (MQ=255)
                                                                                                                             gcgcCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGGCTTCCAcgg                      <  2:608149/151‑3 (MQ=255)
                                                                                                                              cgcCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGAGCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCAccgg                     <  1:1169624/151‑1 (MQ=255)
                                                                                                                               gcCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCAccggg                    >  5:1668127/1‑151 (MQ=255)
                                                                                                                                ccGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCAccgggc                   >  8:1211749/1‑151 (MQ=255)
                                                                                                                                 cGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCAccgggcc                  <  3:135534/151‑1 (MQ=255)
                                                                                                                                 cGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCAccgggcc                  <  4:163147/151‑1 (MQ=255)
                                                                                                                                  gACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCAccgggccg                 <  5:1223024/151‑1 (MQ=255)
                                                                                                                                   aCCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCAccgggccgg                <  5:849173/151‑1 (MQ=255)
                                                                                                                                     cATGACGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGACATCAGCAGTGGCATCACCTGTGGCAGCACGCCGAAGAGGATCTCTTCGAGCTTCTTGGCAACGGTGCCGCGAAAGCCTTCCACccggccgggc              >  5:915151/1‑151 (MQ=255)
                                                                                                                                        gCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTc            >  8:112662/1‑150 (MQ=255)
                                                                                                                                         ccGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCaa          <  8:1536823/151‑1 (MQ=255)
                                                                                                                                         ccGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCa           >  3:359422/1‑150 (MQ=255)
                                                                                                                                          cgacgaCGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAAt         >  3:6148/1‑151 (MQ=255)
                                                                                                                                           gacgacGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATc        >  3:1252629/1‑151 (MQ=255)
                                                                                                                                            acgacgGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCg       <  8:503714/151‑1 (MQ=255)
                                                                                                                                              gacgGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCt     <  7:928652/151‑1 (MQ=255)
                                                                                                                                              gacgGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCt     >  7:1541339/1‑151 (MQ=255)
                                                                                                                                              gacgGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCt     <  6:388757/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCC  >  NC_000913/4322457‑4322752

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