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

CATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGCTCG  >  NC_000913/1505890‑1506188
                                                                                                                                                      |                                                                                                                                                       
cATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑cc                                                                                                                                                      >  2:913954/1‑151 (MQ=255)
   aaTGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCg                                                                                                                                                   <  2:1214689/151‑1 (MQ=255)
    aTGTCAGAGTCAGCGTACTGACGCCCATTGTCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGc                                                                                                                                                  >  2:96943/1‑151 (MQ=255)
    aTGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGc                                                                                                                                                  >  1:1431856/1‑151 (MQ=255)
      gTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCa                                                                                                                                                <  1:1494595/151‑1 (MQ=255)
       tCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGcc                                                                                                                                                 >  1:702274/1‑149 (MQ=255)
       tCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCa                                                                                                                                                <  2:1055887/150‑1 (MQ=255)
        cAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGt                                                                                                                                              >  2:216329/1‑151 (MQ=255)
            gTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCg                                                                                                                                          >  2:663063/1‑151 (MQ=255)
                 cGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCgtggg                                                                                                                                      >  2:565719/1‑150 (MQ=255)
                     cTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCgtgggtgg                                                                                                                                   >  2:1394585/1‑149 (MQ=255)
                      tGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTggagg                                                                                                                                >  2:695648/1‑151 (MQ=255)
                      tGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTggagg                                                                                                                                >  2:1056203/1‑151 (MQ=255)
                       gACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTggagg                                                                                                                                >  2:858979/1‑150 (MQ=255)
                        aCGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTggagga                                                                                                                               >  1:1317345/1‑150 (MQ=255)
                           cccttggcacgCCCAAGCGCCGTCATCCGCCCAGAGTTTTGGGCGGGGGGGGCCCGGGCAAAATCCGCTCTTTGCCAGATTTTGCCCGCGGACTTGGGGT‑AGCCATTT‑AATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTggagg                                                                                                                                <  2:407008/136‑1 (MQ=255)
                           ccctgccccccccccGGCCCCGTCACCCAGCGAGAGTTTTGGGGGGGGGGGCCCCCGGCGACAACCGCTCCTTGCCAGATTTTGGCGGCTGAGGGGGGGG‑AGCAAATT‑AACCCCGCAAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAAt                                                                                                                             <  1:595309/138‑1 (MQ=255)
                           cccATTGCCACCCCCGGCCCCGTCACCCAGCCAGAGTTTTGTGCGGGGGGGGCCCGGGCGAAAATCGCTTCTTCCCAGTTTATTCCCGCTGACCGGGCGT‑AGCAAATA‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATa                                                                                                                            <  2:986093/150‑1 (MQ=255)
                               ttGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGaaa                                                                                                                        >  1:500644/1‑150 (MQ=255)
                                  ccAGCCCCGGCGCCGTCACCCAGCCAGAGATTGGTGCAGGGGTGGCTCCGGCGACAACCGCGGCTTGCCAGTTTATTCCCGCTGAACGGGCGT‑AGCCAATT‑AATCCCCCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAg                                                                                                                    <  2:1123797/151‑1 (MQ=255)
                                    aGCCCCAGCGCCGTCATCCAGCAAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGa                                                                                                                   >  1:979823/1‑150 (MQ=255)
                                     gCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACg                                                                                                                 >  2:691192/1‑151 (MQ=255)
                                       cccGGCCCCGTCACCCACCAAGAGTTTTGTGCGGGGGGGGCCCCGGCGAAAATCGCTGTTTGCAAGATTTTTCCCGTTGAACGGGCGT‑GGCAAATA‑AATCCCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCa                                                                                                               <  2:965365/151‑1 (MQ=255)
                                              ccctCACCCAGCCGGAATTTTGCGCGGGGGTGGCCCGGGCGACATTCGCTGCTTCCCAGTGTATTCCCGCTGACCGGGCGT‑GGCAAATA‑AATACCGCCAAAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCAc                                                                                                        <  2:449073/148‑1 (MQ=255)
                                                  cccccagccagAGTTTTGGGCGGGGGTGGCCCCGGCGACATCCGCTCCTTGCCAGTTTATGGCGGTTGAAGGGGCGT‑AGCAAATA‑AATCCCGCAAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAg                                                                                                            <  1:9946/140‑1 (MQ=255)
                                                   acccagcgagAGTTTGGTGCGGTGGGGGCCGGGGCAACATCCGCTGTTGCCCAGTTTATGGCCGTTGAGGGGGGGG‑GGCAAATA‑AACCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCg                                                                                                    <  1:545964/148‑1 (MQ=255)
                                                         caggagtTTTGTGCGGTGGTGCCTCCGGCAAAAACCGCTGCTTCCAAGATTTTTGCCGGAGACGTGGGGA‑AGCAAATA‑AACCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGt                                                                                               <  1:74117/146‑1 (MQ=255)
                                                         caagggaTTTGGGAGGTGGTGGCTCCGGCGAAAACCGCTGCTGGCAAGATTTTTGCCGGGGAGTGGGGGT‑AGCAAATA‑AACCCCGCAAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt                                                                                                      <  1:349479/137‑1 (MQ=255)
                                                         caagatttttGTGGGGTGGGGGCCCCGGCAACAACCGCTGTTTGCAATTTTTTTGCCGTTAAGGGGGCGT‑AGCAAATA‑AATCCCCCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTc                                                                                                     <  1:1171961/136‑1 (MQ=255)
                                                         caagagATTTGTGCAGTGGGGGTCCGGGCGAAAACCGCTGTTTGCCAGATTATGGCCGCGAACCGGGCGT‑AGCAAATA‑AATCCCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt                                                                                              <  2:92033/148‑1 (MQ=255)
                                                           aggtgtttGTGCGGTGGGGGCCCCGGCGAAAACCGCTGCTTGCAATATTATGGCGGCTGACCGGGCGT‑GGCAAATA‑AACCCCGCAAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat                                                                                            <  1:1521299/145‑1 (MQ=255)
                                                           agagTTTTGTGCAGGGGTGGCCCGGGCGACAACCGCTGCTGGCCAGATTATGGCCGCTAAACTGGCGT‑AGCAAATA‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTta                                                                                             <  2:1494737/149‑1 (MQ=255)
                                                                     gcgggggtggCCCGGGCGACAATCGCTGCTTGCCAGTTTTTGGCCGCTGAACTGGCTT‑AGCCAATA‑AATCCCGCCAGAAC‑CGCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGc                                                                                 <  2:1436893/146‑1 (MQ=255)
                                                                       aGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                                               >  2:657025/1‑151 (MQ=255)
                                                                       aGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                                               >  2:729647/1‑151 (MQ=255)
                                                                          gggggCTCCGGCGACAATCGCTGCTTCCAAAATTATTGCCGCTGAAGGGGCGT‑ACCAAATT‑AATCCCCCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCa                                                                                <  1:883587/144‑1 (MQ=255)
                                                                            tggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGt                                                                            >  1:734836/1‑149 (MQ=255)
                                                                                tCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                      >  1:1344179/1‑151 (MQ=255)
                                                                                 ccGGCAAAAATCGCTGTTTCCAAAATTATTGCCGCTGAACGGGCGT‑A‑CCAATTAAACCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCa                                                                        <  2:1456039/148‑1 (MQ=255)
                                                                                   ggCGACAATCGCGGTTTCCAAGATTATCGCCGGGGAAGGGGGGA‑GGCAAATA‑AAACCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGc                                                                         <  1:376922/145‑1 (MQ=255)
                                                                                   ggCGACAACCGCTGCTTGCAATTTTTTTGCCGCTAACCGGGCGT‑AGCAAATA‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCa                                                                        <  2:1442596/146‑1 (MQ=255)
                                                                                      caaaaccGCTGCTTGCCAGATTATCGCCGCGAAACGGTCGA‑AGCCAAAA‑AAACCCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCg                                                                <  2:154502/145‑1 (MQ=255)
                                                                                       aCAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGt                                                               >  1:335858/1‑151 (MQ=255)
                                                                                       aCAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGt                                                               >  1:80863/1‑151 (MQ=255)
                                                                                         aaTCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGtt                                                              >  1:973841/1‑150 (MQ=255)
                                                                                           tCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTtat                                                            >  2:980624/1‑150 (MQ=255)
                                                                                             gctgctTGCCAGATTATTGCCGCTGAACGGGCGG‑AGCAAATT‑AATCCCCCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATaa                                                          <  2:683009/150‑1 (MQ=255)
                                                                                               tttttGCCAGTTTATTGCCGCTGAACGGGCGT‑ACCAAATA‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATaa                                                          <  2:261375/145‑1 (MQ=255)
                                                                                                  ttGCCAGATTATCCCCGGTGAATGGGCGT‑ACCAAATA‑AACCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATaa                                                          <  1:316503/145‑1 (MQ=255)
                                                                                                     caaattttttcccGCGGAACGGGCGT‑AGCAAATT‑AACCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt                                                 <  1:499626/140‑1 (MQ=255)
                                                                                                       atattattGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAATCAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGGCCGTTATAACGCACAAGTg                                                >  1:1163592/3‑150 (MQ=255)
                                                                                                       agtttttggCCGCTGACCGGGCGT‑AGCAAATT‑AATACCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                                                  <  1:522838/140‑1 (MQ=255)
                                                                                                       aattttttGCCGC‑GGACTGGGGTAAGCCAATA‑AATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc                                                       <  1:1347264/137‑1 (MQ=255)
                                                                                                          ttttggCCGTTGAGCGGGCGT‑GGCAAATT‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGa                                               <  1:398960/143‑1 (MQ=255)
                                                                                                          ttattGCCGCTGAACGGGCGT‑AGCAAATA‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGAt                                              <  2:1046077/149‑1 (MQ=255)
                                                                                                             ttGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACt                                         >  1:491979/1‑151 (MQ=255)
                                                                                                              tGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACtt                                        >  1:716429/1‑151 (MQ=255)
                                                                                                                ccGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTcc                                      >  2:868760/1‑151 (MQ=255)
                                                                                                                ccGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTcc                                      >  1:172493/1‑151 (MQ=255)
                                                                                                                  gctaatgtGGGGT‑CGCAAATA‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAAc                                          <  1:256645/138‑1 (MQ=255)
                                                                                                                      gagctggggaaccaaaTT‑AA‑ACCCCCCGAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                                 <  1:784525/137‑1 (MQ=255)
                                                                                                                       aCTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc                               >  2:915353/1‑151 (MQ=255)
                                                                                                                          ggCGT‑GGCAAATA‑AATCCCGCAAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGt                                    <  1:1229930/143‑1 (MQ=255)
                                                                                                                           ggggtaccaaaTA‑AATCCCCCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCg                             <  2:80871/141‑1 (MQ=255)
                                                                                                                           gggggaGCAAATA‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgg                            <  1:274547/146‑1 (MQ=255)
                                                                                                                                  caaaTA‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCggagg                         <  2:28544/145‑1 (MQ=255)
                                                                                                                                   cAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGCGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCg                    >  2:1045824/1‑151 (MQ=255)
                                                                                                                                   cAATT‑AATACAGCCAGAAA‑ACCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTc                     >  1:387214/1‑150 (MQ=255)
                                                                                                                                       t‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATg                >  1:74545/1‑151 (MQ=255)
                                                                                                                                          aTACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATgg               >  1:767850/1‑150 (MQ=255)
                                                                                                                                                  aGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATctcg      <  2:909364/151‑1 (MQ=255)
                                                                                                                                                     aa‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATctcgctc   >  2:35250/1‑151 (MQ=255)
                                                                                                                                                      a‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATctcgctcg  >  2:1098221/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
CATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGT‑AGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGCTCG  >  NC_000913/1505890‑1506188

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