Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,506,0280AC45.0% 83.2 / 18.0 61V18G (GTA→GGA) ydcOBenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base A (31/2);  new base C (0/27);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 4.61e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCA  >  NC_000913/1505880‑1506172
                                                                                                                                                       |                                                                                                                                                
tACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATAc                                                                                                                                                 <  1:22377/150‑1 (MQ=255)
      aTACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGa                                                                                                                                          >  1:722361/1‑151 (MQ=255)
        aCCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGaaa                                                                                                                                        >  1:408398/1‑151 (MQ=255)
          cATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAAc                                                                                                                                       >  2:1036919/1‑150 (MQ=255)
              aTGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTACTACCGCCAGAAACCCCg                                                                                                                                   >  2:940943/1‑150 (MQ=255)
              aTGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCg                                                                                                                                   >  1:1153784/1‑150 (MQ=255)
                        cAGCGTACTGACGACCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGCAGCgtg                                                                                                                        >  2:846890/1‑151 (MQ=255)
                          gCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCgtggg                                                                                                                      >  1:194671/1‑151 (MQ=255)
                             tACTGCCGCCAATTGCCACCCCCGGCGCCGCCACCCAGCAAGAGTTTTGTGCGGTGGTGGCCCCGGCAACAATCGCGGTTTGCAAGATTATGGCCGCGGAACGGGGGTGGCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCgtgggtgg                                                                                                                   <  2:1128226/151‑1 (MQ=255)
                                 gACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggagga                                                                                                               >  2:453806/1‑151 (MQ=255)
                                  ccccccATTCCCACCCCCAGCCCCTTCACCCAGCCAGAGTTTTGTGCGGTGGGGGCCCCGGCGACAACCGCTGCTTCCCAGATTTTGGCCGCGGAACGGGCGTAGCAAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTggagga                                                                                                               <  1:1132015/147‑1 (MQ=255)
                                       cATTGCGAGCCCACCCGCCTTCGTCGAGGCAAAGATTGGTGGAGTGGTGGCTCCGGCGACACCCGCGCGTTGCAAGATTCTTGCCGGTGAAGGGGCGAGGCAAATAAACCCCGCCGGAAACACCGCCAGTAGCGTGGGTGGAGGAATAGaa                                                                                                         <  2:565411/151‑1 (MQ=255)
                                        attcccccccccGGCGCCGTCCCCCGCCCAGAGTTTTGTGCGGGGGTGGCCCGGGCGAAAATCGCTGCTTGCCAGATTATTCCCGCTGAACTGGGGTAGCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGaaa                                                                                                        <  1:656194/143‑1 (MQ=255)
                                           gCCAGCCCCAGCGCGCTCGTCCAGCCAGAGTTTGGGGCTGGGTGGCCCCGGAGTAGAGTCGGCGCTGCCCAGTTTTCCGCCGCGGACGGGGCGGAACAAATAAATCCCACCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAAc                                                                                                      <  2:904784/150‑1 (MQ=255)
                                           gCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCCGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACa                                                                                                     >  2:965589/1‑151 (MQ=255)
                                               cccacggcgcCGACTCCCGGCCGGAGATTGGTGGGTGGGTGGGTCGGGCGACAACCGCTGCTGGCCAGATTAGTGCCGTTGTACGGGGGGAGAAAATAAAGCCCCCCAAAACCACCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACg                                                                                                 <  2:246345/145‑1 (MQ=255)
                                                ccccAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACg                                                                                                 >  2:947557/1‑150 (MQ=255)
                                                  ccggcgcCGCCACCCAGCCAGAGTTTTGTGCAGGGGTGGCCCGGGCGACAATCGCTGCTTGCCAGATTATGGCCGCTGAACGGGGGTAGCAAATTAATCCCCCCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCAt                                                                                              <  1:56417/148‑1 (MQ=255)
                                                      ccccGTCACCCACCAAGAGTTTTGGGCAGGGGGGGCTCCGGCAACAACCGCTGTTTCCAAAATTATGGCCGCTAAACGGGCGTAGCAAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                                                                                          <  1:941812/149‑1 (MQ=255)
                                                          gTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCAc                                                                                        >  1:968588/1‑149 (MQ=255)
                                                          ctCATCCGGCGAAGGATTTGGGGGGTGTTGGCCCCGGCCACATTGGCGTCTCCCCAGGTTTTTCCCGCGGAATTGGCGAACAAATTAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGAGGCATAGTCACCt                                                                                      <  1:848992/150‑1 (MQ=255)
                                                           ccatcccgcagaagaTTGGGGCAGTGGTGGCCCGGGAAACAATCCCTCTTTGCAGGATTATCGCCGTGAACGGGGCAAAAACAATAAACCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt                                                                                      <  2:867997/138‑1 (MQ=255)
                                                                 agccagAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGt                                                                               >  2:732642/1‑151 (MQ=255)
                                                                       agttttGGGCAGGGGGGGCTCCGGCAAAAATCGCTGCTTCCAAGATTATCGCCGATGAACTGGCGTAGCAAATTAACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAg                                                                          <  2:495163/147‑1 (MQ=255)
                                                                       agtttggtgCGGTGGGGGCCCCGGCAACAATCGCTGCTTGCCAGATTTTGCCCGCTAAAGGGGGGTAGCAAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAg                                                                          <  1:567427/144‑1 (MQ=255)
                                                                        gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                                        >  2:352678/1‑151 (MQ=255)
                                                                        gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                                        >  1:405440/1‑151 (MQ=255)
                                                                         ttttGTGCGGTGGGGGCTCCGGCGACATTCGCGGCTTGCCAGTTTATGGCCGCTGACCGGGTGTACCAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGt                                                                       <  2:968178/150‑1 (MQ=255)
                                                                          tttGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGCAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCCTAGTCACCTCGTGCGTTATAGCGTa                                                                      >  2:867109/1‑151 (MQ=255)
                                                                           tgttgCGGTGGTGGCCCCGGCGAAAACCGCTCCTTGCCAATTTTTGGCCGCTGAACTGGCGTAGCAAATTAAACCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc                                                                     <  1:804452/148‑1 (MQ=255)
                                                                                cggggggggcccgggCGAAAACCGCTGTTTGCAAAATTTTGGCCGCTAAACGGGGGTAGCAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCa                                                                <  1:690832/138‑1 (MQ=255)
                                                                                  gtggtgGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                               >  2:692854/1‑150 (MQ=255)
                                                                                  gtggtgGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                               >  2:586015/1‑150 (MQ=255)
                                                                                  gtggtgGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                               >  1:433120/1‑150 (MQ=255)
                                                                                  gtggtgGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                               >  1:394959/1‑150 (MQ=255)
                                                                                  gaggtggtGCCGGCCACAGTCGCTGCTTGCCATTTTTGTGCGGGTGACGTGGGGGACAAAATTAATCCCGCCAAAAACACCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                               <  2:436768/148‑1 (MQ=255)
                                                                                         ccccGGCGACAATCGCTGCTGGCAAGATTATTGCCGCTGACCTGGCGTACCAAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAt                                                       <  1:526705/149‑1 (MQ=255)
                                                                                          ccgggCAACACTCGTTGTTGGCAAGTTTATTGCCGCTGAACGGGCGTAGCCAATTAATCCCGGCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                      <  1:979897/148‑1 (MQ=255)
                                                                                            cGGCGACAACCCCTGCTTGCCAGTTTTTTCCCGCTGAACTGGCGTAGCAAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTg                                                    <  1:732979/151‑1 (MQ=255)
                                                                                               cGACAGTCGCTGCTCCCAATATTGTCGCCGGTGACGGGGCGTAGCAAATAAACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTc                                                  <  2:595639/150‑1 (MQ=255)
                                                                                                gACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCg                                                >  1:186759/1‑151 (MQ=255)
                                                                                                     accgctctttGCCAAATTTTGGCCGCGAACCGGGGGTAGCCAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTtat                                            <  1:195091/143‑1 (MQ=255)
                                                                                                           gctTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGca                                      >  1:363888/1‑151 (MQ=255)
                                                                                                            tttGCCAGATTATGCCCGATAGAGTGGCGGAGCCAATAAACGCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcac                                     <  2:1059730/150‑1 (MQ=255)
                                                                                                               cccAGTTTTTTGCCGCTAAACGGGCGTAGCAAATTAATCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                                  <  1:918435/150‑1 (MQ=255)
                                                                                                               cccAGATTATGGCCGTGGAACGGGGTTAGCAAATAAATCCCGCCGGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                                  <  2:81434/150‑1 (MQ=255)
                                                                                                               cccAGATTATGCCCGCTAAACTGGCGAGGCCAATTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                                  <  1:501809/150‑1 (MQ=255)
                                                                                                                   gATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGCGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGAt                              >  2:327328/1‑151 (MQ=255)
                                                                                                                   gATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGAt                              >  2:77317/1‑151 (MQ=255)
                                                                                                                   gATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGAt                              >  1:568786/1‑151 (MQ=255)
                                                                                                                       aagagagaggAACTGGCATAGCAAATAAATACCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAAc                          <  2:38874/142‑1 (MQ=255)
                                                                                                                              ctaaaCGGGGGTAGCAAATAAACCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGt                    <  1:1191127/147‑1 (MQ=255)
                                                                                                                              cggAACGGGGGTACCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGt                    <  1:51206/148‑1 (MQ=255)
                                                                                                                              cTGACCTGGCGTAGCCATTTATTCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt                   <  2:236640/151‑1 (MQ=255)
                                                                                                                                gAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAAc                          >  1:38881/1‑142 (MQ=255)
                                                                                                                                gAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                 >  2:7017/1‑151 (MQ=255)
                                                                                                                                      ggCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgg            >  1:1027222/1‑151 (MQ=255)
                                                                                                                                       gCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCATAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgga           >  1:737579/1‑151 (MQ=255)
                                                                                                                                            gCCAATTAATACTGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGt      >  2:334393/1‑151 (MQ=255)
                                                                                                                                                aTTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCa  >  1:945537/1‑151 (MQ=255)
                                                                                                                                                aTTAATACCGCCAGAAAAACCGCCAGCAGCGTGGGTGGAGGAATAGAAAACAGACGCATACTCCCCTCGTGCGTTACAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACTCACAAGTGATAAACTTCCGTTTTGCCGGAGGCGACGc   >  2:260306/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
TACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCA  >  NC_000913/1505880‑1506172

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