Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,350,8440AC38.8% 120.4 / 16.3 67L728V (TTA→GTA) arcBaerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base A (33/8);  new base C (0/26);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 8.30e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATC  >  NC_000913/3350700‑3350983
                                                                                                                                                  |                                                                                                                                           
tGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCAc                                                                                                                                       >  1:324808/1‑151 (MQ=255)
 ggCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCAc                                                                                                                                       <  1:135011/150‑1 (MQ=255)
   cTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTg                                                                                                                                     >  1:46814/1‑150 (MQ=255)
     tttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGa                                                                                                                                    <  1:515824/149‑1 (MQ=255)
        gCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccg                                                                                                                               >  2:1085902/1‑151 (MQ=255)
        gCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCAGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                <  1:904321/150‑1 (MQ=255)
         ccacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccg                                                                                                                               <  2:96937/150‑1 (MQ=255)
          cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                             >  1:977179/1‑151 (MQ=255)
                 gCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCtt                                                                                                                      >  1:321126/1‑151 (MQ=255)
                   tttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                     >  2:710033/1‑150 (MQ=255)
                      cAGCACTTCTACGTGGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCCCCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTCACCCATCTGTTGCAGGTGGCGTACCCCCACTGCCCCCGCCGCACCTTTAAtt                                                                                                                 <  1:982948/151‑1 (MQ=255)
                       aGCACTTCTCCTTCGTGACGCCCCTCTTCTTTCATCTCTTCAACCCATCCCCCGACTTTTTCTTCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGACCCAGCTTTTGCAGAGGGCGTACCCCCACTGACCCCGCCGCACCTTTAAtt                                                                                                                 <  1:3874/150‑1 (MQ=255)
                          aCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTAt                                                                                                             >  1:32160/1‑151 (MQ=255)
                                 cgtcggGCCGCCCCTCTTCTTTCATCTCTCCAACCCTTCCCCCGTCTTTATCTCCCCGGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCACCGTTTGCAGTTGGCGTACCCCCACTAACCCCGCCGCACCTTTAATTTTATGt                                                                                                           <  1:455105/146‑1 (MQ=255)
                                  gtcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcc                                                                                                     >  1:722451/1‑151 (MQ=255)
                                   tcgtGACGCCCCTCTTCTTTCTTCCCTTCAATCCACTCCCCGACTTTATCTCCCCAGGCCGAAATTTCAGTATATTAAATTTGCTGACCCATTGTTTTCAGATGGCGTAACCCCACGGACCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                       <  2:788043/148‑1 (MQ=255)
                                         ccccccttttcttTCATCTCTTCATCCCATTCCCCTACTTTATCTCCCCTGGCCGGAAGGACGGTGTACTAAATTGGCCCACCCATTTTTTGTAGTTGCCGCAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa                                                                                              <  2:46801/143‑1 (MQ=255)
                                            cccttttcttTCCCCTCTTAAACCCATCCCCCGCTTTTTTCTCCCCGGGCGGGAGGTCAAGGGGACTGAATTGCCTGCCCCATTTTTTGCAGAGGGCGTACCCCAACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca                                                                                               <  2:126746/142‑1 (MQ=255)
                                            ccctctTCTTTCTTCTCTTAAACCCATCCCCCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTTATTTTGCTGCCCCATTTGTTGCAGATGGCGTAACCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                           <  2:378754/149‑1 (MQ=255)
                                            cACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                           >  1:121516/1‑151 (MQ=255)
                                             aCTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                                                                          >  1:216685/1‑151 (MQ=255)
                                              ctctTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCt                                                                                         >  1:134316/1‑151 (MQ=255)
                                                 tttttcctcctttacAACCCATCCCCCGCCGTTATCTCCCCGGGCGGGAAGGTCAGGAGACGGATTTGCCTGCCCCAGTTGTTGCAGATGGCGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                           <  2:922035/132‑1 (MQ=255)
                                                 ttcttTCATCTCTTCAATCCATTCACCGACGATATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt                                                                                      >  1:242519/1‑151 (MQ=255)
                                                     ttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                                  >  2:650919/1‑151 (MQ=255)
                                                         tctcTTCACTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                              >  2:696259/1‑151 (MQ=255)
                                                          ctctTCAATCCATCCCCCGACGTTATCTCCCCAGGCGGGAGGGTCAGGAGACTAAATTTGCTGCCCCAGCTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATcc                                                                                <  1:702297/148‑1 (MQ=255)
                                                          ctctTCAACCCATCCCCCGACTTTATCTCCCCGGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGAGGGCATACCCCCACGGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                                  <  1:113209/146‑1 (MQ=255)
                                                               cAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTCTCCTGCGCCg                                                                         >  2:130723/1‑150 (MQ=255)
                                                                aacccATTCCCCGACTTTATCTCCCCAGCCCGGAAGGTCAGGAGCCTGATTTGCCTACCCCATTTTTTGCAGATGCCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCt                                                                               <  2:359589/140‑1 (MQ=255)
                                                                aacccATCCCCCGCCGTTATCTCCCCAGCCCGGAGGGTCAGGAGACAGAATTTGCTGCCCCAGCTGTTGAAGATGGCGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                             <  1:742261/142‑1 (MQ=255)
                                                                aacccATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTTTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGt                                                                        <  1:1086255/147‑1 (MQ=255)
                                                                 aTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                       >  2:132249/1‑150 (MQ=255)
                                                                  tCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                       >  2:13201/1‑149 (MQ=255)
                                                                   ccttccACCGACTTTCTTTCCCCAGGCCGGAAGGTCAGGAGACTAATTTTCCTGCCCCAGGTGTTGAGGAGGGCATACCCCCCCGACCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                       <  2:799312/143‑1 (MQ=255)
                                                                   ccATCCCCCGACTTTATCTCCCCGGGCCGGAGGGTCAGGAGACTGAATTGGCTGCCCCGGTTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAg                                                                     <  1:101116/150‑1 (MQ=255)
                                                                   ccATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAg                                                                     <  1:907782/150‑1 (MQ=255)
                                                                    cATCCACCGACTTTATCTTCCCAGGCCGGAAGGTCGGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAg                                                                     <  1:931533/149‑1 (MQ=255)
                                                                     aTTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                                  >  1:1047754/1‑151 (MQ=255)
                                                                      ttCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                 >  1:331091/1‑151 (MQ=255)
                                                                        cACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                >  2:163501/1‑150 (MQ=255)
                                                                          ccgccgTTTTCTCCCCAGGCCGGAAGGGCAGGAGCCTGAATTTGCTGCCCCAGTTTTTGCAGTTGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                 <  2:674354/147‑1 (MQ=255)
                                                                           cgcctttttCTTCCGAGGCGGGAAGGTGAGGAGACTGAATTTGCGGACCCAGCTGTTGCAGAGGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACt                                                             <  1:769114/142‑1 (MQ=255)
                                                                               gTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                        >  2:1033602/1‑151 (MQ=255)
                                                                                ttctctcCCAAGGCAGAAGGGCCAGGGCACGGAATTTGCTACCCAAGTTTTGGAAGTGGTCGTACCCCCATTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCa                                                          <  2:333781/145‑1 (MQ=255)
                                                                                 tATCTCCCCGGGCCGAAGGGTCGGGCAAATGAATTTGCTCCCCCATCTGTTGCAGATGGCGTAACCCCCCGCACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAg                                                         <  2:192700/148‑1 (MQ=255)
                                                                                     ttCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg                                                  >  1:1001701/1‑151 (MQ=255)
                                                                                       cccAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt                                                   <  1:41112/148‑1 (MQ=255)
                                                                                       cccAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                                >  2:603912/1‑151 (MQ=255)
                                                                                         cAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                              >  1:846039/1‑151 (MQ=255)
                                                                                               ggAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCg                                        >  1:1072332/1‑151 (MQ=255)
                                                                                                   ggTCGGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcat                                    <  1:407019/151‑1 (MQ=255)
                                                                                                        aGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCatcatc                                >  2:909247/1‑151 (MQ=255)
                                                                                                         ggAGACTGAATTTGCTGCCCCAGTTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat                                 <  2:32151/149‑1 (MQ=255)
                                                                                                         ggAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtct                               >  2:563522/1‑151 (MQ=255)
                                                                                                           agaCTGATTTTGCTCCCCCAGCTGTTGCAGATGCCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt                              <  2:121481/150‑1 (MQ=255)
                                                                                                            gaCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct                            >  1:854069/1‑151 (MQ=255)
                                                                                                                gAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa                         >  2:724909/1‑150 (MQ=255)
                                                                                                                   tttGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCa                          <  1:489315/146‑1 (MQ=255)
                                                                                                                   tttGCTGCCCCAGTTGTTGAAGATGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttc                             <  2:553501/143‑1 (MQ=255)
                                                                                                                         gccccAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                    <  1:240529/144‑1 (MQ=255)
                                                                                                                           cccAGCTGTTGCAGATGCCTTACCCCCCCGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAcc             <  2:1058474/151‑1 (MQ=255)
                                                                                                                             cAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCg           >  1:1102438/1‑151 (MQ=255)
                                                                                                                                cTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg        >  1:249128/1‑151 (MQ=255)
                                                                                                                                     gCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGAt   >  1:338222/1‑151 (MQ=255)
                                                                                                                                      cAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGa    <  1:434486/149‑1 (MQ=255)
                                                                                                                                       cgATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGCGATc  >  2:621535/2‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
TGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATC  >  NC_000913/3350700‑3350983

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