Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,858,701 A→T 7.8% intergenic (‑207/+18) PP_mr50 ← / ← ligA Bacteria_large_SRP ncRNA/DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,858,7010AT 7.8% 250.5 / 7.8 86intergenic (‑207/+18)PP_mr50/ligABacteria_large_SRP ncRNA/DNA ligase
Reads supporting (aligned to +/- strand):  ref base A (36/42);  new base T (3/4);  total (40/46)
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.09e-01

GATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATGGGTCTT  >  NC_002947/4858556‑4858847
                                                                                                                                                  |                                                                                                                                                   
gATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAg                                                                                                                                                  <  4:356637/149‑1 (MQ=255)
  tACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAggg                                                                                                                                                <  3:85792/149‑1 (MQ=255)
   aCAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAgggg                                                                                                                                               >  7:77486/1‑149 (MQ=255)
    cAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAggggg                                                                                                                                              <  4:303647/149‑1 (MQ=255)
    cAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAgggg                                                                                                                                               <  4:381891/148‑1 (MQ=255)
       ccACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGCCCGCGCCGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCCCAGGGGGGCa                                                                                                                                           >  6:422186/1‑149 (MQ=255)
          cTGCGGCAGAAAACGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCTGTGGGGGGCaact                                                                                                                                        <  5:131396/149‑4 (MQ=39)
          cTGCGGCAGAAAAAGTCCTCCTTTCTCAATGACCTGTCAGCGTCGAGATGGCCGGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGGTGCCCCCCTGTGGGGGGCaact                                                                                                                                        <  8:147884/149‑4 (MQ=37)
          cTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAgggg                                                                                                                                               <  1:88530/142‑1 (MQ=255)
          cTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAgggg                                                                                                                                               >  2:88530/1‑142 (MQ=255)
          cTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACgg                                                                                                                                        <  6:161000/149‑1 (MQ=255)
              ggCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGc                                                                                                                                    <  3:146147/149‑1 (MQ=255)
                 aGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTa                                                                                                                                <  4:121292/149‑1 (MQ=255)
                  gAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGTACGGATGCCTAg                                                                                                                               <  5:105199/149‑1 (MQ=255)
                  gAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGcc                                                                                                                                   >  1:193623/1‑146 (MQ=255)
                   aaaaaGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACACGCTCCTACGGGCTACCGCGCATGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGCGGGGCACGGATGCCTAGa                                                                                                                              >  5:38105/1‑149 (MQ=255)
                     aaaGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCACTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGAc                                                                                                                             >  1:308604/1‑149 (MQ=255)
                        gTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCt                                                                                                                         <  3:152860/149‑1 (MQ=255)
                         tcatcaTTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTa                                                                                                                        >  1:180237/1‑149 (MQ=255)
                           atcatTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATg                                                                                                                      >  8:300832/1‑149 (MQ=255)
                               tttCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCt                                                                                                                  >  5:42128/1‑149 (MQ=255)
                                 tctcAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGt                                                                                                                >  3:46118/1‑149 (MQ=255)
                                   tcAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTc                                                                                                              >  2:180507/1‑149 (MQ=255)
                                    cAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCg                                                                                                             >  2:31205/1‑149 (MQ=255)
                                     aaTGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAAGCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCg                                                                                                             <  4:361357/148‑1 (MQ=255)
                                     aaTGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCg                                                                                                             <  4:355347/148‑1 (MQ=255)
                                     aaTGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCg                                                                                                             <  4:361944/148‑1 (MQ=255)
                                      aTGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGc                                                                                                           <  2:180237/149‑1 (MQ=255)
                                         aCTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCa                                                                                                         <  4:46118/148‑1 (MQ=255)
                                         aCTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAg                                                                                                        <  8:87730/149‑1 (MQ=255)
                                          gttGTCTGGGTCGAGTTGGCCTGTACGGGCGTCTTCGCGGTTAAACCGGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAg                                                                                                        <  6:42128/147‑1 (MQ=39)
                                                aGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAAtgt                                                                                                  <  4:133559/148‑1 (MQ=255)
                                                aGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAAtgt                                                                                                  <  2:308604/148‑1 (MQ=255)
                                                  ggTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGcc                                                                                                                                   <  3:224245/114‑1 (MQ=255)
                                                  ggTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGcc                                                                                                                                   >  4:224245/1‑114 (MQ=255)
                                                        gaTGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAgg                                                                                          >  4:273296/1‑148 (MQ=255)
                                                                 gACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTc                                                                                <  4:129557/149‑1 (MQ=255)
                                                                       ccTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTAt                                                                                                                       <  3:72540/104‑1 (MQ=255)
                                                                       ccTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTAt                                                                                                                       >  4:72540/1‑104 (MQ=255)
                                                                       ccTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCAc                                                                          >  5:277998/1‑149 (MQ=255)
                                                                            ttCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGAcc                                                                                                                            <  4:61539/95‑1 (MQ=255)
                                                                            ttCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGAcc                                                                                                                            >  3:61539/1‑95 (MQ=255)
                                                                            ttCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAAtgt                                                                                                  <  2:326849/121‑1 (MQ=255)
                                                                            ttCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAAtgt                                                                                                  >  1:326849/1‑121 (MQ=255)
                                                                               gcgGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTcaca                                                                   <  5:26688/149‑1 (MQ=255)
                                                                               gcgGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTcac                                                                    <  6:203844/148‑1 (MQ=255)
                                                                                       aCCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTc                                                           >  5:122733/1‑149 (MQ=255)
                                                                                        cccGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCg                                                          >  2:236193/1‑149 (MQ=255)
                                                                                         ccGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCgg                                                         <  1:242415/149‑1 (MQ=255)
                                                                                              ccTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGgccgcc                                                    >  4:8961/1‑149 (MQ=255)
                                                                                               cTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCt                                                   >  2:264361/1‑149 (MQ=255)
                                                                                                tACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCtt                                                  >  5:212645/1‑149 (MQ=255)
                                                                                                 aCGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTg                                                 >  6:202600/1‑149 (MQ=255)
                                                                                                  cGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGATGACAGTTGCCCCCCTGAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCTGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCGTTg                                                 <  6:38105/148‑1 (MQ=255)
                                                                                                       gACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCa                                                                                                         <  3:357021/87‑1 (MQ=255)
                                                                                                       gACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCa                                                                                                         >  4:357021/1‑87 (MQ=255)
                                                                                                        aCCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCt                                          >  4:37009/1‑149 (MQ=255)
                                                                                                        aCCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCt                                          <  6:122733/149‑1 (MQ=255)
                                                                                                           gcgcAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCg                                       <  1:264361/149‑1 (MQ=255)
                                                                                                                        aaGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTGAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAg                           >  1:297166/1‑149 (MQ=255)
                                                                                                                        aaGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                          <  7:321794/149‑1 (MQ=255)
                                                                                                                        aaGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                          <  2:297166/149‑1 (MQ=255)
                                                                                                                                 gaagaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGccac                 >  5:329212/1‑149 (MQ=255)
                                                                                                                                 gaagaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGcca                  <  1:180507/148‑1 (MQ=255)
                                                                                                                                  aagaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAg                                     <  8:349781/128‑1 (MQ=255)
                                                                                                                                  aagaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAg                                     >  7:349781/1‑128 (MQ=255)
                                                                                                                                  aagaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCCGGCCCAGccgc                 >  4:304313/1‑146 (MQ=255)
                                                                                                                                  aagaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                          >  4:194503/1‑139 (MQ=255)
                                                                                                                                  aagaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                          <  3:194503/139‑1 (MQ=255)
                                                                                                                                   agaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGccacc                <  2:193623/148‑1 (MQ=255)
                                                                                                                                       ccgtGCCCCCCTGTGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGc            >  5:219078/4‑148 (MQ=255)
                                                                                                                                        cgtGCCCCCCTGTGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCaa                                                              <  1:255525/95‑1 (MQ=255)
                                                                                                                                        cgtGCCCCCCTGTGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCaa                                                              >  2:255525/3‑97 (MQ=255)
                                                                                                                                         tttccccccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAAt         <  8:77486/146‑1 (MQ=255)
                                                                                                                                         gtGCCCCCCTGTGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAAt         >  4:293878/2‑149 (MQ=255)
                                                                                                                                           gCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATg        <  3:37009/148‑1 (MQ=255)
                                                                                                                                            cggcccACAGGGGGGCACGGATGCATAGTCCGCTATGTCCTGTTCGGCCAGGAATGTGACATAGGCATCTTCGTCCAGGTCCTTGACTGCGATCTGGGCCTGCTTGGGCAGCTTGGTGCCTGCGCCAGGGGCAGGGACGTCGCAATgg       <  2:195471/145‑1 (MQ=255)
                                                                                                                                            ccccccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCa                      <  7:348849/133‑1 (MQ=255)
                                                                                                                                            ccccccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCa                      >  8:348849/1‑133 (MQ=255)
                                                                                                                                            ccccccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATggc      >  5:75317/1‑148 (MQ=255)
                                                                                                                                            ccccccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATgg       <  3:273296/148‑1 (MQ=255)
                                                                                                                                             cccccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTc                                                                                <  8:116571/74‑1 (MQ=255)
                                                                                                                                             cccccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTc                                                                                >  7:116571/1‑74 (MQ=255)
                                                                                                                                             cccccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCCCCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATGGGt     >  3:224114/1‑149 (MQ=255)
                                                                                                                                               cccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATGGGTCt   >  8:341805/1‑149 (MQ=255)
                                                                                                                                                ccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATGGGTCtt  >  4:30497/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATGGGTCTT  >  NC_002947/4858556‑4858847

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: