Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,528,194	A→T	8.1%	pseudogene (293/2037 nt)	rhsE →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,528,194	0	A	T	8.1%	226.3 / 8.2	86	pseudogene (293/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base A (47/31);  new base T (5/2);  total (53/33)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

TACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACAACGGCGA  >  NC_000913/1528049‑1528338                                                                                                                                                  |                                                                                                                                                 ttcaccGACCAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCGGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTGCACGCGCGACCGCCTGCACCGTGAGACg                                                                                                                                              <  1:329456/145‑1 (MQ=255)    aaCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTg                                                                                                                                           >  1:149336/1‑150 (MQ=255)     aCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGAGCGACCGCCTGCACCGTGAGACGGTgc                                                                                                                                          >  1:240391/1‑150 (MQ=255)       gAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTCTCACGGTgcag                                                                                                                                        <  1:540877/150‑3 (MQ=255)        aaCAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCa                                                                                                                                       >  2:1655097/1‑150 (MQ=255)         aCAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTCTCACGGTgcaggc                                                                                                                                      <  1:1421202/150‑5 (MQ=255)            ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCtt                                                                                                                                   >  1:651712/1‑150 (MQ=255)            ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCtt                                                                                                                                   >  1:1371325/1‑150 (MQ=255)              ggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCg                                                                                                                                 >  2:1438729/1‑150 (MQ=255)                   gCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgcagc                                                                                                                            >  1:1389410/1‑150 (MQ=255)                     aaaCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCAt                                                                                                                          <  2:226682/150‑1 (MQ=255)                       aCCGCGTCACGCCGGACAGTCTGCCGCCAGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATgg                                                                                                                        >  1:1471828/1‑150 (MQ=255)                        ccGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGTATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggc                                                                                                                       >  1:408287/1‑150 (MQ=255)                        ccGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCACTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggc                                                                                                                       <  2:829528/150‑1 (MQ=255)                         cgcgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCCCCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                                                                                                      >  1:1679939/1‑150 (MQ=255)                           cggcACGCCGGACAGGCTGCGGCCGGTGGAGTGGCGGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcagg                                                                                                                    <  1:89833/147‑1 (MQ=255)                           cgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcagg                                                                                                                    >  1:1481645/1‑150 (MQ=255)                             tCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggca                                                                                                                  >  2:610405/1‑150 (MQ=255)                              cACGCCGGACAGTCTGCCGCCGGTGTAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGgcaggcag                                                                                                                 <  1:441473/150‑1 (MQ=255)                              cACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGgcaggcag                                                                                                                 >  1:346072/1‑150 (MQ=255)                               aCGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGt                                                                                                                <  1:922195/150‑1 (MQ=255)                               aCGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGt                                                                                                                >  1:1079229/1‑150 (MQ=255)                                cGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAATa                                                                                                               >  2:16921/1‑150 (MQ=255)                                     gACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGcc                                                                                                          <  2:1371325/150‑1 (MQ=255)                                       cAGTCTGCCGACGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATCGCAGGCAGTACTGCCGc                                                                                                        >  1:146170/1‑150 (MQ=255)                                         gTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCat                                                                                                      >  2:475836/1‑150 (MQ=255)                                            tGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCAtata                                                                                                   <  2:1481645/150‑1 (MQ=255)                                             gccgccGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaa                                                                                                  <  1:1655097/150‑1 (MQ=255)                                              ccgccgGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaaa                                                                                                 <  2:689673/150‑1 (MQ=255)                                                  cgGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGa                                                                                             <  2:1337245/150‑1 (MQ=255)                                                    gTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAcc                                                                                           <  2:583160/150‑1 (MQ=255)                                                     tGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCa                                                                                          >  2:1350848/1‑150 (MQ=255)                                                      ggAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAg                                                                                         >  1:1487179/1‑150 (MQ=255)                                                       gAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGc                                                                                        <  2:1528260/150‑1 (MQ=255)                                                       gAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGc                                                                                        <  1:632974/150‑1 (MQ=255)                                                        aGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGcc                                                                                       >  2:1076403/1‑149 (MQ=255)                                                         gTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCTGTAATGCCGCATATAAACTGACCAGcac                                                                                      <  2:651712/150‑1 (MQ=255)                                                           ggCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACAt                                                                                    >  2:803890/1‑150 (MQ=255)                                                             cTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATac                                                                                  <  1:569892/150‑1 (MQ=255)                                                                  cTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGGGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACAcccc                                                                             <  2:1084016/150‑1 (MQ=255)                                                                      ggCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGCGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACCCCCCCGcc                                                                         >  1:553535/1‑149 (MQ=255)                                                                      ggCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCa                                                                         >  1:421992/1‑150 (MQ=255)                                                                             gCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAAAGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGt                                                                  >  2:259109/1‑150 (MQ=255)                                                                              cTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGtt                                                                 <  2:149336/150‑1 (MQ=255)                                                                                  cTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACag                                                             >  1:148534/1‑150 (MQ=255)                                                                                    tGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACagag                                                           <  2:1487179/150‑1 (MQ=255)                                                                                      cGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGcc                                                         <  2:202771/150‑1 (MQ=255)                                                                                      cGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGcc                                                         >  1:712741/1‑150 (MQ=255)                                                                                          cATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGCCCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCCGCCCAAACACCCCCGCAGGCCAGTTACAGAGCCAGCa                                                     >  2:136591/1‑150 (MQ=255)                                                                                          cATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGcc                                                     >  2:1347463/1‑149 (MQ=255)                                                                                           ttGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCAc                                                    <  2:350938/149‑1 (MQ=255)                                                                                                          aCGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGgtg                                     <  1:475836/150‑1 (MQ=255)                                                                                                          aCGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGgtg                                     <  2:1471828/150‑1 (MQ=255)                                                                                                           cgccgcTGCAGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCAGCACCTGAACAGCCTGgtgg                                    >  2:266528/1‑149 (MQ=255)                                                                                                            gccgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGCCGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCCGCACCTGAACAGCCTGGTGTa                                   >  1:761924/1‑150 (MQ=255)                                                                                                            gccgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTa                                   <  1:1650099/150‑1 (MQ=255)                                                                                                            gccgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTa                                   >  1:881202/1‑150 (MQ=255)                                                                                                             ccgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGGGTAt                                  >  1:432529/1‑150 (MQ=255)                                                                                                              cgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCAGCACCTGACCAGCCTGGTGTATg                                 >  1:1533348/1‑150 (MQ=255)                                                                                                              cgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATg                                 >  1:455524/1‑150 (MQ=255)                                                                                                               gctgctGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCCGGCCAGTTACAGAGCCAGCACCTGAACAGCTTGGTGTATGa                                >  1:809596/1‑150 (MQ=255)                                                                                                                ctgctgGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGAc                               <  2:346072/150‑1 (MQ=255)                                                                                                                ctgctgGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGAc                               <  1:1350848/150‑1 (MQ=255)                                                                                                                 tgctgGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGAcc                              >  2:1401914/1‑150 (MQ=255)                                                                                                                 tgctgGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGAcc                              >  1:1275710/1‑150 (MQ=255)                                                                                                                   ctgGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCCGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTCTGCCCGt                            >  2:268616/1‑150 (MQ=255)                                                                                                                     ggAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACCCCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGa                          >  2:464341/1‑150 (MQ=255)                                                                                                                      gAGTTCACGCGCGACCGCCTGCACCGTTAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACAAACACCCCCGCAGGCCAGTTACAGAGCCAGCACCAGAACATCCTGGTGTATGACCGAGAc                         >  2:1034889/1‑150 (MQ=17)                                                                                                                          tCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACg                     >  1:1060920/1‑150 (MQ=255)                                                                                                                            aCGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCTGCAGCATGTCATGCAGTAATGCCTCATATATACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACggg                   <  2:263195/150‑1 (MQ=14)                                                                                                                             cgcgcgACCGCCTGCTCCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGt                  >  1:981901/1‑150 (MQ=255)                                                                                                                             cgcgcgACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACCTACACCCCCGCAGGCCAGTTAGAGAGCCCGCACCTGAACAGCCTGGTGTATGACCGTGACTACgggg                  >  1:1413385/1‑149 (MQ=255)                                                                                                                               cgcgACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTgg                <  2:633948/150‑1 (MQ=255)                                                                                                                                 cgACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGaa              >  2:442326/1‑150 (MQ=255)                                                                                                                                   aCCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATg            >  2:880963/1‑150 (MQ=255)                                                                                                                                   aCCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATg            <  2:421992/150‑1 (MQ=255)                                                                                                                                     cGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGAc          >  2:605760/1‑150 (MQ=255)                                                                                                                                      gCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:338111/1‑150 (MQ=40)                                                                                                                                      ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:577686/5‑150 (MQ=255)                                                                                                                                      ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:422716/5‑150 (MQ=255)                                                                                                                                      ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:584393/5‑150 (MQ=255)                                                                                                                                      ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:1247999/5‑150 (MQ=255)                                                                                                                                      ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:1506779/5‑150 (MQ=255)                                                                                                                                       ccTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACaa        <  1:1401914/150‑1 (MQ=40)                                                                                                                                            aCCGTGAGACGGTGCGCAGATTCGGCAGCATGGAAGTCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCAAGTTACAGAGCCAGAACCTGAACATCCTGGTGTATGACCGTTACTACGGGTGGAATGACAACGGCg   <  2:1046170/150‑1 (MQ=14)                                                                                                                                             ccGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACAACGGCGa  >  1:744543/1‑150 (MQ=39)                                                                                                                                                  |                                                                                                                                                 TACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACAACGGCGA  >  NC_000913/1528049‑1528338

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