Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,993,337 Δ1 bp 54.1% coding (2185/2547 nt) cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,993,3370T.54.1% 198.1 / 278.4 111coding (2185/2547 nt)cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base . (34/26);  total (58/53)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 112 total aligned reads displayed.

CGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGA  >  NC_000913/3993190‑3993482
                                                                                                                                                   |                                                                                                                                                 
cgcgTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAAt                                                                                                                                                   >  2:961748/1‑148 (MQ=255)
     ttCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACcg                                                                                                                                              >  1:24724/1‑148 (MQ=255)
           tCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAGATTACCGGCggtgg                                                                                                                                       >  3:392198/1‑149 (MQ=255)
           tCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtg                                                                                                                                        >  1:216376/1‑148 (MQ=255)
           tCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCggtggt                                                                                                                                      <  2:770693/148‑1 (MQ=255)
            cAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtgg                                                                                                                                       >  1:909586/1‑148 (MQ=255)
             aaaCCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGgtggtg                                                                                                                                     >  4:35101/1‑149 (MQ=255)
             aaaCCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGa                                                                                                                                   >  3:888421/1‑149 (MQ=255)
              aaCCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGtggtgg                                                                                                                                    >  4:635710/1‑149 (MQ=255)
              aaCCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGa                                                                                                                                   <  1:304737/148‑1 (MQ=255)
                   gggggTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCt                                                                                                                              >  2:387800/1‑148 (MQ=255)
                      ggTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTg                                                                                                                           <  1:620026/148‑1 (MQ=255)
                       gttgttCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTg                                                                                                                           >  3:660977/1‑149 (MQ=255)
                          gttCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGcc                                                                                                                         <  1:961748/148‑1 (MQ=255)
                          gttCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCa                                                                                                                        >  3:686083/1‑149 (MQ=255)
                            tcttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGc                                                                                                                      >  3:325605/1‑149 (MQ=255)
                              ttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGa                                                                                                                    >  4:476440/1‑149 (MQ=255)
                               tcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAg                                                                                                                  >  2:945987/1‑148 (MQ=255)
                               tcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAg                                                                                                                  <  2:83755/148‑1 (MQ=255)
                                   aCGCCTTAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAggg                                                                                                                >  2:518743/1‑148 (MQ=255)
                                   aCGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGa                                                                                                               <  3:241770/149‑1 (MQ=255)
                                   aCGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGat                                                                                                              >  2:863539/1‑148 (MQ=255)
                                    cGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGat                                                                                                              <  4:877477/149‑1 (MQ=255)
                                      ccTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatc                                                                                                             <  1:518743/148‑1 (MQ=255)
                                      ccTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatcat                                                                                                           >  1:685734/1‑148 (MQ=255)
                                       cTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATcc                                                                                                         >  3:112562/1‑149 (MQ=255)
                                        tGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatcat                                                                                                           <  2:664097/148‑1 (MQ=255)
                                        tGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGCAGGGATCATcc                                                                                                         >  2:23010/1‑148 (MQ=255)
                                        tGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATcc                                                                                                         >  2:523365/1‑148 (MQ=255)
                                         gAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAg                                                                                                       >  4:639297/1‑149 (MQ=255)
                                           aTGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGt                                                                                                      >  2:695433/1‑148 (MQ=255)
                                               aTCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                                                 >  3:36383/1‑149 (MQ=255)
                                                tCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATcc                                                                                                         >  1:643000/1‑140 (MQ=255)
                                                tCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATcc                                                                                                         <  2:643000/140‑1 (MQ=255)
                                                tCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGCTCttt                                                                                                 <  2:685734/148‑1 (MQ=255)
                                                 cGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCAGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttttt                                                                                               <  4:33034/149‑1 (MQ=255)
                                                  ggTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                                                 <  2:606211/148‑1 (MQ=255)
                                                  ggTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                                                 <  1:835198/148‑1 (MQ=255)
                                                  ggTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTc                                                                                              <  4:652950/149‑1 (MQ=255)
                                                      cAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCga                                                                                            >  3:745923/1‑149 (MQ=255)
                                                             tGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACg                                                                                     >  4:863580/1‑149 (MQ=255)
                                                                         tCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATg                                                                        <  1:387800/148‑1 (MQ=255)
                                                                          cGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAAt                                                                         >  2:230153/1‑148 (MQ=255)
                                                                          cGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAAt                                                                         >  2:169824/1‑148 (MQ=255)
                                                                             gTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCttt                                                                   <  4:602874/149‑1 (MQ=255)
                                                                               tttATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCtt                                                                    <  1:487569/148‑1 (MQ=255)
                                                                                 tATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTaa                                                                 <  3:476440/149‑1 (MQ=255)
                                                                                  aTGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTaa                                                                 <  2:842036/148‑1 (MQ=255)
                                                                                    ggCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCt                                                            <  3:346711/149‑1 (MQ=255)
                                                                                     gcgcGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAatat                                                              <  1:169824/148‑1 (MQ=255)
                                                                                        cgATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTa                                                           >  2:767879/1‑148 (MQ=255)
                                                                                           tttCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATtct                                                     >  3:866463/1‑149 (MQ=255)
                                                                                           tttCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATtct                                                     <  4:788493/149‑1 (MQ=255)
                                                                                              cGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTcgcg                                                   <  2:604510/148‑1 (MQ=255)
                                                                                                cATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGGCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcga                                                  <  4:779981/149‑1 (MQ=255)
                                                                                                cATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcg                                                   >  2:921296/1‑148 (MQ=255)
                                                                                                 aTAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcgacga                                               >  4:375443/1‑149 (MQ=255)
                                                                                                  tAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcgac                                                 >  2:958141/1‑148 (MQ=255)
                                                                                                  tAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcgacga                                               <  1:290970/148‑1 (MQ=255)
                                                                                                  tAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaa                                              >  4:391739/1‑149 (MQ=255)
                                                                                                   aaCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaaa                                             >  3:201879/1‑149 (MQ=255)
                                                                                                   aaCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaa                                              <  1:292639/148‑1 (MQ=255)
                                                                                                    aCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                                                 >  1:214963/1‑96 (MQ=255)
                                                                                                    aCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                                                 <  2:214963/96‑1 (MQ=255)
                                                                                                     cAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaa                                              <  2:909586/148‑1 (MQ=255)
                                                                                                     cAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaa                                              <  1:965835/148‑1 (MQ=255)
                                                                                                       aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAg                                            <  2:623586/148‑1 (MQ=255)
                                                                                                       aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGaa                                                                                                                   >  1:945310/1‑75 (MQ=255)
                                                                                                       aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGaa                                                                                                                   <  2:945310/75‑1 (MQ=255)
                                                                                                       aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATc                                                             <  4:476753/129‑1 (MQ=255)
                                                                                                       aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATc                                                             >  3:476753/1‑129 (MQ=255)
                                                                                                       aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCaa                                         >  3:948764/1‑149 (MQ=255)
                                                                                                          tGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCaa                                         <  1:767879/148‑1 (MQ=255)
                                                                                                           gCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCgg                                     <  3:289436/149‑1 (MQ=255)
                                                                                                           gCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCg                                      <  1:23010/148‑1 (MQ=255)
                                                                                                                  cTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTg                                 <  1:940647/148‑1 (MQ=255)
                                                                                                                  cTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTg                                 <  2:24724/148‑1 (MQ=255)
                                                                                                                   tGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAAACGGGTTGAGGt                             >  4:769357/1‑149 (MQ=255)
                                                                                                                       aGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGta                            <  1:132379/148‑1 (MQ=255)
                                                                                                                       aGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGta                            <  2:216376/148‑1 (MQ=255)
                                                                                                                           cAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcacc                     <  4:360474/149‑1 (MQ=255)
                                                                                                                           cAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcac                      >  1:39385/1‑148 (MQ=255)
                                                                                                                              gTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcacca                    <  4:366712/149‑1 (MQ=255)
                                                                                                                              gTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACt                  >  4:668331/1‑149 (MQ=255)
                                                                                                                              gTTGAAACCAATCACGTCAAA‑‑ACCGGC‑‑‑GGTGGAGGGCTTTGTCACCGAAGAGAGCATCCAGCACTTTTTCGAAGAAACGCAAGACGACAAGAGCGTCAATATCTACATTCTCGTCGAAAGAAACAAGCGTGAGGTATATTACCACTGct               >  4:900308/1‑148 (MQ=255)
                                                                                                                               ttGAAACCGATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAg                                                                                                                  >  2:612070/1‑52 (MQ=255)
                                                                                                                               ttGAAACCGATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAg                                                                                                                  <  1:612070/52‑1 (MQ=255)
                                                                                                                               ttGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACt                  >  1:917494/1‑148 (MQ=255)
                                                                                                                               ttGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACt                  >  2:321469/1‑148 (MQ=255)
                                                                                                                                  aaaCCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGa              >  4:615731/1‑149 (MQ=255)
                                                                                                                                     ccAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAgg           <  4:147478/149‑1 (MQ=255)
                                                                                                                                      cAATCACGTCAAA‑‑ACCGGCGGTGGTGTACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGgc          >  4:30780/1‑149 (MQ=255)
                                                                                                                                       aaTCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGgca         <  4:875850/149‑1 (MQ=255)
                                                                                                                                        aTCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCg                                      <  2:798930/119‑1 (MQ=255)
                                                                                                                                        aTCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCg                                      >  1:798930/1‑119 (MQ=255)
                                                                                                                                             gTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAgg                                                                                                                 <  4:467159/41‑1 (MQ=255)
                                                                                                                                              tCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCaaa    <  3:664919/149‑1 (MQ=255)
                                                                                                                                               cAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCaaa    >  1:87065/1‑148 (MQ=255)
                                                                                                                                                 aaTTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAaga  <  2:71327/148‑1 (MQ=255)
                                                                                                                                                 aaTTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAaga  <  2:763742/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGA  >  NC_000913/3993190‑3993482

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