Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA48.3% 39.4 / 34.9 58L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (22/8);  new base A (0/28);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 1.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA  >  NC_000913/357740‑358034
                                                                                                                                                      |                                                                                                                                                    
cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                      >  1:106046/1‑151 (MQ=255)
  gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                                    <  2:233440/151‑1 (MQ=255)
   aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                   <  1:935417/151‑1 (MQ=255)
     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                                  <  1:824791/150‑1 (MQ=255)
     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                >  2:825818/1‑151 (MQ=255)
                  cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                   >  2:678474/1‑151 (MQ=255)
                       aaTTCCAGCCGCAGACCTGTGTCGGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTCCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                               <  2:330478/150‑1 (MQ=18)
                       aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                >  1:685403/1‑149 (MQ=255)
                         ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                            >  1:782403/1‑151 (MQ=255)
                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                         >  1:264038/1‑151 (MQ=255)
                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                         >  2:738789/1‑151 (MQ=255)
                             aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                        >  1:438529/1‑151 (MQ=255)
                               ccGCAGCCCTGTGTCAGCGGCTACCGTGATTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                           <  2:812628/146‑1 (MQ=18)
                                cGCAGCCTTTTGCCAGCGTCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTATCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                                     <  2:359209/151‑1 (MQ=255)
                                           gccAGCGGCTCCCGTGATTCTTTCCCCCCAACACCCGCGCATTTACCCACCGCCATTTGAAAAAATCCCTTCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                                                                                         <  1:1006419/148‑1 (MQ=255)
                                              aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaa                                                                                                           >  1:832733/1‑145 (MQ=255)
                                              aGCGGCTACCGTGATTCATTCCCCCCACAACCCGCGCATTTACCCACCGCCATGTGCAAAAATCCTTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                                                                                                       <  2:934878/149‑1 (MQ=25)
                                                 gccTACCGGTATTCTTTCCCCCCACCACCCGCCCTTTTACCCAACCCTTTTTCAAAAAATGCCTCCGGGGCGGCTGTCGCCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                                  <  1:529722/149‑1 (MQ=255)
                                                   cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAtgaacaaaaaaagttaa                                                                                                   >  1:618991/1‑136 (MQ=25)
                                                      ccGTGATTCTTTCCCGCCACCACCCGCGATTTTACCCACCGCCATGTGAAAAAATCCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                                                 <  2:601661/147‑1 (MQ=255)
                                                           aTTCTTTCCCGCAAACACCCGCTCTTTTACCCACCCCCGTGTAAAAAATTCCCTCCGCGGCGGTTTTCTCCCAGTTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt                                                                                         <  2:523610/150‑1 (MQ=255)
                                                                ttcccGCCAACACCCCCGCATTTCCCCACCGCCATGTCAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTCCCCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt                                                                                      <  2:695945/145‑1 (MQ=255)
                                                                  cccGCCAACACCCGCGCTTTTCCCCAACGCCATGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                                 <  1:272990/150‑1 (MQ=255)
                                                                   ccccaaccacccGCCCATTAACCCACCCCCATGTAAAAAAATGCCTCCGGGGCGGTTGTCGGCCACTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg                                                                                   <  1:187623/138‑1 (MQ=255)
                                                                    cGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAACCCCCCCTGTTGTGTgg                                                                                >  2:84432/1‑150 (MQ=25)
                                                                         acaccCGCGCTTTTACCCAACGCCATGTGAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                          <  2:398284/151‑1 (MQ=255)
                                                                            cccGCGCTTTTATCCACCGCCATGAGAAAAACTGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                         <  1:696502/148‑1 (MQ=255)
                                                                                cgcATTTACCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                                    <  2:763253/150‑1 (MQ=255)
                                                                                     taccccacctccATGTGAAAAAAGCCCTTCGCGGCGGCTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                               <  1:696805/140‑1 (MQ=255)
                                                                                      tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGAAAAAAATTTTAACCCCCCATGTTGTGTGGGCATCGCTGCTggggaaa                                                              >  2:40661/1‑145 (MQ=25)
                                                                                         ccAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAAACCCCCATGTTGTGTGGGGAACGCTGCTGGGAACaa                                                            >  2:1005114/1‑148 (MQ=25)
                                                                                         ccAACCCCATTTGCAAAATTCCTTTCGCGGGGGCTGTCTGCCACCTATTTTTCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                                 <  2:984284/144‑1 (MQ=255)
                                                                                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAGCTTTTAACCCCCCCTGTTTTGTGGGGATCGCTGCTGGGTACAAagtg                                                        >  2:1012915/1‑149 (MQ=37)
                                                                                                cATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGAAAAAACTTTAAATCCGTCCTTTTGTGTGGCAAAAGCTGCTGGGTACaa                                                            >  1:763907/1‑141 (MQ=37)
                                                                                                cATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAACCCCCCCTGTTTTTTGGGGATCGCTGCTGGGAaaaaaatgggaa                                                    >  2:821808/1‑141 (MQ=25)
                                                                                                   gtgAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                                 <  2:827136/149‑1 (MQ=255)
                                                                                                       aaaaaTGCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                             <  2:676421/149‑1 (MQ=255)
                                                                                                       aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                                <  2:522369/146‑1 (MQ=255)
                                                                                                       aaaaaTGCCTCCGCGGCGTCTTCCTCCCAGTTTTTTTCCCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                                <  1:514051/146‑1 (MQ=255)
                                                                                                       aaaaaTGCCTCCGCGGCGGCGGTCTGCCGGCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                             <  2:782480/149‑1 (MQ=255)
                                                                                                        aaatTGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                          <  2:125560/151‑1 (MQ=255)
                                                                                                          aaTGCCTCCGCGGCGGCTGTCTCCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                        <  2:629024/151‑1 (MQ=255)
                                                                                                              ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGAGGAAAAACCTTTAAACCCGCCTGTTGTGGGGGGAACGCTGCTGGTGACAATGAGGGAAggggggggggggg                                       >  2:246475/1‑138 (MQ=255)
                                                                                                                ttCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAACTTTAACCCCCCATGGTGTGTGGGGAACGCTGCTTGTGAAAAAATGGGAAAGGaaaggggggg                                       >  2:378283/1‑138 (MQ=255)
                                                                                                                   gcggcgGCTGTCTGCGAGCTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg                                 <  1:104186/149‑1 (MQ=255)
                                                                                                                        ggttGTCTGCCAGTTTTTTTCCCCCCGAAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                            <  1:878376/146‑1 (MQ=255)
                                                                                                                           tGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACTTTTAACCCGCCCTGTTTTGTGGGGATCGCTGCTGGGGAAAATGTGGGAAGGGGAGGggggcggaaaagcacc                             >  2:274721/1‑135 (MQ=255)
                                                                                                                                 ccAGTTTTTTTCCCCCCCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg                    <  1:151999/148‑1 (MQ=255)
                                                                                                                                 ccAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGGGAAATAAATTTAAACCCGTCCTGTTTTGTGGGGAACGCTGCTGGTGACAATATGGGAAGGGGAGGGGGGCGGATCAGCTCCAGAACcgc                     >  1:730047/1‑147 (MQ=255)
                                                                                                                                    gCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGCCAGAGAAAAAAACTTTAACCCCCCACGTTTTGTGGGGATTGCTGCTGGGGAAAAAGTGGGAAGGGGAGGGGGGGGGATCACCCCCAGAACCGCGCTa                 >  1:298083/1‑148 (MQ=255)
                                                                                                                                     cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAACGAAAAAACTTTTAAACCGCCCTCTTTTTTGGCAAGCGCTGCTCGTTACAATGGGGGAAGGGGAGGGCGGGGGATCAACTCCCGAACCCCGCCAAtt              >  2:464315/1‑150 (MQ=255)
                                                                                                                                      ttttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt               <  1:399471/146‑1 (MQ=255)
                                                                                                                                        ttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg           <  1:711346/150‑1 (MQ=255)
                                                                                                                                        ttttcCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt              <  2:741399/147‑1 (MQ=255)
                                                                                                                                          ttcccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt         <  2:87382/147‑1 (MQ=255)
                                                                                                                                             ccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc        <  2:513697/148‑1 (MQ=255)
                                                                                                                                             cccccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc     <  2:792076/148‑1 (MQ=255)
                                                                                                                                              cgccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGAAAAAAAATTTAATCCCCCCTGTTTTGTGGGCATAGCTTCTGGGAACAAAATTGGAAAGGGAGGGGGGGGGAACAACTCCACAACCGCGCCAATTTAGGTCGcc     >  2:747329/1‑150 (MQ=255)
                                                                                                                                                ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa  <  2:26159/151‑1 (MQ=255)
                                                                                                                                                ccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc   <  1:1005014/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA  >  NC_000913/357740‑358034

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