Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,482,761 T→C 6.9% A229A (GCT→GCC dusB → tRNA dihydrouridine synthase DusB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,482,7610TC 6.9% 205.7 / 2.6 71A229A (GCT→GCCdusBtRNA dihydrouridine synthase DusB
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base C (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGGCGAAGCCGAGTACGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATGGGC  >  NC_002947/5482614‑5482908
                                                                                                                                                   |                                                                                                                                                   
cacCGGCGAAGCCGAGTACGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCt                                                                                                                                                     <  1:242420/148‑1 (MQ=255)
cacCGGCGAAGCCGAGTACGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCt                                                                                                                                                     <  2:229342/148‑1 (MQ=255)
cacCGGCGAAGCCGAGTACGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTg                                                                                                                                                    >  7:198086/1‑149 (MQ=255)
       gAAGCCGAGTACGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAggg                                                                                                                                             >  5:41805/1‑149 (MQ=255)
                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                                                                                                                                  <  4:76897/143‑1 (MQ=255)
                        cATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAGGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGa                                                                                                                                  >  3:76897/1‑143 (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                                                                                                                                        >  3:262077/1‑95 (MQ=255)
                                                                  cAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGggcggc                                                                                                                                        <  4:262077/95‑1 (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                                                                    >  7:200658/1‑148 (MQ=255)
                                                                                 cTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGAc                                                                    >  5:130329/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                                                           <  5:150409/149‑1 (MQ=255)
                                                                                         aaaaGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGaa                                                           <  6:163531/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           >  1:329669/1‑149 (MQ=255)
                                                                                                                                         ttGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGAt           <  4:50971/149‑1 (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)
                                                                                                                                                   |                                                                                                                                                   
CACCGGCGAAGCCGAGTACGACACCATCGGTGCCATCAAGCAGGCGGTGTCAATCCCGGTTTTTGCCAACGGCGATATCACCTCGCCAGAAAAGGCCCGGGCGGTGCTGGAAACCACCGGGGTCGATGGCCTGTTGATTGGCCGGGCTGCCCAGGGGCGGCCATGGATCTTTCGCGAGATCGAGCATTACCTGCGCACTGGCGAACATCTGCCAGCGCCGCAACTGGACGAAGTGGAACGCATCCTGCTGGAGCATCTGGCCGCGCTGCATGCCTTCTATGGCGATGTGATGGGC  >  NC_002947/5482614‑5482908

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: