Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	714,588	C→T	A344V (GCC→GTC)	pgm →	phosphoglucomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	714,588	0	C	T	97.4%	252.1 / ‑2.2	80	A344V (GCC→GTC)	pgm	phosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (44/32);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

TCGTCACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGCCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGT  >  NC_000913/714492‑714683                                                                                                 |                                                                                                tCGTCACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGt                                                                                              <  1:1272968/100‑1 (MQ=255) tCGTCACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGt                                                                                              >  2:577160/1‑100 (MQ=255)  cGTCACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTgtcgtc                                                                                             >  2:1298119/1‑100 (MQ=255)  cGTCACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTgtcgtc                                                                                             >  1:2083248/1‑100 (MQ=255)   gTCACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCATCATCGTCCGCAGTGGGGCAAAGATGTTGtcgtcg                                                                                            >  2:1235574/1‑100 (MQ=255)   gTCACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGtcgtcg                                                                                            >  2:504977/1‑100 (MQ=255)    tCACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCgg                                                                                           >  1:470942/1‑100 (MQ=255)     cACTCCGGGAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGt                                                                                          >  2:2484613/1‑100 (MQ=255)     cACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGt                                                                                          <  1:927552/100‑1 (MQ=255)     cACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGt                                                                                          >  2:2043626/1‑100 (MQ=255)       cTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTaa                                                                                        >  1:855695/1‑100 (MQ=255)          cGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTCCCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAAc                                                                                     <  1:2517883/100‑1 (MQ=255)           ggCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACg                                                                                    >  2:2568245/1‑100 (MQ=255)           ggCAGG‑TTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCAGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGc                                                                                   >  2:323084/1‑100 (MQ=255)            gCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGc                                                                                   >  1:1648135/1‑100 (MQ=255)               ggTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCggg                                                                                >  2:1921371/1‑100 (MQ=255)                 tttGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGtt                                                                              <  1:1571399/100‑1 (MQ=255)                 tttGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGtt                                                                              <  2:728127/100‑1 (MQ=255)                    gatgaATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCa                                                                           >  2:1677649/1‑100 (MQ=255)                    gatgaATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCa                                                                           >  2:729648/1‑100 (MQ=255)                      tgaATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATc                                                                         <  1:1804584/100‑1 (MQ=255)                      tgaATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATc                                                                         >  1:1254754/1‑100 (MQ=255)                      tgaATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATc                                                                         <  1:530034/100‑1 (MQ=255)                       gaATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCt                                                                        <  1:1989028/100‑1 (MQ=255)                       gaATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCt                                                                        >  2:1407984/1‑100 (MQ=255)                         aTCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGCTGTTGTCGTCGGTCCGCCGCTGGTCTCaccggc                                                                      >  2:1754696/1‑95 (MQ=255)                         aTCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGc                                                                      <  1:719558/100‑1 (MQ=255)                          tCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCg                                                                     <  2:1498020/100‑1 (MQ=255)                           ccGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCga                                                                    <  1:508934/100‑1 (MQ=255)                            cGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCgat                                                                   >  1:1154904/1‑100 (MQ=255)                            cGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCgat                                                                   >  2:26683/1‑100 (MQ=255)                               aCCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCgatgat                                                                >  1:1790089/1‑100 (MQ=255)                                  aCTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGGGATCATCGa                                                             >  2:2222089/1‑100 (MQ=255)                                     aCCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCg                                                          <  1:1887050/100‑1 (MQ=255)                                       cTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCgtg                                                        <  2:2569620/100‑1 (MQ=255)                                          gCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTgg                                                     <  1:1918249/100‑1 (MQ=255)                                          gCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGCCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTgg                                                     <  2:2018317/100‑1 (MQ=255)                                           cGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGt                                                    >  1:321230/1‑100 (MQ=255)                                               ggCAGTCAATTACCTGTTCCAGCATCGTGCGCAGTGGGGCAAATATGTTGTCGTCGGTAAAACGCTTGTTTCATCTGCGATGATCGACCGTGTGGTCAAc                                                <  2:1055299/100‑1 (MQ=255)                                               ggCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAAc                                                <  2:286443/100‑1 (MQ=255)                                                 caatcaatTACCTGTTCCAGCATCGGCCGCATTGTGGCCAAGGTGTTGTCGTCGGTAAAACACTGGTTTCATCTGCGATGATCGACCGGGTGGTCAACGa                                              >  2:273045/1‑100 (MQ=255)                                                 caatcaatTACCTGTTCCAGCATAGTCCGCAGTGGGGCAAAGATGTTGACGTCGGTAAAAAGCTGGTTTCATCTACGATGATCGACCGTGTGGCCAACGa                                              >  1:1736320/1‑100 (MQ=255)                                                      aatTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTgg                                         >  2:1107826/1‑100 (MQ=255)                                                      aatTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTgg                                         >  1:811602/1‑100 (MQ=255)                                                      aatTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGAGATGATCGACCGTGTGGTCAACGACTTgg                                         <  2:414205/100‑1 (MQ=255)                                                       atTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTggg                                        <  2:2112123/100‑1 (MQ=255)                                                       atTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTggg                                        >  2:2111917/1‑100 (MQ=255)                                                        ttACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGCCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACCTGGGc                                       <  1:1063762/100‑1 (MQ=255)                                                         tACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGcc                                      >  2:1581971/1‑100 (MQ=255)                                                          aCCTGTTCCAGCATCGTCCTCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCg                                     >  2:2577180/1‑100 (MQ=255)                                                           ccTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGt                                    >  1:1299424/1‑100 (MQ=255)                                                            cTGTTCCAGCGTCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTa                                   <  2:2215303/100‑1 (MQ=255)                                                             tGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTaa                                  >  1:2560454/1‑100 (MQ=255)                                                             tGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTaa                                  <  2:1400690/100‑1 (MQ=255)                                                             tGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGAGGATCGACCGTGTGGTCAAAGACTTGGGCCGTaa                                  >  1:2612456/1‑100 (MQ=255)                                                             tGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTAGCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTaa                                  <  2:1401023/100‑1 (MQ=255)                                                              gTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTaaa                                 <  1:833964/100‑1 (MQ=255)                                                              gTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTaaa                                 >  1:2027469/1‑100 (MQ=255)                                                                tCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACt                               >  1:2455127/1‑100 (MQ=255)                                                                 ccAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTg                              >  2:520688/1‑100 (MQ=255)                                                                  cAGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTgg                             >  2:209441/1‑100 (MQ=255)                                                                   aGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGt                            <  2:840955/100‑1 (MQ=255)                                                                   aGCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGt                            <  2:17162/100‑1 (MQ=255)                                                                    gCATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTa                           <  1:2339354/100‑1 (MQ=255)                                                                     cATCGTCCGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAg                          >  2:1852605/1‑100 (MQ=255)                                                                       tCGTCCGCAGTGGGGCAAAGATGTTGTCATCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGaa                        >  2:1913409/1‑100 (MQ=255)                                                                            cGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTccc                   <  1:1434766/100‑1 (MQ=255)                                                                            cGCAGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTccc                   >  1:1806422/1‑100 (MQ=255)                                                                               aGTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGt                <  1:2531371/100‑1 (MQ=255)                                                                                gTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTa               >  1:1883386/1‑100 (MQ=255)                                                                                gTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTa               <  1:1690692/100‑1 (MQ=255)                                                                                gTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTa               <  2:52683/100‑1 (MQ=255)                                                                                gTGGGGCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTa               >  1:1882574/1‑100 (MQ=255)                                                                                  ggggCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAgg             <  2:1556630/100‑1 (MQ=255)                                                                                  ggggCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAgg             >  2:2355886/1‑100 (MQ=255)                                                                                   gggCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGt            >  2:2391815/1‑100 (MQ=255)                                                                                    ggCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGtt           >  2:185667/1‑100 (MQ=255)                                                                                    ggCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGtt           >  2:337922/1‑100 (MQ=255)                                                                                    ggCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGtt           <  1:213059/100‑1 (MQ=255)                                                                                    ggCAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGtt           >  2:781887/1‑100 (MQ=255)                                                                                      cAAAGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTc         >  1:2202476/1‑100 (MQ=255)                                                                                       aaaGATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCa        <  1:2255378/100‑1 (MQ=255)                                                                                          gATGTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAAt     >  1:1498277/1‑100 (MQ=255)                                                                                             gTTGTCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGt  <  1:41472/100‑1 (MQ=255)                                                                                                 |                                                                                                TCGTCACTCCGGCAGGTTTGATGAATCCGAACCACTACCTGGCGGTGGCAATCAATTACCTGTTCCAGCATCGTCCGCAGTGGGGCAAAGATGTTGCCGTCGGTAAAACGCTGGTTTCATCTGCGATGATCGACCGTGTGGTCAACGACTTGGGCCGTAAACTGGTAGAAGTCCCGGTAGGTTTCAAATGGT  >  NC_000913/714492‑714683

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: ‑