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

CCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGC  >  NC_000913/1505889‑1506185
                                                                                                                                                       |                                                                                                                                                   
ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGaaa                                                                                                                                                     <  2:608145/150‑1 (MQ=255)
ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAAc                                                                                                                                                    >  1:931886/1‑151 (MQ=255)
ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAAc                                                                                                                                                    >  2:239852/1‑151 (MQ=255)
    aaTGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAAcccc                                                                                                                                                 <  2:421441/150‑1 (MQ=255)
       gTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGcc                                                                                                                                              <  1:940100/150‑1 (MQ=255)
        tCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGcc                                                                                                                                              >  1:588545/1‑149 (MQ=255)
            agTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGc                                                                                                                                        >  1:1021639/1‑151 (MQ=255)
                    tACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCgtgggtgg                                                                                                                                >  1:790726/1‑151 (MQ=255)
                    tACGGACGCCCTTTGCCAGCCCCAGCCCCGTCACCCACCAAGAGTTTTGGGCGGTGGGGGCTCCGGCGACAATCGCTGTTTGCCAAATTATTGCCGCTGAATGGGCGTAGCA‑AATAAATC‑CCGCAAGAAACCCCGCCAGTAGCgtgggtg                                                                                                                                 <  1:820165/150‑1 (MQ=255)
                     aCTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCgtgggtgg                                                                                                                                >  2:451805/1‑150 (MQ=255)
                          cccccatccccccccccAGCGCCTTAACCCAGCAAGAGTTTGGTGCAGTGGTGCCTCGGGAAAAAATCGCTGTTTCCCAGATTATTGCCGCGGAAGTGGGTTAGCC‑AATTAATA‑CCGCCAAAACCCCCGCCAGTAGCGTGGGTggag                                                                                                                              <  2:341784/134‑1 (MQ=255)
                            cctaTTGCCGGCCCCCGCGCCGCCACCAACCGAGAGATTTGGGGAGGGGTGGCTCGGACGAAAATCGCTGCTTGCCAGATTTTTGCCGCTGAGCGGGCGTGCCA‑AATTATTC‑CCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAg                                                                                                                        <  2:791879/148‑1 (MQ=255)
                             cctttGCCAGCCCCACCGCCGTCACCCAGCCAGAGTTTTGTGCAGTGGTGGCCCGGGCGACAATCCCTGCTTGCCAAATTATTGCCGCTGACCTGGCGTGGCC‑AATTAATC‑CCCCCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGa                                                                                                                       <  2:308822/148‑1 (MQ=255)
                               aTTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGaa                                                                                                                      >  1:266987/1‑150 (MQ=255)
                                ttcccACCCCCGGCGCCGTCACCCAGCCAGATTTTTGTGCGGTGGTGGCCCCGGCAACAATCGCTGCTTCCCAGATTATGGCCGCTGACCGGGCGTAGCA‑AATAAATC‑CCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGaaa                                                                                                                     <  2:642077/147‑1 (MQ=255)
                                tcgCCACCCCCAGCGCCGTCACCCAGCAAGAGATTTGTGCGGTGGTGCCTCGGGCGACAACCCCTGTTTGCCAGTTTTTCCCCGCGAAACTGGCGAAGCA‑AATTAAAC‑CCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGaa                                                                                                                      <  2:650826/147‑1 (MQ=255)
                                  gCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAAc                                                                                                                   >  1:401099/1‑150 (MQ=255)
                                   ccACCCCCAGCGCCGTCACCCACCAAGAGTTTTGTGCAGTGGTGGCCCCGGCGAAAATCGCTGTTTGCCAGATTATTGCCGCTGAACGGGCGTACCA‑AATAAATC‑CCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGaaaa                                                                                                                    <  2:504272/148‑1 (MQ=255)
                                             cgcCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                                                                                                       >  1:237150/1‑151 (MQ=255)
                                              gcCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                                                                                                       >  1:383545/1‑150 (MQ=255)
                                               ccttCACCCAGCCAGAGTTTTGGGCGGTGGGGCCCCCGGCGACAACCGCTGCTTGCAAGATTATTGCCGCTGACCGGGGGTAGCA‑AATAAACC‑CCCCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAg                                                                                                         <  1:363056/144‑1 (MQ=255)
                                                      ccagccagAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGc                                                                                              >  1:826048/1‑151 (MQ=255)
                                                       cagccagAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCg                                                                                             >  2:337228/1‑151 (MQ=255)
                                                            agatttttGGGCAGTGGTGGCCCCGGCAACAATCGTTTCTTGCCATATTATTGCCGCTGAACGGGCGTAGAC‑AATTAATC‑CCGCCGGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat                                                                                         <  1:658898/145‑1 (MQ=255)
                                                             gattttggtgCAGGGGGGGCTCCGGCAACACTCGCTTCTCGCAAGATTATGGCCGATGAAGGGGCGTAGAA‑AATAAAAC‑CCCCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtata                                                                                        <  1:366544/143‑1 (MQ=255)
                                                                    gggcgggggtggCTCCGGCAACATCCGGTGCTTGCCAAATTATTCCCGCTAAACGGGGGT‑GCC‑AAATAAAACCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAg                                                                                <  1:702709/144‑1 (MQ=255)
                                                                       cgggggtgcCTCCGCCGACAATCGCTCCTTCCCAGATTATTGCCGCTAAACGGGTGTACCA‑AATTAACC‑CCCCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc                                                                                  <  1:684868/142‑1 (MQ=255)
                                                                        ggggggggcccgggCGAAAATCGCTGCTGCCAAGATTATGCCCGCGGACCGGGGGAAGCA‑AATAAAAC‑CCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                                            <  1:940838/139‑1 (MQ=255)
                                                                          tggtggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCg                                                                          >  2:113402/1‑151 (MQ=255)
                                                                            gtggCCCCGGCCACAACCGCTGCTTGCCAGTTTTTTCCCGCT‑AATGTGCGTAACCAAATAAATC‑CCCCAATAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGt                                                                         <  1:145509/150‑1 (MQ=255)
                                                                                 tCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                   >  1:142653/1‑151 (MQ=255)
                                                                                  ccGGCGACAATCGCTGCTTGCCAGATTATGGCCGCTGAGCTGGGGTAGCA‑AATAAATC‑CCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                   <  2:833104/150‑1 (MQ=255)
                                                                                    gggaaaatccccggctTGCCAGATTATGGCCGGTGAAGGGGGGTAGCA‑AATAAACC‑CCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                   <  1:417483/135‑1 (MQ=255)
                                                                                          aaTCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGtt                                                           >  1:791029/1‑150 (MQ=255)
                                                                                          aaTCGCGGTTTCCCAGATTATTGCCGCTGAGCTGGCGTAGCC‑AATAAACC‑CCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGtt                                                           <  1:639818/150‑1 (MQ=255)
                                                                                                tgctTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc                                                    >  1:997481/1‑151 (MQ=255)
                                                                                                tgctTCCCAGATTTTTGCCGCTAAACGGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATaa                                                       <  2:199010/148‑1 (MQ=255)
                                                                                                  ctTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcac                                                  >  2:845607/1‑151 (MQ=255)
                                                                                                      cgagtttatccccGTTGAACGGGCGTGGCA‑AATAAATC‑CCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt                                              <  1:259963/140‑1 (MQ=255)
                                                                                                      ccAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt                                              >  2:249681/1‑151 (MQ=255)
                                                                                                      caaGATTATTGCCGTTAAACGGGCGTAGCA‑AATTAATC‑CCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                                               <  2:364176/148‑1 (MQ=255)
                                                                                                       cAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTg                                             >  1:819618/1‑151 (MQ=255)
                                                                                                        agtttattGCCGCTGAACTGGCGTAGCC‑AATTAATC‑CCGCAAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt                                              <  2:398088/146‑1 (MQ=255)
                                                                                                            tattGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATaa                                         >  2:517229/1‑150 (MQ=255)
                                                                                                            tattGCCGCTAAACTGGCGTAGCA‑AATTAATC‑CCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATaa                                         <  2:829340/150‑1 (MQ=255)
                                                                                                             attGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATaaa                                        >  1:294263/1‑150 (MQ=255)
                                                                                                                 ccGCTGAACGGGCGTAGCA‑AATTAATC‑CCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACtt                                     <  2:264078/149‑1 (MQ=255)
                                                                                                                    cTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt                                >  1:523427/1‑151 (MQ=255)
                                                                                                                    cTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt                                >  1:717489/1‑151 (MQ=255)
                                                                                                                     tGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGTATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGttt                               >  2:905606/1‑151 (MQ=255)
                                                                                                                       acggggggTAGCA‑AATAAATC‑CCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                              <  1:66953/143‑1 (MQ=255)
                                                                                                                       aaCGGGCGTGGCA‑AATAAACC‑CCGCAAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                              <  1:910464/150‑1 (MQ=255)
                                                                                                                              gTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCggagg                      >  2:977586/1‑151 (MQ=255)
                                                                                                                              gTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCggagg                      >  1:821855/1‑151 (MQ=255)
                                                                                                                                accc‑AATAAACC‑CCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGgaggag                    <  1:404432/149‑1 (MQ=255)
                                                                                                                                     aaataaacccccccAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTc                  <  1:462124/138‑1 (MQ=255)
                                                                                                                                     aaTAAATC‑CCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTc                  <  1:50314/149‑1 (MQ=255)
                                                                                                                                       taaaTC‑CCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTc                  <  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)
                                                                                                                                                       |                                                                                                                                                   
CCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCC‑AATTAATA‑CCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGC  >  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: ‑