Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,482,772 G→A 7.1% R233Q (CGG→CAG)  dusB → tRNA dihydrouridine synthase DusB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,482,7720GA 7.1% 203.1 / 3.8 70R233Q (CGG→CAG) dusBtRNA dihydrouridine synthase DusB
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base A (3/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

TACGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATGGGCGTACGTATCG  >  NC_002947/5482630‑5482918
                                                                                                                                              |                                                                                                                                                  
tACGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCCGCCCCTGGGCAGcccgg                                                                                                                                              >  1:329318/1‑146 (MQ=255)
 aCGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCCGCCCCTGGGCAGcccgg                                                                                                                                              <  2:329318/148‑4 (MQ=255)
  cGACACCATCGGTTCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATgg                                                                                                                                             <  2:172112/148‑1 (MQ=255)
  cGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATgg                                                                                                                                             <  8:187850/148‑1 (MQ=255)
   gACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGAt                                                                                                                                           >  6:301277/1‑149 (MQ=255)
      acCATCGGTGCCATCAAGCAGGCGGTGCCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCtt                                                                                                                                        <  5:148285/149‑1 (MQ=255)
        cATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAGGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGa                                                                                                                                            >  3:76897/1‑143 (MQ=255)
        cATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAGGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGa                                                                                                                                            <  4:76897/143‑1 (MQ=255)
        cATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTc                                                                                                                                      >  3:250000/1‑149 (MQ=255)
         aTCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTcg                                                                                                                                     >  5:195206/1‑149 (MQ=255)
                 cATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATc                                                                                                                              >  6:208300/1‑148 (MQ=255)
                  aTCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATc                                                                                                                              >  6:270900/1‑147 (MQ=255)
                  aTCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATc                                                                                                                              >  6:150409/1‑147 (MQ=255)
                  aTCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCCGCCCCTGGGCAGcccgg                                                                                                                                              <  8:240438/131‑4 (MQ=255)
                  aTCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCCGCCCCTGGGCAGcccgg                                                                                                                                              >  7:240438/1‑128 (MQ=255)
                   tCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAg                                                                                                                           <  3:17094/149‑1 (MQ=255)
                     aaGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGc                                                                                                                          <  5:301277/148‑1 (MQ=255)
                       gCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCAtt                                                                                                                       <  6:100267/149‑1 (MQ=255)
                         aGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTAc                                                                                                                     <  6:41805/149‑1 (MQ=255)
                          ggCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTAcc                                                                                                                    >  3:226655/1‑149 (MQ=255)
                           gCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCt                                                                                                                   >  5:174860/1‑149 (MQ=255)
                              gtgtCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTgcg                                                                                                                >  3:50971/1‑149 (MQ=255)
                               tgtCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTgcgc                                                                                                               <  3:20498/149‑1 (MQ=255)
                                 tCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAg                                                                                                                           >  2:311164/1‑135 (MQ=255)
                                 tCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAg                                                                                                                           <  1:311164/135‑1 (MQ=255)
                                       ccGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCg                                                                                                        >  5:271588/1‑148 (MQ=255)
                                       ccGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCg                                                                                                        >  5:254775/1‑148 (MQ=255)
                                              ttGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCt                                                                                                 <  6:271588/148‑1 (MQ=255)
                                              ttGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCt                                                                                                 <  4:250000/148‑1 (MQ=255)
                                                 ccAACGGCGATATCACCTCGCCTGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCAGGGATCTTTCGCGAGATCGAGCATTACCTGCGCCCTGGCGACCAGCTGccc                                                                                             >  7:50542/1‑148 (MQ=255)
                                                 ccAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGcc                                                                                              <  8:198086/148‑1 (MQ=255)
                                                 ccAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGc                                                                                               >  3:208827/1‑147 (MQ=255)
                                                  cAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGggcggc                                                                                                                                                  <  4:262077/95‑1 (MQ=255)
                                                  cAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGggcggc                                                                                                                                                  >  3:262077/1‑95 (MQ=255)
                                                  cAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAg                                                                                            >  3:132650/1‑149 (MQ=255)
                                                          atatCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCa                                                                                     >  6:48804/1‑148 (MQ=255)
                                                               aCCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGa                                                                               <  5:48804/149‑1 (MQ=255)
                                                               aCCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGa                                                                               <  4:213047/149‑1 (MQ=255)
                                                                 cTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGAc                                                                              >  5:130329/1‑148 (MQ=255)
                                                                 cTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGAc                                                                              >  7:200658/1‑148 (MQ=255)
                                                                     ccAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGt                                                                         <  6:195206/149‑1 (MQ=255)
                                                                     ccAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGt                                                                         <  5:208300/149‑1 (MQ=255)
                                                                      cAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTg                                                                        >  2:70122/1‑149 (MQ=255)
                                                                         aaaaGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGaa                                                                     <  6:163531/149‑1 (MQ=255)
                                                                         aaaaGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGaa                                                                     <  5:150409/149‑1 (MQ=255)
                                                                          aaaGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAAc                                                                    <  5:270900/149‑1 (MQ=255)
                                                                               cccGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATc                                                               <  4:132650/149‑1 (MQ=255)
                                                                                     cGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCctgctg                                                         <  6:254775/149‑1 (MQ=255)
                                                                                            ggAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATc                                                  <  1:70122/149‑1 (MQ=255)
                                                                                            ggAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATc                                                  >  3:278639/1‑149 (MQ=255)
                                                                                             gAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCt                                                 >  3:301753/1‑149 (MQ=255)
                                                                                             gAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCt                                                 >  7:199545/1‑149 (MQ=255)
                                                                                              aaaCCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTg                                                <  5:132703/149‑1 (MQ=255)
                                                                                                    ccGGGGTCGATGGCCTATTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCcgcg                                          >  2:309881/1‑149 (MQ=255)
                                                                                                        ggTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCTTCTGGCCGCGCTGc                                      <  4:226655/149‑1 (MQ=255)
                                                                                                        ggTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGc                                      <  8:199545/149‑1 (MQ=255)
                                                                                                            gATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGc                                  >  3:150086/1‑149 (MQ=255)
                                                                                                            gATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGc                                  <  4:301753/149‑1 (MQ=255)
                                                                                                                  cTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTa                            >  7:148376/1‑149 (MQ=255)
                                                                                                                   tGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTAt                           >  3:15447/1‑149 (MQ=255)
                                                                                                                     ttGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATg                          <  4:208827/148‑1 (MQ=255)
                                                                                                                       gATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGc                        <  8:200658/148‑1 (MQ=255)
                                                                                                                         ttGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGAt                     <  4:50971/149‑1 (MQ=255)
                                                                                                                         ttGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGAt                     >  1:329669/1‑149 (MQ=255)
                                                                                                                           ggCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGAtg                    <  6:174860/148‑1 (MQ=255)
                                                                                                                            ccatggCCGCCCCTGGGCAGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGAtgt                   >  1:197584/5‑148 (MQ=255)
                                                                                                                               gggCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATg               >  2:225876/1‑149 (MQ=255)
                                                                                                                                ggCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATgg              <  1:309881/149‑1 (MQ=255)
                                                                                                                                  cTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATGGGc            <  4:278639/149‑1 (MQ=255)
                                                                                                                                      ccAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATGGGcgtac        >  4:233437/1‑149 (MQ=255)
                                                                                                                                        aGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATGGGcgtacgt      <  4:150086/149‑1 (MQ=255)
                                                                                                                                         ggggCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATGGGcgtacgta     <  1:225876/149‑1 (MQ=255)
                                                                                                                                            gcggcCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATGGGCGTACGTATCg  <  8:148376/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TACGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATGGGCGTACGTATCG  >  NC_002947/5482630‑5482918

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: ‑