Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,226,115	A→C	8.4%	intergenic (‑35/‑485)	minC ← / → ycgJ	Z‑ring positioning protein MinC/PF05666 family protein YcgJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,226,115	0	A	C	8.4%	237.8 / 12.2	83	intergenic (‑35/‑485)	minC/ycgJ	Z‑ring positioning protein MinC/PF05666 family protein YcgJ
Reads supporting (aligned to +/- strand):  ref base A (36/40);  new base C (5/2);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

GAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAAC  >  NC_000913/1225969‑1226262                                                                                                                                                   |                                                                                                                                                    gAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAgcg                                                                                                                                                  >  1:1564668/1‑150 (MQ=255) gAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAgcg                                                                                                                                                  <  1:1477694/150‑1 (MQ=255)       ttttGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCgat                                                                                                                                            <  1:1153985/150‑1 (MQ=255)         ttGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCgatga                                                                                                                                          >  2:476288/1‑150 (MQ=255)           gTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgatt                                                                                                                                        <  2:454629/150‑2 (MQ=255)             cTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGtt                                                                                                                                      <  2:169077/150‑1 (MQ=255)              ttCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTc                                                                                                                                     <  1:33809/150‑1 (MQ=255)               tCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGCCATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTcc                                                                                                                                    >  2:1613381/1‑150 (MQ=255)               tCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTcc                                                                                                                                    <  1:493744/150‑1 (MQ=255)               tCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTcc                                                                                                                                    <  2:1564668/150‑1 (MQ=255)                ccAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCg                                                                                                                                   >  2:73585/1‑150 (MQ=255)                ccAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCg                                                                                                                                   <  2:1024736/150‑1 (MQ=255)                ccAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGAATTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCg                                                                                                                                   >  2:1126364/1‑150 (MQ=255)                   gcgcCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAg                                                                                                                                <  2:1335424/150‑1 (MQ=255)                   gcgcCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAg                                                                                                                                <  2:1521625/150‑1 (MQ=255)                     gcCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGAc                                                                                                                              >  1:1161523/1‑150 (MQ=255)                     gcCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGAc                                                                                                                              >  2:1498264/1‑150 (MQ=255)                      ccTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACt                                                                                                                             >  1:1073546/1‑150 (MQ=255)                            ggATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAgg                                                                                                                       >  1:1412953/1‑150 (MQ=255)                                aaCCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATg                                                                                                                   <  1:282614/150‑1 (MQ=255)                                 aCCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGt                                                                                                                  <  1:1389530/150‑1 (MQ=255)                                 aCCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGt                                                                                                                  >  1:215592/1‑150 (MQ=255)                                 aCCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGt                                                                                                                  <  2:1446412/150‑1 (MQ=255)                                       ggTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGtt                                                                                                                 >  1:413456/1‑145 (MQ=255)                                       ggTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGtt                                                                                                                 <  2:413456/145‑1 (MQ=255)                                        gTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                                                           <  2:1471114/150‑1 (MQ=255)                                            tGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAACATACCTTATAGTCttcg                                                                                                       <  2:635152/150‑4 (MQ=255)                                             gCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAt                                                                                                      <  1:1232969/150‑1 (MQ=255)                                               cTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGATATAGTCTGGATta                                                                                                    <  1:1441233/150‑1 (MQ=255)                                                tCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATtat                                                                                                   >  1:956480/1‑150 (MQ=255)                                                  aTGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTatat                                                                                                 <  1:796625/150‑1 (MQ=255)                                                   tGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATAtt                                                                                                >  2:960870/1‑150 (MQ=255)                                                    gCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTg                                                                                               >  1:329070/1‑150 (MQ=255)                                                     cAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGa                                                                                              <  1:1116048/150‑1 (MQ=255)                                                     cAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGa                                                                                              >  1:1148341/1‑150 (MQ=255)                                                     cAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGa                                                                                              >  2:178340/1‑150 (MQ=255)                                                     cAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGa                                                                                              <  1:465055/150‑1 (MQ=255)                                                           aaCCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAg                                                                                        <  2:1039872/150‑1 (MQ=255)                                                               acaGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAc                                                                                    >  1:764195/1‑150 (MQ=255)                                                                caGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                                                   >  1:717372/1‑150 (MQ=255)                                                                  gATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCt                                                                                 <  1:1178517/150‑1 (MQ=255)                                                                   aTAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTc                                                                                <  2:627361/150‑1 (MQ=255)                                                                      aaGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCa                                                                             >  2:386365/1‑150 (MQ=255)                                                                        gTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCAtt                                                                           <  2:898032/150‑1 (MQ=255)                                                                           aaGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTAt                                                                        <  1:745243/150‑1 (MQ=255)                                                                               tACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCag                                                                    <  1:73585/150‑1 (MQ=255)                                                                                aCTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCaga                                                                   >  2:1591002/1‑150 (MQ=255)                                                                                       ttAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAg                                                            <  2:553490/150‑1 (MQ=255)                                                                                       ttAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAg                                                            <  2:764195/150‑1 (MQ=255)                                                                                         aaGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTc                                                          <  1:1085249/150‑1 (MQ=255)                                                                                           gCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCtt                                                        >  2:504917/1‑149 (MQ=255)                                                                                            cTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCAGGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTAt                                                       >  2:393899/1‑150 (MQ=255)                                                                                              cGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTc                                                     <  2:215592/150‑1 (MQ=255)                                                                                                    gCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGAt                                               <  2:30197/150‑1 (MQ=255)                                                                                                      gtgtTTGACATCCTGGCCTTACTCAATAAGCTATAAACCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATaa                                             <  2:1021491/150‑1 (MQ=255)                                                                                                        gtTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAAt                                           >  1:779368/1‑150 (MQ=255)                                                                                                          ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCa                                                                                          <  2:298501/101‑1 (MQ=255)                                                                                                          ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCa                                                                                          >  1:298501/1‑101 (MQ=255)                                                                                                             aCATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgcgc                                      >  1:180033/1‑150 (MQ=255)                                                                                                              cATCCTTGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCt                                     <  2:956480/150‑1 (MQ=255)                                                                                                                     gCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCCAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCAGTGATAAATGGCGCTGATTCAt                              >  2:504522/1‑150 (MQ=255)                                                                                                                      ccTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATa                             <  2:395504/150‑1 (MQ=255)                                                                                                                       cTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAg                            >  1:157343/1‑150 (MQ=255)                                                                                                                           cTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTa                        >  2:481460/1‑150 (MQ=255)                                                                                                                               aTTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTaaaaa                    >  1:72131/1‑150 (MQ=255)                                                                                                                                 tAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTAAAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGAGCTGATTCATAGCTTAAAAAAt                  >  2:859354/1‑150 (MQ=255)                                                                                                                                   gCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCt                          <  2:1181327/140‑1 (MQ=255)                                                                                                                                   gCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCt                          >  1:1181327/1‑140 (MQ=255)                                                                                                                                    cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAcc               <  1:476288/150‑1 (MQ=255)                                                                                                                                    cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAcc               >  2:1554162/1‑150 (MQ=255)                                                                                                                                      gaacatcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATACTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCt             >  1:1028541/5‑150 (MQ=255)                                                                                                                                        tAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTg           <  2:329070/150‑1 (MQ=255)                                                                                                                                         catcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGt          >  1:1009581/2‑150 (MQ=255)                                                                                                                                         catcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGt          >  2:211207/2‑150 (MQ=255)                                                                                                                                         catcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGt          >  2:1050547/2‑150 (MQ=255)                                                                                                                                          atcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACGCTTGta         >  2:625783/1‑149 (MQ=255)                                                                                                                                          atcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtc         <  1:1050547/150‑1 (MQ=255)                                                                                                                                          atcatcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtc         >  1:511680/1‑150 (MQ=255)                                                                                                                                            catcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtcaa       >  1:1388484/1‑150 (MQ=255)                                                                                                                                            catcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCATTGtcaa       >  1:300011/1‑150 (MQ=255)                                                                                                                                              tcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtcaatc     <  1:1498264/150‑1 (MQ=255)                                                                                                                                                gCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtcaatcaa   >  2:441039/1‑150 (MQ=255)                                                                                                                                                 ccAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAAc  <  1:1591002/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAAC  >  NC_000913/1225969‑1226262

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