Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,147,347 A→T 7.7% 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.7% 245.4 / 9.8 78A1072A (GCT→GCArneribonuclease E
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base T (4/2);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

CGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAG  >  NC_002947/2147198‑2147489
                                                                                                                                                       |                                                                                                                                              
cGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgca                                                                                                                                                <  2:440769/150‑1 (MQ=255)
 gggTTGACCCGCGAAGAGGCACGTCAGGAAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTGTTTCATTCCGCTTTGGCTATCAGCCGCGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGATGGTTCCTGGGAAGCTgcag                                                                                                                                               <  2:742775/150‑1 (MQ=255)
        cccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                                                                                                                        <  2:114213/150‑1 (MQ=255)
          cgcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCt                                                                                                                                      >  1:1951378/1‑150 (MQ=255)
                      cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGCGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                                                                                                          >  1:1250216/1‑150 (MQ=255)
                      cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCCTCGTCGAGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                                                                                                          <  1:1955586/150‑1 (MQ=255)
                           ggCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                                                                                                                     >  2:2150138/1‑150 (MQ=255)
                                aCCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTc                                                                                                                >  1:1974297/1‑150 (MQ=255)
                                 ccAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCg                                                                                                               >  2:2271440/1‑150 (MQ=255)
                                  cAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCgg                                                                                                              >  2:2324659/1‑150 (MQ=255)
                                      aaaaGCCCCGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                                                                                                          <  2:1509860/134‑1 (MQ=39)
                                      aaaaGCCCCGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                                                                                                          >  1:1509860/1‑134 (MQ=39)
                                         aGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTCATGAGATTTGTGCTCATCGGCGGCTTCCAGGGCTCGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCg                                                                                                       <  1:1060863/150‑1 (MQ=255)
                                            cccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                                                                                                             >  2:1225556/1‑109 (MQ=255)
                                            cccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                                                                                                             <  1:1225556/109‑1 (MQ=255)
                                            cccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtg                                                                                                                                     <  1:734909/117‑1 (MQ=255)
                                            cccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtg                                                                                                                                     >  2:734909/1‑117 (MQ=255)
                                            cccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTATGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt                                                                                                    >  2:810725/1‑150 (MQ=255)
                                            cccGCCAGGTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt                                                                                                    >  1:1177414/1‑150 (MQ=255)
                                             ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTg                                                                                                   >  2:1936512/1‑150 (MQ=255)
                                              cGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTg                                                                                                                           <  2:1103907/125‑1 (MQ=255)
                                              cGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTg                                                                                                                           >  1:1103907/1‑125 (MQ=255)
                                              cGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgc                                                                                                  <  1:2324659/150‑1 (MQ=255)
                                               gCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCGTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcg                                                                                                 >  1:1226435/1‑150 (MQ=255)
                                               gCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcag                                                                                                                       >  2:1987598/1‑128 (MQ=255)
                                               gCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcag                                                                                                                       <  1:1987598/128‑1 (MQ=255)
                                               gCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcg                                                                                                 <  2:1974297/150‑1 (MQ=255)
                                                        cACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCAc                                                                                        >  2:1258118/1‑150 (MQ=255)
                                                        cACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCAc                                                                                        <  2:1177414/150‑1 (MQ=255)
                                                        cACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATAGACGTTTATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCAc                                                                                        <  1:126798/150‑1 (MQ=255)
                                                         aCTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACt                                                                                       >  1:1715778/1‑150 (MQ=255)
                                                         aCTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACt                                                                                       <  1:2150138/150‑1 (MQ=255)
                                                         aCTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACt                                                                                       >  2:1084666/1‑150 (MQ=255)
                                                           tGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTc                                                                                     >  1:1784841/1‑150 (MQ=255)
                                                                gggCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTg                                                                                <  1:136339/150‑1 (MQ=255)
                                                                 ggCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTgg                                                                               <  2:571308/150‑1 (MQ=255)
                                                                      tttCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGt                                                                          >  1:2153297/1‑150 (MQ=255)
                                                                           ttcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                     <  1:1321309/150‑1 (MQ=255)
                                                                                cTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACg                                                                <  1:2271440/150‑1 (MQ=255)
                                                                                cTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACg                                                                >  2:2013401/1‑150 (MQ=255)
                                                                                 tttGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtgt                                                                                                                                    >  1:1738726/1‑81 (MQ=255)
                                                                                 tttGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtgt                                                                                                                                    <  2:1738726/81‑1 (MQ=255)
                                                                                 tttGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACgg                                                               >  1:1457899/1‑150 (MQ=255)
                                                                                    ggCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCa                                                                                         >  2:1700209/1‑121 (MQ=255)
                                                                                    ggCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCa                                                                                         <  1:1700209/121‑1 (MQ=255)
                                                                                     gCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCAAGGGCTGGTTCCTGGGAATCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGt                                                           >  1:1500894/1‑150 (MQ=255)
                                                                                       tATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTgg                                                                               <  1:1435849/128‑1 (MQ=255)
                                                                                       tATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTgg                                                                               >  2:1435849/1‑128 (MQ=255)
                                                                                        aTCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTgg                                                        >  1:1484137/1‑150 (MQ=255)
                                                                                          cAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGcc                                                      <  2:1484137/150‑1 (MQ=255)
                                                                                           cgCCGTGATGACTTTTGTGCTCATCTGCGACTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTTGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCa                                                     <  1:1244154/149‑1 (MQ=255)
                                                                                            gCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAg                                                    <  1:336578/150‑1 (MQ=255)
                                                                                            gCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAg                                                    <  2:1250216/150‑1 (MQ=255)
                                                                                                         gTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTgg                                         >  1:2126469/1‑150 (MQ=255)
                                                                                                               gCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGtt                                   >  2:1942588/1‑150 (MQ=255)
                                                                                                                cTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                                  <  1:2013401/150‑1 (MQ=255)
                                                                                                                 tCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCg                                 >  1:749104/1‑150 (MQ=255)
                                                                                                                      ggcggcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCATTGGCCAGCATGACAGTGGCAGGTTCGACCAc                            <  1:1084666/150‑1 (MQ=255)
                                                                                                                         ggcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                     >  1:575479/1‑106 (MQ=255)
                                                                                                                         ggcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                     <  2:575479/106‑1 (MQ=255)
                                                                                                                          gcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgca                        <  1:1936512/150‑1 (MQ=255)
                                                                                                                            ttCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAg                                                    >  2:134290/1‑120 (MQ=255)
                                                                                                                            ttCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAg                                                    <  1:134290/120‑1 (MQ=255)
                                                                                                                             tCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgcaggt                     <  2:1457899/150‑1 (MQ=255)
                                                                                                                               cAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgcaggtgc                   <  2:2039719/150‑1 (MQ=255)
                                                                                                                                    cTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgca                        >  2:446484/1‑140 (MQ=255)
                                                                                                                                    cTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgca                        <  1:446484/140‑1 (MQ=255)
                                                                                                                                        ttCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGac          >  1:832714/1‑150 (MQ=255)
                                                                                                                                           cTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCgg                                                                      <  2:407/87‑1 (MQ=255)
                                                                                                                                           cTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCgg                                                                      >  1:407/1‑87 (MQ=255)
                                                                                                                                           cTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccac       >  2:2010796/1‑150 (MQ=255)
                                                                                                                                           acaggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccac       >  1:1081615/4‑150 (MQ=255)
                                                                                                                                           acaggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAAAAGCATCGGCCGCCCGCTTGCGGCTACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccac       >  2:84392/4‑150 (MQ=255)
                                                                                                                                            tGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGa                                               <  1:1442857/109‑1 (MQ=255)
                                                                                                                                            tGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGa                                               >  2:1442857/1‑109 (MQ=255)
                                                                                                                                            tGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccacc      >  1:2084336/1‑150 (MQ=255)
                                                                                                                                            caggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccacc      >  2:2250284/3‑150 (MQ=255)
                                                                                                                                             gggAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCg     >  2:2294321/1‑150 (MQ=255)
                                                                                                                                             aggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCa                                                  >  1:105355/2‑105 (MQ=255)
                                                                                                                                             aggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCa                                                  <  2:105355/104‑1 (MQ=255)
                                                                                                                                               gAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCa   <  2:1081615/150‑1 (MQ=255)
                                                                                                                                                aaGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCATCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAg  <  1:871459/150‑1 (MQ=255)
                                                                                                                                                     gcagcaGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCg                                                                       >  2:2291765/1‑76 (MQ=255)
                                                                                                                                                     gcagcaGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCg                                                                       <  1:2291765/76‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                              
CGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAG  >  NC_002947/2147198‑2147489

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