Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,506,0380AC41.6% 152.8 / 11.7 77F15V (TTT→GTT) ydcOBenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base A (40/5);  new base C (0/32);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 9.45e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGC  >  NC_000913/1505889‑1506185
                                                                                                                                                      |                                                                                                                                                    
ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGaaa                                                                                                                                                      <  2:951118/150‑1 (MQ=255)
   tAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAAcccc                                                                                                                                                 >  1:1100968/1‑151 (MQ=255)
     aTGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGc                                                                                                                                               >  2:70938/1‑151 (MQ=255)
      tGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGcc                                                                                                                                              >  1:879880/1‑151 (MQ=255)
      tGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGcc                                                                                                                                              >  1:1079413/1‑151 (MQ=255)
       gtaagagTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGcc                                                                                                                                              >  2:863286/4‑150 (MQ=255)
       gTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGcc                                                                                                                                              <  2:243928/150‑1 (MQ=255)
              tCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCg                                                                                                                                       >  1:1085770/1‑150 (MQ=255)
                      cggccGCCCTTTGCCAGCCCCAGCCCCGTCACCCAGCCAGAGATTTGTGAAGGGGGGGCTCCGGCAACAACCGTTGCTTGCCAGATTATTCCGCTTGACCTGGGTTACCAAATTAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTgga                                                                                                                               <  1:1054256/146‑1 (MQ=255)
                         aCGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggagga                                                                                                                            >  2:280600/1‑150 (MQ=255)
                           gcccttggccccccccGGCGCCGTCACCCGGCCAGAGTTTGGTGGGGGGGTGCCTCGGGCGACAATCGTTGCTTGCCAGATTATGCCCGCTGAACGGGGGTAGCAAATAAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGaa                                                                                                                           <  2:706726/137‑1 (MQ=255)
                              cATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGa                                                                                                                       >  1:47029/1‑150 (MQ=255)
                                ttcccggCCCCAGCGCCGTCACCCAGCCAGAGATTTGTGCAGTGGTGGCCCCGGCAACAACCGCTGTTTGCCAGATTATTGCCGCTGAACGGGGGTAGCCAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGaa                                                                                                                      <  1:435344/143‑1 (MQ=255)
                                   ccacccccgccccggTCACCCAGCGAGATTTTTGTGAGGGGGGGGCTCCGGGAAAAACCGCTGTTTCCAAGATTTTTCCGGCGGACCGGGGGTAGCAAATTAATACCCCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGaa                                                                                                                      <  2:1194497/132‑1 (MQ=255)
                                   ccAGCCCCGGCGCCCTCACCCAGCCAGATTTTTGTGCAGTGGGGGCTCCGGCGAAAACCGCTGTTTGCCAGATTATTGCCGTTGAACGGGGGTAGCAAATAAATCCCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACa                                                                                                                  <  2:1162878/150‑1 (MQ=255)
                                    cAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAAc                                                                                                                   >  1:700506/1‑148 (MQ=255)
                                    cAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAAc                                                                                                                   >  2:11148/1‑148 (MQ=255)
                                        cccAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGc                                                                                                             >  2:842001/1‑150 (MQ=255)
                                         ccggcgcCGTCACCCAGCAAGAGTTTTGGGCGGTGGTGGCTCCGGCGACATTCTTTGTTGGCCAGATTATTGCCGCTAAACGGGCGTACCAAATTAATCCCCCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCAt                                                                                                           <  1:843919/148‑1 (MQ=255)
                                             ccccGTCACCCAGCCAGAGTTTGGTGCAGGGGTGGCTCCGGCGAAAATCGCTGTTTCCAAAATTATGGCCGCGGAACTGGCGTAGCAAATAAACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCAt                                                                                                           <  1:331259/145‑1 (MQ=255)
                                               ccttCACCCAGCCAGAGTTTTGTGCTGTGGGGGCTCGGGCAACATTCGCTGCTTGCCAGTTTTTTGCCGCTGAACGGGGGTAGCAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAg                                                                                                         <  2:627848/144‑1 (MQ=255)
                                               ccGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCAc                                                                                                     >  1:607803/1‑151 (MQ=255)
                                                  tCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTc                                                                                                  >  1:400363/1‑151 (MQ=255)
                                                    acccgcccagAATTTTGTGGAGTGGTGGCTCGGGCAAAAATCGCTTCTTGCAAGATTATTGCCGCTGAATGGGCGTACCAAATTAACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTc                                                                                                  <  1:71443/143‑1 (MQ=255)
                                                    aTCCAGCCAGAGTTTTGTGCAGTGGGGGCTCCGGCAACAACCGCTGCTTGCCAGATTATTCCCGCGGAACGGGGGTAGCAAATAAATCCCGCAAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt                                                                                                   <  1:366175/148‑1 (MQ=255)
                                                      ccggcaagggttttGTGGGGTGGTGGCTCGGGCAAAAATCGCTGCTTCCAATATTATTGCCGTGGAACGGGCGTAGCAAATTAATCCCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCg                                                                                                 <  2:950810/137‑1 (MQ=255)
                                                      ccagccagAGATTTGTGCAGGGGTGGCCCCGGCGACAATCGCTGTTTGCCAGATTATTGCCGCTGAACTGGCGTAGCAAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGt                                                                                                <  1:202707/149‑1 (MQ=255)
                                                            agagTTTTGTGCAGTGGTGGCCCCGGCGAAAATCGCTGCTTGCCAAATTATTGCCGCGGAACGGGCGTACCCAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat                                                                                         <  1:174744/150‑1 (MQ=255)
                                                               gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                                                     >  1:1045746/1‑151 (MQ=255)
                                                               gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                                                     >  1:495176/1‑151 (MQ=255)
                                                                 tttGTGCAGGGGTGGCTCCGGCGAAAATCGCTGCTTGCAAGATTATTGCCGCTGAACTGGGGTAGCAAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAg                                                                                       <  1:1037114/147‑1 (MQ=255)
                                                                    gtgCGGTGGGGGCTCCGGCGACAATCGCTGTTTGCCAGATTATTGCCGCGGAACTGGCGTAGCAAATAAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc                                                                                  <  1:796490/149‑1 (MQ=255)
                                                                      gcggtggtgGCCCCGGCGACAATCGCTGCTTGCCAGATTATGGCCGCTGACCGGGCGTAGCAAATAAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAgg                                                                               <  1:21441/147‑1 (MQ=255)
                                                                         gaggtggCCCCGGCAAAAACCGCGGCTCCCAAATTTATCGCCGATGACTGGGCGTAGCAAATAAACCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc                                                                                  <  2:252858/142‑1 (MQ=255)
                                                                          tggtggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCg                                                                          >  1:676532/1‑151 (MQ=255)
                                                                                 tCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                   >  1:695800/1‑151 (MQ=255)
                                                                                 tCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                   >  1:718521/1‑151 (MQ=255)
                                                                                    ggCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt                                                                >  1:778302/1‑151 (MQ=255)
                                                                                    ggCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt                                                                >  1:785840/1‑151 (MQ=255)
                                                                                    ggCAACAATCGCTGCTTGCCAGATTATTGCCGCTGAACGGGCGTAGCAAATTAATCCCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAt                                                                    <  1:1063661/147‑1 (MQ=255)
                                                                                     ggaaaaactggcgctTGCCATATTTTTGCCGGTAACCGGGGGGAGCAATTTAAACGCGCCAGACCCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAttt                                                                  <  2:488826/136‑1 (MQ=255)
                                                                                      cGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTcc                                                              >  1:595212/1‑151 (MQ=255)
                                                                                        aaaaTCGCTGTTGCCCATATTATTGCCGCTGACCGGGCGTAGCAATTAAATCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTcc                                                              <  1:460225/147‑1 (MQ=255)
                                                                                        aCAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGt                                                            >  2:1005757/1‑151 (MQ=255)
                                                                                           aTCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTtat                                                         >  2:843703/1‑151 (MQ=255)
                                                                                              gctctttGCCAGATTTTTGCCGCGAAAGGGGCGTAGCAAATTAATCCCCCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATaa                                                       <  2:329462/145‑1 (MQ=255)
                                                                                                   ttGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcac                                                  >  1:302278/1‑150 (MQ=255)
                                                                                                   ttGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcac                                                  >  2:798982/1‑150 (MQ=255)
                                                                                                    tGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca                                                 >  1:67004/1‑150 (MQ=255)
                                                                                                       cAGTTTATTCCCGCTAACCGGGCGTAGCAAATAAATCCCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca                                                 <  1:773861/147‑1 (MQ=255)
                                                                                                        aGATTATTGCCGCTGAACTGGCGTAGCCAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGa                                            <  1:1191429/151‑1 (MQ=255)
                                                                                                           ttattGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATa                                          >  2:1074317/1‑150 (MQ=255)
                                                                                                           ttattCCCGCTGACCGGGCGGAACAAATAAATACCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGAt                                           <  1:369884/149‑1 (MQ=255)
                                                                                                            tattGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATaa                                         >  1:783115/1‑150 (MQ=255)
                                                                                                             attGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAAc                                       >  2:281737/1‑151 (MQ=255)
                                                                                                                 ccGCTGACCGGGCGTAGCAAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTcc                                   <  1:1050542/151‑1 (MQ=255)
                                                                                                                   gCTGAACGGGGGTAGCAAATTAATCCCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCg                                  <  2:378414/150‑1 (MQ=255)
                                                                                                                    cTGAACGGGCGTAGCAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACtt                                     <  1:698924/146‑1 (MQ=255)
                                                                                                                       gAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                              >  1:1168758/1‑151 (MQ=255)
                                                                                                                        accgggggtaccaaaTAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                              <  1:717342/137‑1 (MQ=255)
                                                                                                                        aacggggttAGCAAATAAACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCg                                  <  2:123894/138‑1 (MQ=255)
                                                                                                                        aaCTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTg                             >  1:337426/1‑151 (MQ=255)
                                                                                                                        aaCTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTg                             >  1:667541/1‑151 (MQ=255)
                                                                                                                         aCTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc                            >  1:906714/1‑151 (MQ=255)
                                                                                                                          cTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGcc                           >  2:505755/1‑151 (MQ=255)
                                                                                                                            gggttAGCAAATAAACCCCACCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc                            <  2:667287/144‑1 (MQ=255)
                                                                                                                                 aGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCggagg                      >  2:1112873/1‑149 (MQ=255)
                                                                                                                                   caaaTTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGt                   <  1:196195/148‑1 (MQ=255)
                                                                                                                                     aaTTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCg                 <  1:317241/150‑1 (MQ=255)
                                                                                                                                     aaTTAATACCGCCAGAAACCCCGCCAGAAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCg                 >  2:672622/1‑150 (MQ=255)
                                                                                                                                     aaTTAATACCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGt                   <  1:1032526/148‑1 (MQ=255)
                                                                                                                                         aaacccGCAAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCAt              <  2:450061/145‑1 (MQ=255)
                                                                                                                                          aacccccaaGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCAt              <  2:929617/140‑1 (MQ=255)
                                                                                                                                              ccccAGAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGaaa         <  2:943617/148‑1 (MQ=255)
                                                                                                                                              cGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAt       >  2:842679/1‑151 (MQ=255)
                                                                                                                                               gCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAtc      >  1:931686/1‑151 (MQ=255)
                                                                                                                                                 cAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAtctc    >  1:367154/1‑151 (MQ=255)
                                                                                                                                                  aGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATctcg   >  1:542970/1‑151 (MQ=255)
                                                                                                                                                    aaaCCCCCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATctcgc  <  1:810414/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGC  >  NC_000913/1505889‑1506185

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