Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	734,227	C→T	5.1%	pseudogene (95/1521 nt)	rhsO →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	734,227	0	C	T	5.1%	238.0 / 2.4	78	pseudogene (95/1521 nt)	rhsO	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base T (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAG  >  NC_000913/734080‑734369                                                                                                                                                    |                                                                                                                                               gtaagCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGcc                                                                                                                                                <  6:345374/148‑1 (MQ=255)      cgcATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATa                                                                                                                                           <  7:213071/148‑1 (MQ=255)       gcATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATaa                                                                                                                                          <  6:306702/148‑1 (MQ=255)           tataAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGACTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCg                                                                                                                                       >  1:415715/1‑147 (MQ=255)           tataAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACACGGCAGGTAACCGACTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGTCCGGATAAACg                                                                                                                                       >  1:421899/1‑147 (MQ=255)            ataAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGt                                                                                                                                      <  6:124307/147‑1 (MQ=255)            ataAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGt                                                                                                                                      <  8:96367/147‑1 (MQ=255)             taAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTa                                                                                                                                     <  7:76554/147‑1 (MQ=255)              aaaCAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCCCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTa                                                                                                                                     >  1:318097/1‑146 (MQ=255)              aaaCAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCCCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTa                                                                                                                                     <  2:318097/146‑1 (MQ=255)                   aataaatTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCg                                                                                                                              >  7:210966/1‑148 (MQ=255)                                 aTTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATc                                                                                                                <  3:302650/148‑1 (MQ=255)                                     ctctTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCttt                                                                                                             <  8:210966/147‑1 (MQ=255)                                     ctctTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCttt                                                                                                             <  5:158915/147‑1 (MQ=255)                                     ctctTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTa                                                                                                            <  2:51623/148‑1 (MQ=255)                                        ttattattaTATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCTTGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTAAATTTACCg                                                                                                         >  8:120096/1‑148 (MQ=255)                                          attattatATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTAcc                                                                                                          <  3:112670/145‑1 (MQ=255)                                          attattatATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTAcc                                                                                                          <  1:90226/145‑1 (MQ=255)                                             attatATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTAt                                                                                                     >  5:358298/1‑147 (MQ=255)                                                   tGTAACCTGGGCATTGATATCCCGTATGCCACAGCCCCGGCCGGTAACCGCCTGCCCGACCCGGCGCTGCACCCGGCCAGCACCCTCAGCCTGTGGCCGGATGACCGTATCGCCCGTGCCGCGCACTCTCTTTACCGTTATGACCGt                                                                                               >  5:190235/1‑147 (MQ=14)                                                   tGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCTGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGt                                                                                               <  8:399609/147‑1 (MQ=255)                                                      aaCCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACg                                                                                           >  6:171826/1‑148 (MQ=255)                                                        ccTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGc                                                                                         <  1:188004/148‑1 (MQ=40)                                                             gCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt                                                                                    <  7:212340/148‑1 (MQ=35)                                                                 tGATATCCCGTATGCCACAGACCCGGCAGGTACCCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACgcgc                                                                                                                      <  4:392875/110‑1 (MQ=18)                                                                 tGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTg                                                                                                                            >  7:250241/1‑104 (MQ=33)                                                                 tGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTg                                                                                                                            <  8:250241/104‑1 (MQ=33)                                                                 tGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACgcgc                                                                                                                      >  3:392875/1‑110 (MQ=33)                                                                 tGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACg                                                                                >  5:358173/1‑148 (MQ=33)                                                                 tGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACg                                                                                >  3:329953/1‑148 (MQ=33)                                                                  gATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTc                                                                                              >  6:342909/1‑133 (MQ=32)                                                                  gATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTc                                                                                              <  5:342909/133‑1 (MQ=32)                                                                                 acaGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGCCCTCa                                                                >  7:229701/1‑148 (MQ=17)                                                                                  caGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt                                                               >  4:410111/1‑148 (MQ=32)                                                                                  caGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCGTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt                                                               >  1:17694/1‑148 (MQ=17)                                                                                   aGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATc                                                              >  5:232618/1‑148 (MQ=32)                                                                                      cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGgag                                                                             <  1:311880/130‑1 (MQ=32)                                                                                      cccGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGgag                                                                             >  2:311880/1‑130 (MQ=32)                                                                                        cGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                                                         <  6:358173/148‑1 (MQ=32)                                                                                           cAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAgg                                                      <  4:329953/148‑1 (MQ=32)                                                                                            aGGTAACCGCCGGCCCGACCCGGAGCGGCACCCGGACAGCACCCTCAGCCTGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTAACCGGTATGCCCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAggg                                                     <  6:190235/148‑1 (MQ=11)                                                                                                       tGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                                          >  7:166170/1‑148 (MQ=32)                                                                                                        gCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGa                                                                     >  2:59620/1‑120 (MQ=32)                                                                                                        gCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGa                                                                     <  1:59620/120‑1 (MQ=32)                                                                                                        gCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                                                         >  8:97159/1‑132 (MQ=32)                                                                                                        gCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                                                         <  7:97159/132‑1 (MQ=32)                                                                                                           cGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGgat                                      >  1:391289/1‑148 (MQ=32)                                                                                                            gACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGgat                                      >  4:100297/1‑147 (MQ=17)                                                                                                            gACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGgat                                      >  2:320268/1‑147 (MQ=17)                                                                                                             aCCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAggg                                                     <  8:387670/131‑1 (MQ=32)                                                                                                             aCCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAggg                                                     >  7:387670/1‑131 (MQ=32)                                                                                                             aCCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGgatga                                    <  6:358298/148‑1 (MQ=32)                                                                                                              cccGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGgatgat                                   >  5:409300/1‑148 (MQ=32)                                                                                                               ccGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGttgatg                                  <  6:232618/148‑1 (MQ=17)                                                                                                                cGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGgatgat                                   >  8:224461/1‑146 (MQ=32)                                                                                                                cGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGAtgatga                                 <  7:120096/148‑1 (MQ=32)                                                                                                                cGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGAtgatga                                 <  6:409300/148‑1 (MQ=32)                                                                                                                 ggAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAg                                <  2:415715/148‑1 (MQ=32)                                                                                                                 ggAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACAGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAg                                <  2:421899/148‑1 (MQ=17)                                                                                                                  gAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgc                               >  1:229351/1‑148 (MQ=32)                                                                                                                    gCTGCACCCGGACAGCACCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgcgc                             <  1:320268/148‑1 (MQ=17)                                                                                                                    gCTGCACCCGGACAGCACCCTCAGCATGTGGCTGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgcgc                             <  3:100297/148‑1 (MQ=17)                                                                                                                           ccGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcacccacc                      <  2:229351/148‑1 (MQ=34)                                                                                                                                cAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTAc                 <  2:17694/148‑1 (MQ=34)                                                                                                                                  gCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa               >  6:432079/1‑148 (MQ=34)                                                                                                                                    aCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCCCGGATGATGAGCGCACCCCCCGGTACCAt              >  8:384825/1‑147 (MQ=17)                                                                                                                                    aCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCCCCGGTACCAt              >  2:444908/1‑147 (MQ=17)                                                                                                                                    aCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCAt              >  2:219035/1‑147 (MQ=32)                                                                                                                                      ccTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                                                         <  7:406382/102‑1 (MQ=32)                                                                                                                                      ccTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                                                         >  8:406382/1‑102 (MQ=32)                                                                                                                                      ccTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc           >  8:18834/1‑148 (MQ=32)                                                                                                                                       cTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACg          <  8:229701/148‑1 (MQ=32)                                                                                                                                         cAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGa         >  5:371046/1‑147 (MQ=32)                                                                                                                                         cAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGa         >  6:173436/1‑147 (MQ=32)                                                                                                                                            cATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGt     <  2:391289/148‑1 (MQ=32)                                                                                                                                             aTGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTc    <  7:18834/148‑1 (MQ=32)                                                                                                                                               gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAg  >  4:70523/1‑148 (MQ=32)                                                                                                                                               gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAg  >  2:218777/1‑148 (MQ=32)                                                                                                                                                tgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAg  <  1:219035/147‑1 (MQ=32)                                                                                                                                                    |                                                                                                                                               GTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATTGATATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAG  >  NC_000913/734080‑734369

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