Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,597,272	T→A	16.3%	pseudogene (716/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,272	0	T	A	16.3%	253.8 / 40.4	100	pseudogene (716/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base T (37/46);  new base A (8/8);  total (46/54)
Fisher's exact test for biased strand distribution p-value = 7.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

TTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTG  >  NC_000913/1597176‑1597365                                                                                                 |                                                                                               ttGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGc                                                                                             >  2:432767/1‑100 (MQ=255)     cTGTGTCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAcc                                                                                         <  1:618648/100‑1 (MQ=255)     cTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACCGAGTCTTGATTCTCCACATTCAGGGTTACGATACAGGTGGCACCGCCAc                                                                                          >  2:2046096/1‑99 (MQ=255)     cTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAcc                                                                                         <  1:1678757/100‑1 (MQ=255)     cTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCAAcc                                                                                         >  2:2046097/1‑100 (MQ=255)        tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGAt                                                                                      <  2:1125553/100‑1 (MQ=255)        tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                      <  2:1865507/100‑3 (MQ=255)        tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                      <  1:1345627/100‑3 (MQ=255)          ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATaa                                                                                    <  2:1409088/100‑1 (MQ=255)           cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAt                                                                                   <  1:1898077/100‑1 (MQ=255)           cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatc                                                                                   <  2:1774258/100‑3 (MQ=255)            aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATATCGGTGGCAATGCCACCGGTatcg                                                                                  <  2:2181999/100‑4 (MQ=39)            aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggtatcg                                                                                  <  2:569479/100‑7 (MQ=255)            aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatc                                                                                   <  1:284298/99‑3 (MQ=255)            aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCAGATACCGGTGGCAATGCCACCggtatcg                                                                                  <  2:2184088/100‑7 (MQ=255)             ttACCGAAATATCCATTATAGAGTCTGCGTATCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                 <  2:313397/100‑1 (MQ=255)             ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                 <  1:941044/100‑1 (MQ=255)             ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                 <  1:873433/100‑1 (MQ=255)               aCCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTa                                                                               >  2:1733324/1‑100 (MQ=255)               aCCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTa                                                                               >  1:294570/1‑100 (MQ=255)                ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAt                                                                              >  1:1242705/1‑100 (MQ=255)                ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAt                                                                              >  1:1221778/1‑100 (MQ=255)                 cGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAtt                                                                             >  1:1800776/1‑100 (MQ=255)                 cGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAtt                                                                             >  1:2020764/1‑100 (MQ=255)                  gAAATATCCATTATAGTGTCTGCGTACTTTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCTGTTGCAATGCCACCGGTatcg                                                                                  <  2:719265/94‑4 (MQ=255)                  gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                            <  2:2207909/100‑1 (MQ=255)                   aaaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAtttt                                                                           <  1:845056/100‑1 (MQ=255)                      tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAg                                                                        >  1:726140/1‑100 (MQ=255)                      tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAg                                                                        <  1:1896250/100‑1 (MQ=255)                      tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAg                                                                        >  2:1805613/1‑100 (MQ=255)                       atCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                       >  2:123465/1‑100 (MQ=255)                         ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                     <  2:669119/100‑1 (MQ=255)                          cATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                    >  1:1683973/1‑100 (MQ=255)                          cATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                    >  2:2110756/1‑100 (MQ=255)                             tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaac                                                                 <  2:1606108/100‑1 (MQ=255)                              ataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaac                                                                <  1:1519636/100‑3 (MQ=255)                               taGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca                                                               <  1:2035107/100‑1 (MQ=255)                                  agTCTGCTTACCGTGATGACAGAGTTTTCATCCTCTACATTCCGGGTTTCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTtaccaa‑cc                                                           <  2:1175937/100‑1 (MQ=255)                                  agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa‑cc                                                           <  1:1150674/100‑1 (MQ=255)                                  agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa‑cc                                                           <  2:930369/100‑1 (MQ=255)                                  agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa‑cc                                                           <  2:528692/100‑1 (MQ=255)                                     cTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑ccccc                                                        >  1:688388/1‑100 (MQ=255)                                       gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑ccccccc                                                      <  1:233348/100‑1 (MQ=255)                                        cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑cccccc                                                       >  1:46094/1‑98 (MQ=255)                                          tACCGTGATGACAGATTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTc                                                    <  2:840330/99‑1 (MQ=255)                                          tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTcg                                                   <  1:101296/100‑1 (MQ=255)                                          tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTcg                                                   >  1:1103750/1‑100 (MQ=255)                                            cgatcttgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTcgcg                                                 <  1:475250/94‑1 (MQ=255)                                            ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTcgcg                                                 >  1:598068/1‑100 (MQ=255)                                            ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTcgcg                                                 <  1:1860372/100‑1 (MQ=255)                                            ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATAGTATTCTCAGCCGTAACCAA‑ACCACCCTcg                                                   >  2:175877/1‑98 (MQ=255)                                              gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGaa                                               <  1:2148319/100‑1 (MQ=255)                                               tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAAt                                              <  1:996038/100‑1 (MQ=255)                                                 atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAg                                            <  1:2048049/100‑1 (MQ=255)                                                 atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAg                                            >  2:254532/1‑100 (MQ=255)                                                  tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGagatcg                                           >  2:475251/1‑95 (MQ=255)                                                    aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCg                                         <  1:106347/100‑1 (MQ=255)                                                     cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGTATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CTCCCCCTCGCGAATAGTCGc                                        <  2:252846/100‑1 (MQ=255)                                                     cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGc                                        >  1:1880493/1‑100 (MQ=255)                                                       gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCt                                      >  1:1690215/1‑100 (MQ=255)                                                         gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTcc                                    >  2:487302/1‑100 (MQ=255)                                                          tCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTccc                                   >  2:1878171/1‑100 (MQ=255)                                                           cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑ACCCCCCTCGCGAATAGTCGCCTcccc                                  >  1:1949591/1‑100 (MQ=255)                                                             tGATCCTGCACATTCAGGGTTAAGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAg                                >  2:1063826/1‑100 (MQ=255)                                                                 ccTGCAGATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGGGAATAGTCGCCTCCCCAGtt                              >  2:1666124/1‑98 (MQ=255)                                                                  cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                            <  1:161711/99‑1 (MQ=255)                                                                   tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCTCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                          <  1:1098341/100‑1 (MQ=255)                                                                     cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                        <  1:1812844/100‑1 (MQ=255)                                                                      acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                       <  1:1860650/100‑1 (MQ=255)                                                                       caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                      >  1:1610555/1‑100 (MQ=255)                                                                        aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                     <  2:1749552/100‑1 (MQ=255)                                                                        aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                     <  2:1242700/100‑1 (MQ=255)                                                                         ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                    <  1:444075/100‑1 (MQ=255)                                                                          tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                   <  1:2156709/100‑1 (MQ=255)                                                                            aGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                 >  1:695145/1‑100 (MQ=255)                                                                               atacgattatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt              >  2:180103/11‑100 (MQ=38)                                                                                ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg             <  2:1880481/100‑1 (MQ=255)                                                                                 tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt              >  1:1237000/1‑98 (MQ=255)                                                                                  cgattatcGGTGGCAATGCCACCGATAATCGTATTTTTATTCTTAATCAA‑CCCCCCCACCCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa           >  2:748947/8‑100 (MQ=255)                                                                                  ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa            <  1:1490280/100‑1 (MQ=255)                                                                                   gattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa          >  1:1219325/7‑100 (MQ=255)                                                                                   gattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa          >  2:825820/7‑100 (MQ=255)                                                                                   cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa          >  1:196870/1‑100 (MQ=255)                                                                                   cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa          >  2:289260/1‑100 (MQ=255)                                                                                   cGATAACGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTTCCCAGTTCCTTGATTACCAATTTGaaa          <  2:146962/100‑1 (MQ=255)                                                                                    gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt         <  1:1996061/100‑1 (MQ=255)                                                                                    gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt         >  1:2216478/1‑100 (MQ=255)                                                                                    gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt         <  2:1600187/100‑1 (MQ=255)                                                                                    gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt         <  2:1503361/100‑1 (MQ=255)                                                                                    attatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa          >  1:1940958/6‑99 (MQ=39)                                                                                     aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt        <  2:625/100‑1 (MQ=255)                                                                                     aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt        >  1:1168937/1‑100 (MQ=255)                                                                                      tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc       >  1:2000518/4‑100 (MQ=255)                                                                                      tACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc       <  1:307408/100‑1 (MQ=255)                                                                                       atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa      >  2:1216084/3‑100 (MQ=255)                                                                                       atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCCCGAATAGTCGCCTCCCCAGTTCCTTGAATACCAATTTGAAAtt        >  1:1796835/3‑98 (MQ=255)                                                                                         cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa     >  1:579607/1‑99 (MQ=255)                                                                                          ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTAAAt    >  1:649948/1‑99 (MQ=255)                                                                                            tGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg  <  1:866796/100‑1 (MQ=255)                                                                                             ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCCCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg  >  1:1833704/1‑98 (MQ=255)                                                                                                 |                                                                                               TTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTG  >  NC_000913/1597176‑1597365

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