Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 19,634 (C)5→6 intergenic (+14/+177) nhaR → / ← insB1 DNA‑binding transcriptional activator NhaR/IS1 family protein InsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091319,6341.C85.9% 193.8 / 26.5 64intergenic (+14/+177)nhaR/insB1DNA‑binding transcriptional activator NhaR/IS1 family protein InsB
Reads supporting (aligned to +/- strand):  ref base . (8/1);  new base C (19/36);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 3.15e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCT‑C‑CC‑C‑C‑AAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATA  >  NC_000913/19343‑19819
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                          
ggTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CCAC‑C‑AAA‑GTTa                                                                                                                                                                                    >  2:460364/1‑300 (MQ=255)
     aaCTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTg                                                                                                                                                                               <  1:675545/300‑1 (MQ=255)
            aGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGTTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTggggg                                                                                                                                                                           <  2:580905/297‑1 (MQ=255)
             gggATTAAACGTAGAAATCCTCTGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTTCGAGGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGAGAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTg                                                                                                                                                                               <  2:45555/292‑1 (MQ=255)
              ggATTAAACGTAGAAATCCTCGGCGATGTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACACTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGGTTCAGCACCCTGCGGGACAGCGAATCTGCAATACGGATTATGCTGCGCTTTT‑AGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTg                                                                                                                                                                               <  2:45980/290‑1 (MQ=255)
              ggATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGagata                                                                                                                                                                      >  1:396835/1‑300 (MQ=255)
                aTTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGagataga                                                                                                                                                                    <  1:40877/300‑1 (MQ=255)
                  tAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGAtt                                                                                                                                                                  <  1:563860/300‑1 (MQ=255)
                    aaCGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCANGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAg                                                                                                                                                                <  2:498522/300‑1 (MQ=255)
                       gtataaattctttgCGAGTTTGATGATGCCGCGTTGATGGAAGCTTTTGGTGCGAGGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTg                                                                                                                                                             <  2:395198/287‑1 (MQ=255)
                            aaTCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACAt                                                                                                                                                        <  1:198016/300‑1 (MQ=255)
                            aaTCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACAt                                                                                                                                                        <  2:480596/300‑1 (MQ=255)
                                 tCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGTAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACAtaac                                                                                                                                                     <  2:33872/298‑4 (MQ=255)
                                 tCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACAtaac                                                                                                                                                     <  1:421567/298‑4 (MQ=255)
                                 tCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTNTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACAtaac                                                                                                                                                     <  1:507668/298‑4 (MQ=255)
                                 tCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTGTGGGGCGAGGCACAGTTCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCTGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACAtaac                                                                                                                                                     <  2:363309/298‑4 (MQ=255)
                                      gcgTTTGATGATTCCGCGTTGATGACAGCTTTTTGTGCGATGCACAATGCAATCTTCGTTGCCCCAACTCTTTATGCATATGACTTTTATGCCTAGAAAACTGTCTTAGAAATTGTTCGCGTCNAGAATGTGATGTAAGAATACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGGTTATTCTGCACTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACAt                                                                                                                                                        <  2:213839/288‑1 (MQ=255)
                                       aGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGAttt                                                                                                                                             <  1:605972/300‑1 (MQ=255)
                                            gatgatGCCGCGTTGATGAAAGCTTTTTGTGCGATGCACAATGCAATCTTCTTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGGCGAGAATGTGATGGAGGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACt                                                                                                                                        <  2:374582/300‑1 (MQ=255)
                                            gatgatGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCTGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGTAGATAGATTAGTTGTACATAACCACGATTTTGACt                                                                                                                                        <  1:535113/300‑1 (MQ=255)
                                            gatgatGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACNTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACt                                                                                                                                        <  1:148204/300‑1 (MQ=255)
                                                     gcgttGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTAC‑CC‑C‑C‑AAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTAGGCTCatt                                                                                                                               >  2:639316/4‑300 (MQ=255)
                                                     gcgttGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTAC‑AC‑C‑C‑AAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCatt                                                                                                                               >  1:287874/4‑300 (MQ=255)
                                                     gcgttGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CCAC‑C‑AAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCatt                                                                                                                               >  2:313669/4‑300 (MQ=255)
                                                     gcgttGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCatt                                                                                                                               >  1:287904/4‑300 (MQ=255)
                                                     gcgttGATGAAAGCTCTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCGTGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCatt                                                                                                                               >  2:313744/4‑300 (MQ=255)
                                                       gttGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCattat                                                                                                                             <  1:350239/299‑1 (MQ=255)
                                                       gttGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCattat                                                                                                                             <  2:464824/299‑1 (MQ=255)
                                                            tgaAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGcc                                                                                                                        <  1:33061/300‑1 (MQ=255)
                                                                  cTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACCTAACCACGATTTTGACTCGGCTCATTATTTGCCCGCtt                                                                                                                   >  2:172837/1‑299 (MQ=255)
                                                                  cTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACAACGATTTTGACCCGGCTCATTATTTGCCCGCtt                                                                                                                   >  2:174144/1‑299 (MQ=255)
                                                                      tGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGTAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGAc                                                                                                              <  2:450653/300‑1 (MQ=255)
                                                                          gCGATGCACAATGTAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTAGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTAAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CAAC‑C‑AAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGGCAtta                                                                                                          >  2:630962/1‑299 (MQ=255)
                                                                          gCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑CACC‑C‑C‑AAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACAtta                                                                                                          >  1:59769/1‑299 (MQ=255)
                                                                            gATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCTTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACAttat                                                                                                         <  1:220130/299‑3 (MQ=255)
                                                                             aTGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTAttt                                                                                                       >  1:443753/1‑300 (MQ=255)
                                                                               gCACAATGCAAGCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTTTTAGAAACTGTTCGCTTCGAGAATGTGATGGGAGAGTACTATGCTATTTTTGCTGAGCTGATGATTCAGCACCCGGCGGTACAGCGAAGCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGTGGCGCGTTAATCGGCACCT‑C‑CC‑C‑C‑AAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACAt                                                                                                                                                        <  2:371069/249‑1 (MQ=255)
                                                                                   aaTGCAATCTTTGTTTCCCCAACGTTTTATGCATATGACTTTTATGCCGATAAAACTTTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGTGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTTAGACATTATTTCCata                                                                                                  <  2:276835/299‑1 (MQ=255)
                                                                                                         cGCTTTATGCATATGANTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGNCATTATTTCCATATGTACNCGGGCGAATAAATaga                                                                            >  1:586918/1‑299 (MQ=255)
                                                                                                         cGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATaga                                                                            >  1:649944/1‑299 (MQ=255)
                                                                                                                      tGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCTAGAATGTGATGTAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTGTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTAc                                                               <  1:261653/299‑1 (MQ=255)
                                                                                                                      tGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTAc                                                               <  1:571742/299‑1 (MQ=255)
                                                                                                                         cTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTNTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTActtc                                                            >  2:502440/1‑299 (MQ=255)
                                                                                                                            ttATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTActtccttc                                                        <  1:35412/300‑1 (MQ=255)
                                                                                                                             tATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGNTGATTCAGCACCAGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCa                                                       <  1:495954/300‑1 (MQ=255)
                                                                                                                                             tCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCTGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTTAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAag                                       <  2:119297/300‑1 (MQ=255)
                                                                                                                                                     aTTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTaa                               <  1:555185/300‑1 (MQ=255)
                                                                                                                                                          tCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCACTATTTGCCCGCTTTAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCt                          <  2:608804/300‑1 (MQ=255)
                                                                                                                                                                      tgtgATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCttttaccatact               >  1:543713/1‑288 (MQ=255)
                                                                                                                                                                        tgATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCATCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGTAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATggg                                                   <  1:440288/261‑1 (MQ=255)
                                                                                                                                                                        tgATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCTGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCAGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATggg                                                   <  1:439884/261‑1 (MQ=255)
                                                                                                                                                                        tgATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CCAC‑C‑AAA‑GTTACGGTGGGGTAGGTAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATggg                                                   >  2:439884/1‑261 (MQ=255)
                                                                                                                                                                        tgATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTTAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATggg                                                   <  1:439766/261‑1 (MQ=255)
                                                                                                                                                                        tgATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATggg                                                   >  2:440288/1‑261 (MQ=255)
                                                                                                                                                                        tgATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATggg                                                   >  2:439766/1‑261 (MQ=255)
                                                                                                                                                                        tgATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATggg                                                   >  2:439273/1‑261 (MQ=255)
                                                                                                                                                                        tgATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATggg                                                   <  1:439273/261‑1 (MQ=255)
                                                                                                                                                                           tGGAAGAGTACCATGCNATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCNGAAAAGAGTAGTAATTGcntttaccatactggattt         >  1:585580/1‑282 (MQ=255)
                                                                                                                                                                             gAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGGCATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAGGAGTAGTAATTGCttttaccatactggatttca       >  1:618976/1‑281 (MQ=255)
                                                                                                                                                                             gAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCttttaccatcctggatttca       >  2:100006/1‑281 (MQ=255)
                                                                                                                                                                             gAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCttttaccatactggatttca       >  1:42617/1‑281 (MQ=255)
                                                                                                                                                                             gAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCCGAAAAGAGTAGTAATTGCttttaccatactggatttca       >  2:106894/1‑281 (MQ=255)
                                                                                                                                                                                  gTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑C‑CCAAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCttttaccatactggatttcattcgt  >  1:213056/1‑276 (MQ=255)
                                                                                                                                                                                                              atTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCT‑C‑CC‑CAC‑AAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGaataaat                                                                               >  1:640174/1‑195 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                          
GGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCT‑C‑CC‑C‑C‑AAA‑GTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATA  >  NC_000913/19343‑19819

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