Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,506,0380AC42.9% 99.1 / 12.1 63F15V (TTT→GTT) ydcOBenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base A (28/8);  new base C (0/27);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 5.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGC  >  NC_000913/1505889‑1506185
                                                                                                                                                       |                                                                                                                                                   
ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGaaa                                                                                                                                                     <  2:608145/150‑1 (MQ=255)
ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAAc                                                                                                                                                    >  1:931886/1‑151 (MQ=255)
ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAAc                                                                                                                                                    >  2:239852/1‑151 (MQ=255)
    aaTGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAAcccc                                                                                                                                                 <  2:421441/150‑1 (MQ=255)
       gTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGcc                                                                                                                                              <  1:940100/150‑1 (MQ=255)
        tCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGcc                                                                                                                                              >  1:588545/1‑149 (MQ=255)
            agTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGc                                                                                                                                        >  1:1021639/1‑151 (MQ=255)
                    tACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCgtgggtgg                                                                                                                                >  1:790726/1‑151 (MQ=255)
                    tACGGACGCCCTTTGCCAGCCCCAGCCCCGTCACCCACCAAGAGTTTTGGGCGGTGGGGGCTCCGGCGACAATCGCTGTTTGCCAAATTATTGCCGCTGAATGGGCGTAGCAAATAAATCCCGCAAGAAACCCCGCCAGTAGCgtgggtg                                                                                                                                 <  1:820165/150‑1 (MQ=255)
                     aCTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCgtgggtgg                                                                                                                                >  2:451805/1‑150 (MQ=255)
                          cccccatccccccccccAGCGCCTTAACCCAGCAAGAGTTTGGTGCAGTGGTGCCTCGGGAAAAAATCGCTGTTTCCCAGATTATTGCCGCGGAAGTGGGTTAGCCAATTAATACCGCCAAAACCCCCGCCAGTAGCGTGGGTggag                                                                                                                              <  2:341784/134‑1 (MQ=255)
                            cctaTTGCCGGCCCCCGCGCCGCCACCAACCGAGAGATTTGGGGAGGGGTGGCTCGGACGAAAATCGCTGCTTGCCAGATTTTTGCCGCTGAGCGGGCGTGCCAAATTATTCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAg                                                                                                                        <  2:791879/148‑1 (MQ=255)
                             cctttGCCAGCCCCACCGCCGTCACCCAGCCAGAGTTTTGTGCAGTGGTGGCCCGGGCGACAATCCCTGCTTGCCAAATTATTGCCGCTGACCTGGCGTGGCCAATTAATCCCCCCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGa                                                                                                                       <  2:308822/148‑1 (MQ=255)
                               aTTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGaa                                                                                                                      >  1:266987/1‑150 (MQ=255)
                                ttcccACCCCCGGCGCCGTCACCCAGCCAGATTTTTGTGCGGTGGTGGCCCCGGCAACAATCGCTGCTTCCCAGATTATGGCCGCTGACCGGGCGTAGCAAATAAATCCCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGaaa                                                                                                                     <  2:642077/147‑1 (MQ=255)
                                tcgCCACCCCCAGCGCCGTCACCCAGCAAGAGATTTGTGCGGTGGTGCCTCGGGCGACAACCCCTGTTTGCCAGTTTTTCCCCGCGAAACTGGCGAAGCAAATTAAACCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGaa                                                                                                                      <  2:650826/147‑1 (MQ=255)
                                  gCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAAc                                                                                                                   >  1:401099/1‑150 (MQ=255)
                                   ccACCCCCAGCGCCGTCACCCACCAAGAGTTTTGTGCAGTGGTGGCCCCGGCGAAAATCGCTGTTTGCCAGATTATTGCCGCTGAACGGGCGTACCAAATAAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGaaaa                                                                                                                    <  2:504272/148‑1 (MQ=255)
                                             cgcCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                                                                                                       >  1:237150/1‑151 (MQ=255)
                                              gcCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                                                                                                       >  1:383545/1‑150 (MQ=255)
                                               ccttCACCCAGCCAGAGTTTTGGGCGGTGGGGCCCCCGGCGACAACCGCTGCTTGCAAGATTATTGCCGCTGACCGGGGGTAGCAAATAAACCCCCCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAg                                                                                                         <  1:363056/144‑1 (MQ=255)
                                                      ccagccagAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGc                                                                                              >  1:826048/1‑151 (MQ=255)
                                                       cagccagAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCg                                                                                             >  2:337228/1‑151 (MQ=255)
                                                            agatttttGGGCAGTGGTGGCCCCGGCAACAATCGTTTCTTGCCATATTATTGCCGCTGAACGGGCGTAGACAATTAATCCCGCCGGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat                                                                                         <  1:658898/145‑1 (MQ=255)
                                                             gattttggtgCAGGGGGGGCTCCGGCAACACTCGCTTCTCGCAAGATTATGGCCGATGAAGGGGCGTAGAAAATAAAACCCCCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtata                                                                                        <  1:366544/143‑1 (MQ=255)
                                                                    gggcgggggtggCTCCGGCAACATCCGGTGCTTGCCAAATTATTCCCGCTAAACGGGGGTGCCAAATAAAACCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAg                                                                                <  1:702709/144‑1 (MQ=255)
                                                                       cgggggtgcCTCCGCCGACAATCGCTCCTTCCCAGATTATTGCCGCTAAACGGGTGTACCAAATTAACCCCCCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc                                                                                  <  1:684868/142‑1 (MQ=255)
                                                                        ggggggggcccgggCGAAAATCGCTGCTGCCAAGATTATGCCCGCGGACCGGGGGAAGCAAATAAAACCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                                            <  1:940838/139‑1 (MQ=255)
                                                                          tggtggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCg                                                                          >  2:113402/1‑151 (MQ=255)
                                                                            gtggCCCCGGCCACAACCGCTGCTTGCCAGTTTTTTCCCGCTAATGTGCGTAACCAAATAAATCCCCCAATAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGt                                                                         <  1:145509/150‑1 (MQ=255)
                                                                                 tCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                   >  1:142653/1‑151 (MQ=255)
                                                                                  ccGGCGACAATCGCTGCTTGCCAGATTATGGCCGCTGAGCTGGGGTAGCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                   <  2:833104/150‑1 (MQ=255)
                                                                                    gggaaaatccccggctTGCCAGATTATGGCCGGTGAAGGGGGGTAGCAAATAAACCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                   <  1:417483/135‑1 (MQ=255)
                                                                                          aaTCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGtt                                                           >  1:791029/1‑150 (MQ=255)
                                                                                          aaTCGCGGTTTCCCAGATTATTGCCGCTGAGCTGGCGTAGCCAATAAACCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGtt                                                           <  1:639818/150‑1 (MQ=255)
                                                                                                tgctTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc                                                    >  1:997481/1‑151 (MQ=255)
                                                                                                tgctTCCCAGATTTTTGCCGCTAAACGGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATaa                                                       <  2:199010/148‑1 (MQ=255)
                                                                                                  ctTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcac                                                  >  2:845607/1‑151 (MQ=255)
                                                                                                      cgagtttatccccGTTGAACGGGCGTGGCAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt                                              <  1:259963/140‑1 (MQ=255)
                                                                                                      ccAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt                                              >  2:249681/1‑151 (MQ=255)
                                                                                                      caaGATTATTGCCGTTAAACGGGCGTAGCAAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                                               <  2:364176/148‑1 (MQ=255)
                                                                                                       cAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTg                                             >  1:819618/1‑151 (MQ=255)
                                                                                                        agtttattGCCGCTGAACTGGCGTAGCCAATTAATCCCGCAAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt                                              <  2:398088/146‑1 (MQ=255)
                                                                                                            tattGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATaa                                         >  2:517229/1‑150 (MQ=255)
                                                                                                            tattGCCGCTAAACTGGCGTAGCAAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATaa                                         <  2:829340/150‑1 (MQ=255)
                                                                                                             attGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATaaa                                        >  1:294263/1‑150 (MQ=255)
                                                                                                                 ccGCTGAACGGGCGTAGCAAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACtt                                     <  2:264078/149‑1 (MQ=255)
                                                                                                                    cTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt                                >  1:523427/1‑151 (MQ=255)
                                                                                                                    cTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt                                >  1:717489/1‑151 (MQ=255)
                                                                                                                     tGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGTATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGttt                               >  2:905606/1‑151 (MQ=255)
                                                                                                                       acggggggTAGCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                              <  1:66953/143‑1 (MQ=255)
                                                                                                                       aaCGGGCGTGGCAAATAAACCCCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                              <  1:910464/150‑1 (MQ=255)
                                                                                                                              gTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCggagg                      >  2:977586/1‑151 (MQ=255)
                                                                                                                              gTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCggagg                      >  1:821855/1‑151 (MQ=255)
                                                                                                                                acccAATAAACCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGgaggag                    <  1:404432/149‑1 (MQ=255)
                                                                                                                                     aaataaacccccccAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTc                  <  1:462124/138‑1 (MQ=255)
                                                                                                                                     aaTAAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTc                  <  1:50314/149‑1 (MQ=255)
                                                                                                                                       taaaTCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTc                  <  1:211634/145‑1 (MQ=255)
                                                                                                                                              ccGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGaaa         <  1:200809/150‑1 (MQ=255)
                                                                                                                                               cccaaGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGaaa         <  2:237154/145‑1 (MQ=255)
                                                                                                                                               cGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAt       <  2:1011166/151‑1 (MQ=255)
                                                                                                                                               cGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAt       >  2:834878/1‑151 (MQ=255)
                                                                                                                                                   aaaaaCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAtct     <  2:461098/147‑1 (MQ=255)
                                                                                                                                                   aGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAtcc     <  2:461140/149‑2 (MQ=255)
                                                                                                                                                    gAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATctcgc  >  2:293326/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGC  >  NC_000913/1505889‑1506185

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