Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,147,304 A→T 5.3% 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,3040AT 5.3% 229.8 / 2.3 76A1072A (GCT→GCA) rneribonuclease E
Reads supporting (aligned to +/- strand):  ref base A (32/40);  new base T (2/2);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCG  >  NC_002947/2147159‑2147443
                                                                                                                                                 |                                                                                                                                           
tttACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                                                                                                          >  4:69895/1‑149 (MQ=255)
     ccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCtt                                                                                                                                      >  2:296860/1‑148 (MQ=255)
     ccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                                                                                                                     >  7:238936/1‑149 (MQ=255)
      cgcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTcc                                                                                                                                    >  1:234651/1‑149 (MQ=255)
       gcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCt                                                                                                                                   <  3:160254/149‑1 (MQ=255)
        cgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCTGCAGCTTccca                                                                                                                                  <  3:187756/149‑3 (MQ=255)
         gAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtgt                                                                                                                                 <  1:283192/149‑1 (MQ=255)
                 ccGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCttt                                                                                                                         <  2:234651/149‑1 (MQ=255)
                      atgCAAACCAAAAAAAGCCCCGCCAGTTTTCACTGGCTGGGCTTTTTCATTCAGATTTGTCTATCAGCCGTGATGAGGTTTTTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgca                                                                                                                     <  7:47729/146‑1 (MQ=255)
                      aGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcag                                                                                                                    <  4:177579/149‑1 (MQ=255)
                        gCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                                                                                                                  <  7:7012/149‑1 (MQ=255)
                        gCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                                                                                                                  >  6:262964/1‑149 (MQ=255)
                          aaaCCAAAAAATGCCCCGCCAGTTTGCACTGGCGGGGCTTTCTAATTCAGCTTTTGCTATCAGCCGTTATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCagcagc                                                                                                                <  5:920/149‑1 (MQ=255)
                           aaCCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGcc                                                                                                               >  1:224662/1‑149 (MQ=255)
                              cAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCg                                                                                                            >  7:8407/1‑149 (MQ=255)
                              cAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCg                                                                                                            >  3:295318/1‑149 (MQ=255)
                                   aatgCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGAGGCTTCCAGGGCTGGTTCGTGTGAAGCTGCTGCAGCTTccc                                                                                                                                   <  4:231582/118‑2 (MQ=255)
                                   aaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCGTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCTGCAGCTTccc                                                                                                                                   >  3:227327/1‑120 (MQ=255)
                                   aaaGACCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCATTGGATATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCTGCAGCTTcca                                                                                                                                   >  3:231582/1‑120 (MQ=255)
                                     aGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTccc                                                                                                     <  3:69895/149‑1 (MQ=255)
                                        cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcag                                                                                                                                            <  7:301399/107‑1 (MQ=39)
                                        cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcag                                                                                                                                            >  8:301399/1‑107 (MQ=39)
                                          cGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTg                                                                                                >  1:147588/1‑149 (MQ=255)
                                            ccAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcg                                                                                              <  2:147588/149‑1 (MQ=255)
                                               gtgcaaaCTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcggcg                                                                                           <  4:135799/143‑1 (MQ=255)
                                                     aCTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCAc                                                                                     >  4:151157/1‑149 (MQ=255)
                                                         gCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTc                                                                                  <  3:151157/148‑1 (MQ=255)
                                                              gCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACACACTTCACGTgg                                                                            <  8:238936/149‑1 (MQ=255)
                                                                    tcattcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGt                                                                       <  1:160679/148‑1 (MQ=255)
                                                                         caGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                  <  6:263861/148‑1 (MQ=255)
                                                                          aGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTg                                                                 >  6:101608/1‑148 (MQ=255)
                                                                                 gCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGcc                                                          <  2:288900/148‑1 (MQ=255)
                                                                                     tCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTgg                                                     >  3:209436/1‑149 (MQ=255)
                                                                                            tgatgaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGGTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCAt                                              <  7:131556/149‑1 (MQ=255)
                                                                                              atgaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGa                                            <  5:262964/149‑1 (MQ=255)
                                                                                                 aGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACa                                          <  4:295318/148‑1 (MQ=255)
                                                                                                  ggTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACtt                                                                                   <  4:163053/106‑1 (MQ=255)
                                                                                                  ggTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACtt                                                                                   >  3:163053/1‑106 (MQ=255)
                                                                                                  ggTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGt                                        <  2:224662/149‑1 (MQ=255)
                                                                                                   gTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGtt                                                                      >  5:131686/1‑118 (MQ=255)
                                                                                                   gTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGtt                                                                      <  6:131686/118‑1 (MQ=255)
                                                                                                   gTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGtt                                                       <  4:294977/133‑1 (MQ=255)
                                                                                                   gTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGtt                                                       >  3:294977/1‑133 (MQ=255)
                                                                                                           tCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                               >  6:86423/1‑149 (MQ=255)
                                                                                                             aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTc                                                                         >  3:110006/1‑105 (MQ=255)
                                                                                                             aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTc                                                                         <  4:110006/105‑1 (MQ=255)
                                                                                                             aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAgg                                  >  6:276741/1‑144 (MQ=255)
                                                                                                             aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAgg                                  <  5:276741/144‑1 (MQ=255)
                                                                                                             aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                               <  5:86423/147‑1 (MQ=255)
                                                                                                             aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                               <  8:8407/147‑1 (MQ=255)
                                                                                                             aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                               <  4:209436/147‑1 (MQ=255)
                                                                                                             aTCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGa                             <  6:318548/149‑1 (MQ=255)
                                                                                                               cggcggCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGAc                            >  1:231102/1‑148 (MQ=255)
                                                                                                                      ttCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACg                                                                                        >  8:267518/1‑81 (MQ=255)
                                                                                                                      ttCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACg                                                                                        <  7:267518/81‑1 (MQ=255)
                                                                                                                      ttCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGttg                                      <  7:207471/131‑3 (MQ=255)
                                                                                                                      ttCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTgg                                      >  8:207471/1‑131 (MQ=255)
                                                                                                                      ttCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCACTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgcag                    >  3:50111/1‑149 (MQ=255)
                                                                                                                       tCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGc                                                <  1:155584/120‑1 (MQ=255)
                                                                                                                       tCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGc                                                >  2:155584/1‑120 (MQ=255)
                                                                                                                       tCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                               >  3:302977/1‑137 (MQ=255)
                                                                                                                       tCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                               <  4:302977/137‑1 (MQ=255)
                                                                                                                                 gTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGcc                                                   <  2:276453/107‑1 (MQ=255)
                                                                                                                                 gTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGcc                                                   >  1:276453/1‑107 (MQ=255)
                                                                                                                                 gTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgc                      <  5:286787/136‑1 (MQ=255)
                                                                                                                                 gTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgc                      >  6:286787/1‑136 (MQ=255)
                                                                                                                                 gTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTc          >  3:140346/1‑148 (MQ=255)
                                                                                                                                     cTGGGAAGCTGCAGCAGCTTCCTGTGCCACTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGacca     <  4:50111/149‑1 (MQ=255)
                                                                                                                                      tGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                  <  7:68558/87‑1 (MQ=255)
                                                                                                                                      tGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                  >  8:68558/1‑87 (MQ=255)
                                                                                                                                        ggAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCg                              >  1:114021/1‑121 (MQ=255)
                                                                                                                                        ggAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCg                              <  2:114021/121‑1 (MQ=255)
                                                                                                                                         gAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCg  >  8:99962/1‑148 (MQ=255)
                                                                                                                                          aaGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTTGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                               <  6:185980/118‑1 (MQ=255)
                                                                                                                                          aaGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                               >  5:185980/1‑118 (MQ=255)
                                                                                                                                          aaGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                               <  6:190527/118‑1 (MQ=255)
                                                                                                                                          aaGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTc                               >  5:190527/1‑118 (MQ=255)
                                                                                                                                              tGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGc                                     <  8:211367/108‑1 (MQ=255)
                                                                                                                                              tGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGc                                     >  7:211367/1‑108 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
TTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCG  >  NC_002947/2147159‑2147443

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