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_000913734,2270CT 5.1% 238.0 / 2.4 78pseudogene (95/1521 nt)rhsOpseudogene, 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: