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→TTT) cynRtranscriptional 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

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA  >  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)
     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑c                                                                                                                                                >  2:825818/1‑151 (MQ=255)
                  cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTcc                                                                                                                                   >  2:678474/1‑151 (MQ=255)
                       aaTTCCAGCCGCAGACCTGTGTCGGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTCCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCagta                                                                                                                               <  2:330478/150‑1 (MQ=18)
                       aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt                                                                                                                                >  1:685403/1‑149 (MQ=255)
                         ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTg                                                                                                                            >  1:782403/1‑151 (MQ=255)
                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                         >  1:264038/1‑151 (MQ=255)
                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                         >  2:738789/1‑151 (MQ=255)
                             aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                        >  1:438529/1‑151 (MQ=255)
                               ccGCAGCCCTGTGTCAGCGGCTACCGTGATTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTgg                                                                                                                           <  2:812628/146‑1 (MQ=18)
                                cGCAGCCTTTTGCCAGCGTCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTATCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                                     <  2:359209/151‑1 (MQ=255)
                                           gccAGCGGCTCCCGTGATTCTTTCCCCCCAACACCCGCGCATTTACCCACCGCCATTTGAAAAAATCCCTTCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATa                                                                                                         <  1:1006419/148‑1 (MQ=255)
                                              aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑aaa                                                                                                           >  1:832733/1‑145 (MQ=255)
                                              aGCGGCTACCGTGATTCATTCCCCCCACAACCCGCGCATTTACCCACCGCCATGTGCAAAAATCCTTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGc                                                                                                       <  2:934878/149‑1 (MQ=25)
                                                 gccTACCGGTATTCTTTCCCCCCACCACCCGCCCTTTTACCCAACCCTTTTTCAAAAAATGCCTCCGGGGCGGCTGTCGCCCAGTTATTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTaa                                                                                                  <  1:529722/149‑1 (MQ=255)
                                                   cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAtgaacaaaaaaagttaa                                                                                                   >  1:618991/1‑136 (MQ=25)
                                                      ccGTGATTCTTTCCCGCCACCACCCGCGATTTTACCCACCGCCATGTGAAAAAATCCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAg                                                                                                 <  2:601661/147‑1 (MQ=255)
                                                           aTTCTTTCCCGCAAACACCCGCTCTTTTACCCACCCCCGTGTAAAAAATTCCCTCCGCGGCGGTTTTCTCCCAGTTATTTTTCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAt                                                                                         <  2:523610/150‑1 (MQ=255)
                                                                ttcccGCCAACACCCCCGCATTTCCCCACCGCCATGTCAAAAAATCCCTCCGCGGCGGCTGTCTGCCAG‑TTTTTTCCCCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgtt                                                                                      <  2:695945/145‑1 (MQ=255)
                                                                  cccGCCAACACCCGCGCTTTTCCCCAACGCCATGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                                 <  1:272990/150‑1 (MQ=255)
                                                                   ccccaaccacccGCCCATTAACCCACCCCCATGTAAAAAAATGCCTCCGGGGCGGTTGTCGGCCACTTATTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTtgtg                                                                                   <  1:187623/138‑1 (MQ=255)
                                                                    cGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAAATTTTAACCCCCCCTGTTGTGTgg                                                                                >  2:84432/1‑150 (MQ=25)
                                                                         acaccCGCGCTTTTACCCAACGCCATGTGAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                          <  2:398284/151‑1 (MQ=255)
                                                                            cccGCGCTTTTATCCACCGCCATGAGAAAAACTGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                         <  1:696502/148‑1 (MQ=255)
                                                                                cgcATTTACCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                                    <  2:763253/150‑1 (MQ=255)
                                                                                     taccccacctccATGTGAAAAAAGCCCTTCGCGGCGGCTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                               <  1:696805/140‑1 (MQ=255)
                                                                                      tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAA‑G‑AAAAAAATTTTAACCCCCCATGTTGTGTGGGCATCGCTGCTggggaaa                                                              >  2:40661/1‑145 (MQ=25)
                                                                                         ccAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑G‑AAAAAAATTTTAAACCCCCATGTTGTGTGGGGAACGCTGCTGGGAAC‑aa                                                            >  2:1005114/1‑148 (MQ=25)
                                                                                         ccAACCCCATTTGCAAAATTCCTTTCGCGGGGGCTGTCTGCCACCTATTTTTCCCCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                                 <  2:984284/144‑1 (MQ=255)
                                                                                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAGCTTTTAACCCCCCCTGTTTTGTGGGGATCGCTGCTGGGTAC‑AAagtg                                                        >  2:1012915/1‑149 (MQ=37)
                                                                                                cATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGG‑G‑AAAAAACTTTAAATCCGTCCTTTTGTGTGGCAAAAGCTGCTGGGTAC‑aa                                                            >  1:763907/1‑141 (MQ=37)
                                                                                                cATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAAATTTTAACCCCCCCTGTTTTTTGGGGATCGCTGCTGGGAaa‑aaaatgggaa                                                    >  2:821808/1‑141 (MQ=25)
                                                                                                   gtgAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAggg                                                 <  2:827136/149‑1 (MQ=255)
                                                                                                       aaaaaTGCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGt                                             <  2:676421/149‑1 (MQ=255)
                                                                                                       aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGa                                                <  2:522369/146‑1 (MQ=255)
                                                                                                       aaaaaTGCCTCCGCGGCGTCTTCCTCCCAGTTTTTTTCCCCCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGa                                                <  1:514051/146‑1 (MQ=255)
                                                                                                       aaaaaTGCCTCCGCGGCGGCGGTCTGCCGGCTTTTTTCCCGCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGt                                             <  2:782480/149‑1 (MQ=255)
                                                                                                        aaatTGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTgcg                                          <  2:125560/151‑1 (MQ=255)
                                                                                                          aaTGCCTCCGCGGCGGCTGTCTCCCAGTTATTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTgcggc                                        <  2:629024/151‑1 (MQ=255)
                                                                                                              ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGAGG‑AAAAACCTTTAAACCCG‑CCTGTTGTGGGGGGAACGCTGCTGGTGAC‑AATGAGGGAAggggggggggggg                                       >  2:246475/1‑138 (MQ=255)
                                                                                                                ttCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAACTTTAA‑CCCCCCATGGTGTGTGGGGAACGCTGCTTGTGAA‑AAAATGGGAAAGGaaaggggggg                                       >  2:378283/1‑138 (MQ=255)
                                                                                                                   gcggcgGCTGTCTGCGAGCTTTTTTCCCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAg                                 <  1:104186/149‑1 (MQ=255)
                                                                                                                        ggttGTCTGCCAGTTTTTTTCCCCCCGAAACAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                            <  1:878376/146‑1 (MQ=255)
                                                                                                                           tGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAACTTTTAACCCGCCCTGTTTTGTGGGGATCGCTGCTGGGGAA‑AATGTGGGAAGGGGAGGggggcggaaaagcacc                             >  2:274721/1‑135 (MQ=255)
                                                                                                                                 ccAGTTTTTTTCCCCCCCAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg                    <  1:151999/148‑1 (MQ=255)
                                                                                                                                 ccAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGG‑G‑AAATAAATTTAAACCCGTCCTGTTTTGTGGGGAACGCTGCTGGTGAC‑AATATGGGAAGGGGAGGGGGGCGGATCAGCTCCAGAACcgc                     >  1:730047/1‑147 (MQ=255)
                                                                                                                                    gCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGTGCCAGA‑GAAAAAAACTTTAA‑CCCCCCACGTTTTGTGGGGATTGCTGCTGGGGAA‑AAAGTGGGAAGGGGAGGGGGGGGGATCACCCCCAGAACCGCGCTa                 >  1:298083/1‑148 (MQ=255)
                                                                                                                                     cTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAC‑G‑AAAAAACTTTTAAACCGCCCTCTTTTTTGGCAAGCGCTGCTCGTTAC‑AATGGGGGAAGGGGAGGGCGGGGGATCAACTCCCGAACCCCGCCAAtt              >  2:464315/1‑150 (MQ=255)
                                                                                                                                      ttttttcCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt               <  1:399471/146‑1 (MQ=255)
                                                                                                                                        ttttcCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg           <  1:711346/150‑1 (MQ=255)
                                                                                                                                        ttttcCCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt              <  2:741399/147‑1 (MQ=255)
                                                                                                                                          ttcccgccgCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt         <  2:87382/147‑1 (MQ=255)
                                                                                                                                             ccgccgAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc        <  2:513697/148‑1 (MQ=255)
                                                                                                                                             cccccgaAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc     <  2:792076/148‑1 (MQ=255)
                                                                                                                                              cgccgcAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGG‑GAAAAAAAATTTAATCCC‑CCCTGTTTTGTGGGCATAGCTTCTGGGAACAAAATTGGAAAGGGAGGGGGGGGGAA‑CAACTCCACAACCGCGCCAATTTAGGTCGcc     >  2:747329/1‑150 (MQ=255)
                                                                                                                                                ccgcAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa  <  2:26159/151‑1 (MQ=255)
                                                                                                                                                ccgaAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc   <  1:1005014/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA  >  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: ‑