Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,147,347 A→T 7.1% A1072A (GCT→GCA rne ← ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,147,3470AT 7.1% 299.3 / 10.6 99A1072A (GCT→GCArneribonuclease E
Reads supporting (aligned to +/- strand):  ref base A (52/40);  new base T (4/3);  total (56/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 103 total aligned reads displayed.

GTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAGGTGCAGG  >  NC_002947/2147201‑2147496
                                                                                                                                                  |                                                                                                                                                      
gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                                                                                                                     >  1:1377401/1‑150 (MQ=255)
    aCCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCtt                                                                                                                                                 >  2:1522450/1‑150 (MQ=255)
      ccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                                                                                                                     <  1:1216185/144‑1 (MQ=255)
      ccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                                                                                                                     >  2:1216185/1‑144 (MQ=255)
      ccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTcc                                                                                                                                               >  2:730409/1‑150 (MQ=255)
          gAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtgtg                                                                                                                                           >  2:1544625/1‑150 (MQ=255)
           aaGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGc                                                                                                                                          <  2:451978/150‑1 (MQ=255)
               ggCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCCCTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCGTGGGAAGCTGCAGCAGCTTCCTGTGCCGCt                                                                                                                                      <  2:1403817/150‑1 (MQ=255)
                gCCCTTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCCGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCtt                                                                                                                                     >  2:1233271/1‑150 (MQ=255)
                gCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGc                                                                                                                                          <  2:2364721/145‑1 (MQ=255)
                gCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGc                                                                                                                                          >  1:2364721/1‑145 (MQ=255)
                gCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCtt                                                                                                                                     >  2:1480949/1‑150 (MQ=255)
                gCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCtt                                                                                                                                     >  1:1933124/1‑150 (MQ=255)
                gCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCtt                                                                                                                                     >  2:1712128/1‑150 (MQ=255)
                  ccGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTg                                                                                                                                   >  2:2265346/1‑150 (MQ=255)
                   cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                                                                                                                  <  2:2272832/150‑1 (MQ=255)
                           aaaCCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCTGCAGCTTccca                                                                                                                                             <  1:2336793/131‑3 (MQ=255)
                           aaaCCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCTGCAGCTTccca                                                                                                                                             >  2:2336793/1‑129 (MQ=255)
                                     aaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTAATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTccc                                                                                                               >  2:1462425/1‑150 (MQ=255)
                                        ccccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCttt                                                                                                                                    >  1:93758/1‑127 (MQ=255)
                                        ccccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTAGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCt                                                                                                            >  1:384772/1‑150 (MQ=255)
                                         cccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt                                                                                                           <  2:1377401/150‑1 (MQ=255)
                                          ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt                                                                                                           <  1:705262/149‑1 (MQ=255)
                                          ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt                                                                                                           >  2:705262/1‑149 (MQ=255)
                                           cGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgc                                                                                                         >  2:167814/1‑150 (MQ=255)
                                           cGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgc                                                                                                         >  2:1858240/1‑150 (MQ=255)
                                                 tttGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGAc                                                                                                   <  1:2425472/150‑1 (MQ=255)
                                                  ttGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACg                                                                                                  <  1:2045044/150‑1 (MQ=255)
                                                    gCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTAATGTGCAGCATTGGCAGTAGCCTCGGCCTCCCGCTTGCGGCGACGCa                                                                                                <  1:1233271/150‑1 (MQ=255)
                                                            ggggCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                                                                                                                                <  2:2192274/95‑1 (MQ=255)
                                                            ggggCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                                                                                                                                <  2:2194426/95‑1 (MQ=255)
                                                            ggggCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                                                                                                                                >  1:2194426/1‑95 (MQ=255)
                                                                cTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGt                                                                                    >  2:381994/1‑150 (MQ=255)
                                                                 tttttCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTcc                                                                                                                <  1:91697/121‑1 (MQ=255)
                                                                 tttttCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTcc                                                                                                                >  2:91697/1‑121 (MQ=255)
                                                                 tttttCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTc                                                                                   >  1:6371/1‑150 (MQ=255)
                                                                  ttttCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCAAATGTTTTAGCAGCCGCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCg                                                                                  <  2:1530855/150‑1 (MQ=255)
                                                                      cattcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCg                                                                              <  2:1864023/150‑1 (MQ=255)
                                                                       attcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCgg                                                                             >  1:130655/1‑150 (MQ=255)
                                                                        ttcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                            >  1:1816071/1‑150 (MQ=255)
                                                                        ttcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                            <  1:1480949/150‑1 (MQ=255)
                                                                            gCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTgg                                                                                      <  1:231170/136‑1 (MQ=255)
                                                                            gCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTgg                                                                                      >  2:231170/1‑136 (MQ=255)
                                                                              tttGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTTTCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACgg                                                                      >  2:1891785/1‑150 (MQ=255)
                                                                              tttGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGACGTCTCGTTCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACgg                                                                      <  1:1878465/150‑1 (MQ=255)
                                                                                  gCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGt                                                                  <  1:1544625/150‑1 (MQ=255)
                                                                                   cTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGtt                                                                 >  1:2019361/1‑150 (MQ=255)
                                                                                    tATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTg                                                                <  2:2019361/150‑1 (MQ=255)
                                                                                     aTCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTgg                                                               <  1:1462425/150‑1 (MQ=255)
                                                                                        aGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCa                                                            >  1:489864/1‑150 (MQ=255)
                                                                                        aGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCa                                                            >  1:486140/1‑150 (MQ=255)
                                                                                               atgaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGAc                                                     >  1:12490/1‑150 (MQ=255)
                                                                                                tgaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACa                                                    >  2:1583069/1‑150 (MQ=255)
                                                                                                tgaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACa                                                    >  1:2231545/1‑150 (MQ=255)
                                                                                                 gaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCa                                                                         >  1:2208899/1‑128 (MQ=255)
                                                                                                 gaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCa                                                                         <  2:2208899/128‑1 (MQ=255)
                                                                                                  aGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGt                                                  <  2:2279293/150‑1 (MQ=255)
                                                                                                  aGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGt                                                  <  2:2279160/150‑1 (MQ=255)
                                                                                                    gTTTGTGCTCATCGGCGGCTTCCAGGGATGGTTGCTGGGAAGCTGCAGCAGCTTCCTGTGCGGCGTTGGCAGCAACCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGCCGTTCGGTGCACGGCCGTTGGCCAGCATGGCAGTgg                                                >  2:207961/1‑150 (MQ=255)
                                                                                                     tttGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGc                                               <  1:1858240/150‑1 (MQ=255)
                                                                                                     tttGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGc                                               <  1:2197288/150‑1 (MQ=255)
                                                                                                     tttGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTAGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGc                                               <  1:1712128/150‑1 (MQ=255)
                                                                                                          gCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGa                                                      <  2:1632947/138‑1 (MQ=255)
                                                                                                          gCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGa                                                      >  1:1632947/1‑138 (MQ=255)
                                                                                                           cTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                                         >  2:421329/1‑150 (MQ=255)
                                                                                                           cTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                                         <  1:877213/150‑1 (MQ=255)
                                                                                                            tcctCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCATCGGTAAGCAGCTGGGCCTCCAGCATGCGGCGACGCACTTCCCGTTGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCg                                        <  1:167814/147‑1 (MQ=255)
                                                                                                            tCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCg                                        >  1:477082/1‑150 (MQ=255)
                                                                                                             cATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGa                                       <  1:1522450/150‑1 (MQ=255)
                                                                                                             cATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGa                                       <  2:1816071/150‑1 (MQ=255)
                                                                                                              aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCCTGa                                                      <  2:1603262/134‑1 (MQ=255)
                                                                                                              aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGa                                                      >  1:1607737/1‑134 (MQ=255)
                                                                                                              aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGa                                                      >  1:1603262/1‑134 (MQ=255)
                                                                                                                cggcggCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCa                                    <  2:747473/150‑1 (MQ=255)
                                                                                                                    ggcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCtt                                                                                                                                     >  2:1453825/1‑50 (MQ=255)
                                                                                                                    ggcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCtt                                                                                                                                     <  1:1453825/50‑1 (MQ=255)
                                                                                                                       ttCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGa                                                                                                    >  1:1498396/1‑79 (MQ=255)
                                                                                                                       ttCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGa                                                                                                    <  2:1498396/79‑1 (MQ=255)
                                                                                                                       ttCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgcagg                             >  2:442999/1‑150 (MQ=255)
                                                                                                                           aGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTc                                                                                            >  2:2064765/1‑83 (MQ=255)
                                                                                                                           aGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTc                                                                                            <  1:2064765/83‑1 (MQ=255)
                                                                                                                           aGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgcaggtgca                         >  1:1629208/1‑150 (MQ=255)
                                                                                                                               cTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGtt                     >  1:1705349/1‑150 (MQ=255)
                                                                                                                                tGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTc                    >  1:823111/1‑150 (MQ=255)
                                                                                                                                   ttCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGac                 >  1:799234/1‑150 (MQ=255)
                                                                                                                                   ttCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGac                 >  1:1322227/1‑150 (MQ=255)
                                                                                                                                     cacaggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGacca               >  1:2178144/5‑150 (MQ=255)
                                                                                                                                      acaggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTgg                                                                                      >  2:1988617/4‑78 (MQ=255)
                                                                                                                                      acaggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTgg                                                                                      <  1:1988617/75‑1 (MQ=255)
                                                                                                                                       tGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccacc             >  2:2225466/1‑150 (MQ=255)
                                                                                                                                       tGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccacc             >  1:1987028/1‑150 (MQ=255)
                                                                                                                                        gggAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCg            >  1:1250854/1‑150 (MQ=255)
                                                                                                                                        aggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGc                                                                     >  2:1945854/2‑93 (MQ=255)
                                                                                                                                        aggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGc                                                                     <  1:1945854/92‑1 (MQ=255)
                                                                                                                                           aaGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAg         <  2:1705349/150‑1 (MQ=255)
                                                                                                                                              cTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCAGCGCAGGTg      >  2:1314775/1‑150 (MQ=255)
                                                                                                                                              cTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAGGTg      >  1:1890339/1‑150 (MQ=255)
                                                                                                                                               tGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgcaggtg                           >  1:2113316/1‑128 (MQ=255)
                                                                                                                                               tGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgcaggtg                           <  2:2113316/128‑1 (MQ=255)
                                                                                                                                                  agcagcTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAGGTGCAgg  >  1:2023858/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAGGTGCAGG  >  NC_002947/2147201‑2147496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: