Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	179,327	A→T	100%	I31F (ATC→TTC)	dgt →	deoxyguanosine triphosphate triphosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	179,327	0	A	T	100.0%	319.6 / NA	95	I31F (ATC→TTC)	dgt	deoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (53/42);  total (53/42)

GAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCG  >  NC_000913/179231‑179423                                                                                                 |                                                                                                 gagGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCt                                                                                               >  2:2203845/1‑100 (MQ=255)  agGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCtt                                                                                              >  1:63987/1‑100 (MQ=255)      cATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCgag                                                                                          >  1:1629601/1‑100 (MQ=255)        tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGagag                                                                                        <  2:939944/100‑1 (MQ=255)        tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGagag                                                                                        <  2:434843/100‑1 (MQ=255)        tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGagag                                                                                        >  1:1598326/1‑100 (MQ=255)              aGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATcg                                                                                  <  1:1028736/100‑1 (MQ=255)               gattgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATcgc                                                                                 >  2:1292218/1‑100 (MQ=255)                attgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATcgcg                                                                                <  1:2405405/100‑1 (MQ=255)                  tgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCggg                                                                              <  1:1556943/100‑1 (MQ=255)                   gttttCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACAACCGGGCGTTAAAACCGAACATGAGATCCTGCGGTCCTTCGTGAGCGATcgc                                                                                 >  2:203945/3‑96 (MQ=255)                    attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCg                                                                            >  2:1547252/1‑100 (MQ=255)                     tttCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGt                                                                           >  2:2289844/1‑100 (MQ=255)                      ttCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTa                                                                          >  1:637784/1‑100 (MQ=255)                      ttCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTa                                                                          <  1:1239488/100‑1 (MQ=255)                       tCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTat                                                                         >  1:2205221/1‑100 (MQ=255)                       tCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTat                                                                         >  1:998646/1‑100 (MQ=255)                        ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatc                                                                        >  1:915951/1‑100 (MQ=255)                         cGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatca                                                                       >  1:82067/1‑100 (MQ=255)                         cGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatca                                                                       >  2:268083/1‑100 (MQ=255)                          gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                      <  1:1207219/100‑1 (MQ=255)                          gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                      >  2:1711440/1‑100 (MQ=255)                          gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                      <  1:2203860/100‑1 (MQ=255)                          gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                      >  2:1355365/1‑100 (MQ=255)                          gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                      >  1:509776/1‑100 (MQ=255)                          gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                      >  2:98347/1‑100 (MQ=255)                           aaaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcatc                                                                     <  2:1389736/100‑1 (MQ=255)                           aaaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcatc                                                                     <  2:541613/100‑1 (MQ=255)                              aaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                  <  2:1307639/100‑1 (MQ=255)                              aaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                  <  1:2133949/100‑1 (MQ=255)                              aaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                  <  2:2385753/100‑1 (MQ=255)                              aaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                  <  1:2280787/100‑1 (MQ=255)                               aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAct                                                                 >  1:1193696/1‑100 (MQ=255)                               aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAct                                                                 >  2:1383519/1‑100 (MQ=255)                               aaaaTAAACTGGCATCGTCGTTACAGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                  >  2:1388303/1‑99 (MQ=255)                                   tAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTcc                                                             <  1:2541962/100‑1 (MQ=255)                                   tAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTcc                                                             <  1:2455204/100‑1 (MQ=255)                                    aaaCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAGCATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTcc                                                             <  2:1981345/99‑1 (MQ=255)                                     aaCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCgg                                                           >  1:181552/1‑100 (MQ=255)                                      aCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc                                                          <  1:200696/100‑1 (MQ=255)                                       cTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCa                                                         >  1:1949109/1‑100 (MQ=255)                                         ggCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCaa                                                        >  2:1435424/1‑99 (MQ=255)                                         ggCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAAt                                                       >  2:990365/1‑100 (MQ=255)                                           cATCGTCGTTACCTTTCACCGCAGCGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGTGCGTATCATCAACTCTCCGGCAATtc                                                     <  2:2165073/100‑1 (MQ=255)                                           cATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtc                                                     >  2:2543065/1‑100 (MQ=255)                                           cATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtc                                                     >  2:317031/1‑100 (MQ=255)                                            aTCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcg                                                    >  2:1705108/1‑100 (MQ=255)                                              cgtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtc                                                  >  1:1293814/1‑100 (MQ=255)                                               gtcgtTACCTTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCTATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtcg                                                 <  2:1381006/100‑1 (MQ=255)                                               gtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtcg                                                 >  1:907307/1‑100 (MQ=255)                                                   ttACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTg                                             <  1:1767656/100‑1 (MQ=255)                                                   ttACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTg                                             >  1:1853946/1‑100 (MQ=255)                                                    tACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGc                                            <  2:2382006/100‑1 (MQ=255)                                                      ccGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaa                                          >  1:424326/1‑100 (MQ=255)                                                      ccGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaa                                          >  2:1853898/1‑100 (MQ=255)                                                       cGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCATGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCATCGTCTGcaac                                         >  2:1188198/1‑100 (MQ=255)                                                        gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                        >  1:1837830/1‑100 (MQ=255)                                                          tCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaa                                      >  1:187724/1‑100 (MQ=255)                                                           cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa                                     >  2:200686/1‑100 (MQ=255)                                                           cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaa                                      <  2:270070/99‑1 (MQ=255)                                                             ccGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGa                                   >  2:1198188/1‑100 (MQ=255)                                                             ccGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGa                                   >  1:1717167/1‑100 (MQ=255)                                                             ccGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGa                                   >  2:1952010/1‑100 (MQ=255)                                                              cGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                  <  1:699080/100‑1 (MQ=255)                                                               gCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                  >  2:2412957/1‑99 (MQ=255)                                                                   ggCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAgg                             <  1:975246/100‑1 (MQ=255)                                                                   ggCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAgg                             >  1:2327713/1‑100 (MQ=255)                                                                     cGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGtt                           >  1:29009/1‑100 (MQ=255)                                                                      gTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttt                          >  1:1249396/1‑100 (MQ=255)                                                                      gTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttt                          >  1:1991503/1‑100 (MQ=255)                                                                       ttAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGtttt                         <  2:2274374/100‑1 (MQ=255)                                                                         aaaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                       <  1:1450728/100‑1 (MQ=255)                                                                         aaaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                       <  1:2330693/100‑1 (MQ=255)                                                                          aaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                      <  2:671589/100‑1 (MQ=255)                                                                           aaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCa                     >  1:1385778/1‑100 (MQ=255)                                                                            aCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                    >  2:2075041/1‑100 (MQ=255)                                                                            aCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                    <  2:1736584/100‑1 (MQ=255)                                                                            aCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                    <  1:2456004/100‑1 (MQ=255)                                                                             ccGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACt                   >  2:2413638/1‑100 (MQ=255)                                                                              cGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTg                  <  2:550553/100‑1 (MQ=255)                                                                              cGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTg                  <  2:510280/100‑1 (MQ=255)                                                                               gAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTgg                 <  1:1711452/100‑1 (MQ=255)                                                                                 aCATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg               <  1:2230653/100‑1 (MQ=255)                                                                                 aCATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg               <  1:1392214/100‑1 (MQ=255)                                                                                   aTGAGATCCTTCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg               <  2:1833415/98‑1 (MQ=255)                                                                                    tGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc            >  2:998244/1‑100 (MQ=255)                                                                                    tGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc            <  1:1705120/100‑1 (MQ=255)                                                                                    tGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc            >  1:1285635/1‑100 (MQ=255)                                                                                    tGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc            >  2:1838158/1‑100 (MQ=255)                                                                                     gtgaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCa           <  1:303962/98‑1 (MQ=255)                                                                                     gagaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCa           <  1:2289858/100‑1 (MQ=255)                                                                                      agaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCaa          >  2:1034334/1‑100 (MQ=255)                                                                                       gaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAAt         >  1:2212552/1‑100 (MQ=255)                                                                                           cTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccg     <  1:1937606/100‑1 (MQ=255)                                                                                              cGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGccgccg  <  1:245660/100‑1 (MQ=255)                                                                                                 |                                                                                                 GAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCG  >  NC_000913/179231‑179423

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