Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA45.0% 57.8 / 33.3 60L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (22/11);  new base A (0/27);  total (22/38)
Fisher's exact test for biased strand distribution p-value = 1.44e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357745‑358038
                                                                                                                                                   |                                                                                                                                                        
aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                    <  1:192832/151‑1 (MQ=255)
aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                    >  2:803513/1‑151 (MQ=255)
 aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                   >  2:1002955/1‑151 (MQ=255)
   cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                                  <  2:1088880/150‑1 (MQ=255)
    gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                                >  1:907089/1‑151 (MQ=255)
    gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                                <  2:905316/151‑1 (MQ=255)
    gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                                 <  1:517911/150‑1 (MQ=255)
           aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                                                                                                                                         >  1:1042438/1‑151 (MQ=255)
              aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                       >  2:1115605/1‑150 (MQ=255)
                gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                    >  1:492135/1‑151 (MQ=255)
                          ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                          >  2:1120553/1‑151 (MQ=255)
                          ccGAAGCCCTGTGTCAGGGCCTCCCGTGATCCTTTCCCGCCACCACCCGCGCTTTTACCCAACGCCTTGTGCAAAAATGCTTTCGCGGCGCCTGTCTGCCAGTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                             <  1:346115/148‑1 (MQ=255)
                             cagccttgtgTCGGCGGTTACCGTGATTCATTCCCCCCAACCCCCGCCCTTTTACCCACCGCGTTTTAAAAAATTCCCTGCGCGGGGGCTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                                        <  1:780243/144‑1 (MQ=255)
                               gccttttgtCAGGGGCTCCCGTTTTCCTTTCCCCCCAACCCCCGCGCTTTTACCCACCGCATTGTGAAAAATTCCCTCCGGGGGGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                                       <  2:189059/143‑1 (MQ=255)
                                cccTGTGTCGGCGCCTCCCGTGATTCATTCCCCCCAACACCCGCCAATTTACCAACCCCCATGTGAAAAAATGCTTTGGGGGCGGCTTTCTGCCGGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                   <  1:970074/150‑1 (MQ=255)
                                      gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGaaa                                                                                                               >  1:1163795/1‑148 (MQ=255)
                                       tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcaaaggaaaa                                                                                                               >  2:70114/1‑138 (MQ=255)
                                        cagcgccccccGTGTTTCTTCCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAAGCCCTTCGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                                       <  1:860479/131‑1 (MQ=255)
                                            gcctcccGTGATTCATTTCCCCCAACAACCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                                      <  2:343693/146‑1 (MQ=255)
                                               tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCCTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGaaaaaaattttaaac                                                                                                   >  2:756624/1‑139 (MQ=25)
                                                cccgtttttCATTCCCCCCAACACCCGCGCATTTACCCACCCCCTTGTGAAAAAATGCCTCCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                                                                                                           <  1:634222/135‑1 (MQ=255)
                                                 ccGTGTTTCATTCCCCCAAACAACCGCGCATTTACCCAACGCCTTGTAAAAAACGCCGTTCGCGGGGGGTCTCTGCCGCCTATTTTCCCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                                   <  2:486866/149‑1 (MQ=255)
                                                 ccGTGATTCTTTCCCCCCAACACCCCCGCTTTTACCCACCCCCATGTGAAAAAATCCCTTCGGGGGGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                                   <  2:1147132/149‑1 (MQ=255)
                                                  cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAAccc                                                                                                  >  2:1011073/1‑149 (MQ=25)
                                                    taaTTCTTTCCCCCCAACACCCGCGCTTTTATCCAACGCCATGTAAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                                <  1:66639/147‑1 (MQ=255)
                                                       ttCATTCCCGCCAACACCCGCGTTTTTACCAACCCCAATGTGAAAAAAGCCTTTCGCGGCGTCTTTGTGCAAGCTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                                            <  2:1145357/150‑1 (MQ=255)
                                                         cATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAGAATTAAGacccccagttcg                                                                                         >  1:155178/1‑142 (MQ=25)
                                                         cATTCCCCCCAACACCCGCGCATTTACCAACCGCCATGTGCAAAAATCCTTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                                           <  2:895164/149‑1 (MQ=255)
                                                           ttcccccaaacaccCGCGCTTTTACCCAACGCCATGTGCAAAAATCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt                                                                                        <  2:283485/142‑1 (MQ=255)
                                                           ttcccGCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                                         <  2:1194130/146‑1 (MQ=255)
                                                             cccGCCAACACCCGCGTATTTACCCACCGCCTTGTAAAAAATTGCTTCCGGGGGGGTGGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                                         <  1:1120003/146‑1 (MQ=255)
                                                              ccccaaacaccCGCGCTTTTACCCAACGCCATGTAAAAAAATGCTTGCGCGGCGGCTGTCTGCAAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                                    <  1:931570/145‑1 (MQ=255)
                                                               cGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAAATTTTAACCCCCCCTGTTTTGTggg                                                                                  >  2:1127865/1‑150 (MQ=25)
                                                               cGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAGCTTTAAACCCGCCCTGTTTTGTgg                                                                                   >  1:200887/1‑150 (MQ=255)
                                                                 ccAACACCCGCGCATTTCCCCCACGCCATTGGAAAAAATGCCTTCGCGGCGGCTTTTTGCCCGTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                                 <  2:1009480/150‑1 (MQ=255)
                                                                      aaccgcaattttACCAAACCCCATGTCAAAAAATCCTTGCGGGGCGGTTTTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                            <  1:836903/141‑1 (MQ=255)
                                                                         cgcttttttACCCAACCCCATGTACAAAAATCCCTTCGCGGCGCCTTTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                              <  1:879879/139‑1 (MQ=255)
                                                                            gtaTTTACCCAACGCCATGTGAAAAAGTGCCTGCGCGGCGGCTGTCTGCAAATTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                                     <  2:907730/149‑1 (MQ=255)
                                                                                      aaCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAAATTTTAACCCCGCATGTTGTTTGGGGATCGCTGCgggggaaaaagtgg                                                           >  2:349295/1‑140 (MQ=25)
                                                                                       aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGAAAAAAAATTTATACCCCTCCTTTTTTGGGGGGATCGCTGCTGGGAaaaaaatg                                                            >  1:379811/1‑145 (MQ=19)
                                                                                        cGCCATTTCAAAAATTGCTTCCGCGGGGGCTTCCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                                         <  2:1168304/151‑1 (MQ=255)
                                                                                             aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAAATTAAACCCCTCCTGTTTTGTGGGAATCGCCGCTGGtgaaaaaatgg                                                           >  2:1095333/1‑135 (MQ=25)
                                                                                              tgtaCAAAAATGCTTTGGCGCCGGTTGTCTCCCTTCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                                      <  2:183160/149‑1 (MQ=255)
                                                                                               gtgAAAAAATGCTTGCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                                     <  2:719930/149‑1 (MQ=255)
                                                                                                   aaaatTGCTTGCGCGCCGGTTTTCTGCCATTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                               <  2:699656/151‑1 (MQ=255)
                                                                                                   aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                                 <  2:445070/149‑1 (MQ=255)
                                                                                                     aaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAGTTTTAATCCCGCCTTTTGTGGGGGCAATGCCGCTGGTGACAAAGTGGGAAGGGGAgggggg                                             >  2:245170/1‑146 (MQ=255)
                                                                                                      aaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAAAAAAAAATTTTAACCCCCCCTGTTGTGTGGGGAACGCTGCTGGGGACAAAGTGGGAAGGGGAgggggg                                             >  2:289343/1‑145 (MQ=255)
                                                                                                      aaTCCCTCCGGGGGGGTTGTCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                            <  1:246686/151‑1 (MQ=255)
                                                                                                           ctttCGGGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                            <  2:1066737/145‑1 (MQ=255)
                                                                                                                  ggcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGCGCGAGAGAAAAAAGTTTTAAACCCGCCATTTGTTTTGGCAACTCCGCCTGGTAAAAAATGGGAAGGGGAAGGGGGGGGAAtaactcc                                  >  2:421487/1‑145 (MQ=255)
                                                                                                                     gccTGTCCCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGa                                          <  1:596731/137‑1 (MQ=255)
                                                                                                                               ccattATTTTCCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc                       <  1:585288/145‑1 (MQ=255)
                                                                                                                                   ttttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt                  <  2:609120/147‑1 (MQ=255)
                                                                                                                                   tATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGCGAGGGAAAAAAATTTTAAACCCTCCTTTTTTTTGGGGATCGCTGCTGGGGACAATATGGGAAGGGGAGGGGGGGGGATCAACCCCCCAACCGCGCCAAtt                  >  1:429274/1‑149 (MQ=255)
                                                                                                                                     ttttcCCCCCGAAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg               <  1:1267810/150‑1 (MQ=255)
                                                                                                                                       tttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCGAGAAAAAAATTTAATCCCCCCATGTTGTGTGGGGAACGCTGCTGGGGAAAATGTGGGGAAGGGAGGGGGGGGGAGCAGCTCCACAAcccccccca                    >  2:994692/1‑136 (MQ=255)
                                                                                                                                        tcccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc            <  1:548377/148‑1 (MQ=255)
                                                                                                                                         cccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc            <  2:733242/148‑1 (MQ=255)
                                                                                                                                          ccgccgCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc         <  1:1196774/151‑1 (MQ=255)
                                                                                                                                                 aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAAt    <  2:203455/149‑1 (MQ=255)
                                                                                                                                                 aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa  <  2:182838/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357745‑358038

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