Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,503,465 G→T 6.3% E90* (GAG→TAG) ‡ rplM → 50S ribosomal protein L13

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,503,4650GT 6.3% 251.1 / 6.1 79E90* (GAG→TAG) ‡rplM50S ribosomal protein L13
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base T (3/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

AAACCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAGCCTCAAG  >  NC_002947/1503318‑1503612
                                                                                                                                                    |                                                                                                                                                    
aaaCCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCga                                                                                                                                                     <  7:348902/149‑1 (MQ=255)
    cAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGATCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCgaga                                                                                                                                                   <  4:254252/147‑1 (MQ=255)
    cAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGaaa                                                                                                                                                 >  3:254252/1‑148 (MQ=255)
    cAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAg                                                                                                                                                 >  8:219388/1‑149 (MQ=255)
        aTACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGa                                                                                                                                             <  4:78332/149‑1 (MQ=255)
        aTACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGa                                                                                                                                             <  5:209656/149‑1 (MQ=255)
           cacCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATc                                                                                                                                          <  8:154392/148‑1 (MQ=255)
            acCCCTCACGTTGACGCCGGCGACTACATCGTCGTTATCAACGCCGAGGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCg                                                                                                                                         >  1:90476/1‑149 (MQ=255)
            acCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGc                                                                                                                                        <  5:319905/149‑1 (MQ=255)
              cccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGcc                                                                                                                                       >  6:30485/1‑148 (MQ=255)
               ccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCaa                                                                                                                                     <  8:304357/149‑1 (MQ=255)
                     gTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGcccc                                                                                                                               >  2:38648/1‑149 (MQ=255)
                           acCGGCGTCTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttcctt                                                                                                                                      <  1:99068/137‑7 (MQ=255)
                           acCGGCGTCTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttcctt                                                                                                                                      >  2:99068/1‑131 (MQ=255)
                           acCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCg                                                                                                                         >  3:99908/1‑149 (MQ=255)
                             cGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCgtg                                                                                                                       >  4:191119/1‑149 (MQ=255)
                                 gACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTAt                                                                                                                   <  6:128561/149‑1 (MQ=255)
                                   cTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATc                                                                                                                  >  8:106604/1‑148 (MQ=255)
                                      cATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCCGCCAAGGCCCCTGAGCGTGTTATCGa                                                                                                                <  1:4341/148‑1 (MQ=255)
                                           tcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCGAAGGCCCCTGAGCGTGTTATCGAAACTGCg                                                                                                         <  1:29780/149‑1 (MQ=255)
                                             gtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCgg                                                                                                        <  5:145476/148‑1 (MQ=255)
                                                   aaCGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATc                                                                                                                  <  5:7478/132‑1 (MQ=255)
                                                   aaCGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATc                                                                                                                  >  6:7478/1‑132 (MQ=255)
                                                   aaCGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAgg                                                                                                 >  5:283951/1‑149 (MQ=255)
                                                   aaCGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAg                                                                                                  <  4:99908/148‑1 (MQ=255)
                                                    aCGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTGCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAgg                                                                                                 >  8:301680/1‑148 (MQ=255)
                                                      gCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCAt                                                                                              <  3:214528/149‑1 (MQ=255)
                                                              aGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCt                                                                                       <  2:90476/149‑1 (MQ=255)
                                                                 ttCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAg                                                                                    <  7:106604/149‑1 (MQ=255)
                                                                 ttCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAg                                                                                    <  3:191119/149‑1 (MQ=255)
                                                                  tCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGCCATGCTGCCTAAg                                                                                    >  6:260125/1‑148 (MQ=255)
                                                                        gACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGc                                                                             >  1:93909/1‑149 (MQ=255)
                                                                         aCTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGc                                                                             <  7:183728/148‑1 (MQ=255)
                                                                          cTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTg                                                                           <  1:38648/149‑1 (MQ=255)
                                                                          cTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTg                                                                           <  6:112240/149‑1 (MQ=255)
                                                                           tGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTgg                                                                          >  4:260271/1‑149 (MQ=255)
                                                                           tGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTgg                                                                          >  1:338813/1‑149 (MQ=255)
                                                                             gTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCgg                                                                                                        <  3:239583/117‑1 (MQ=255)
                                                                             gTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCgg                                                                                                        >  4:239583/1‑117 (MQ=255)
                                                                             gTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGt                                                                        >  1:113565/1‑149 (MQ=255)
                                                                                       ttCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATg                                                               >  5:182519/1‑148 (MQ=255)
                                                                                          cGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTAc                                                            <  6:283951/148‑1 (MQ=255)
                                                                                           gACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGc                                                                             <  8:36150/130‑1 (MQ=255)
                                                                                            aCAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGc                                                         <  2:338813/149‑1 (MQ=255)
                                                                                              aaaaTGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCaa                                                       <  7:219388/149‑1 (MQ=255)
                                                                                                      ctactCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAg                                               >  3:253130/1‑149 (MQ=255)
                                                                                                         ctCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTTATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAgt                                              <  5:30485/147‑1 (MQ=255)
                                                                                                         ctCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTa                                                                                                                    <  4:63625/77‑1 (MQ=255)
                                                                                                         ctCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTa                                                                                                                    >  3:63625/1‑77 (MQ=255)
                                                                                                         ctCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGc                                                     >  5:138756/1‑140 (MQ=255)
                                                                                                         ctCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGc                                                     <  6:138756/140‑1 (MQ=255)
                                                                                                          tCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTa                                           >  2:374456/1‑149 (MQ=255)
                                                                                                          tCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTa                                           >  5:48702/1‑149 (MQ=255)
                                                                                                            ccACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAcg                                         >  6:329145/1‑149 (MQ=255)
                                                                                                            ccACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAcg                                         >  1:39553/1‑149 (MQ=255)
                                                                                                            ccACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAcg                                         >  7:255771/1‑149 (MQ=255)
                                                                                                               cTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCgg                                      <  2:93909/149‑1 (MQ=255)
                                                                                                                  cGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGtg                                   >  2:186358/1‑149 (MQ=255)
                                                                                                                    gTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGtgct                                 >  7:239029/1‑149 (MQ=255)
                                                                                                                         ccAGTCGGTATCAAGTAAATCAACTTCGAGAAGTTTATCGCCAAGGCCCCTGAGAGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAATAACCCTCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCa                            <  4:253130/149‑1 (MQ=255)
                                                                                                                         ccAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCa                            <  5:329145/149‑1 (MQ=255)
                                                                                                                          cAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCAc                           <  3:260271/149‑1 (MQ=255)
                                                                                                                          cAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCAc                           >  2:357329/1‑149 (MQ=255)
                                                                                                                          cAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACACGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGCTGCTGCTCa                            >  8:38808/1‑148 (MQ=255)
                                                                                                                           aGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCAcc                          >  6:109545/1‑149 (MQ=255)
                                                                                                                            ggCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCAccc                         <  8:255771/149‑1 (MQ=255)
                                                                                                                               ggTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCcaca                      <  2:39553/149‑1 (MQ=255)
                                                                                                                                     aaGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGCACCCGCTGCGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCCGGTGCTGCTCACGCAAACACTg                 >  2:34262/1‑148 (MQ=255)
                                                                                                                                     aaGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAgt                                              >  2:105570/1‑119 (MQ=255)
                                                                                                                                       ggAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTc              <  5:109545/149‑1 (MQ=255)
                                                                                                                                         aaaTCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTca             <  7:301680/148‑1 (MQ=255)
                                                                                                                                          gaTCAACTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcag            >  1:293537/2‑148 (MQ=255)
                                                                                                                                              aaCTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcagca          >  2:150497/1‑146 (MQ=255)
                                                                                                                                              aaCTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcagca          <  1:150497/146‑1 (MQ=255)
                                                                                                                                               aCTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcagc           <  3:151499/144‑1 (MQ=255)
                                                                                                                                               aCTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcagc           >  4:151499/1‑144 (MQ=255)
                                                                                                                                               aCTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCa                            <  2:178760/127‑1 (MQ=255)
                                                                                                                                               aCTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCa                            >  1:178760/1‑127 (MQ=255)
                                                                                                                                                 ttCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAGCCTc     >  8:108249/1‑148 (MQ=255)
                                                                                                                                                   cGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAc                                          <  7:197766/109‑1 (MQ=255)
                                                                                                                                                   cGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAc                                          >  8:197766/1‑109 (MQ=255)
                                                                                                                                                    gagaAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAGCCTCAAg  >  1:102852/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAACCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAGCCTCAAG  >  NC_002947/1503318‑1503612

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: