Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,492,699 G→A 8.3% E194E (GAG→GAA)  PP_3983 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,492,6990GA 8.3% 198.7 / 7.2 72E194E (GAG→GAA) PP_3983hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/35);  new base A (3/3);  total (34/38)
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.90e-01

GCATGCAGATGACGATCTGCTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGATCAGTACCACCATGCTCGAACCTTCGCTTGG  >  NC_002947/4492555‑4492843
                                                                                                                                                |                                                                                                                                                
gCATGCAGATGACGATCTGCTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCtt                                                                                                                                               <  8:149599/148‑1 (MQ=255)
gCATGCAGATGACGATCTGCTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTg                                                                                                                                              <  8:34767/149‑1 (MQ=255)
 cATGCAGATGTCGTTCTGCTGATTCTATTGCTGGCAATCAAAAGCTCTTCTAAGAGCAATCCAAGCTACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCGTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTg                                                                                                                                              <  7:127993/148‑1 (MQ=255)
  aTGCAGATGACGATCTGCTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCGATCTCTACGACATCCTAGAAGCCCTACAATGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTgg                                                                                                                                             >  4:174416/1‑148 (MQ=255)
  aTGCAGATGACGATCTGCTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTggg                                                                                                                                            <  7:306766/149‑1 (MQ=255)
      aGATGACGATCTGCTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATc                                                                                                                                         >  4:49961/1‑148 (MQ=255)
          gACGATCTGCTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCgaga                                                                                                                                     >  7:247141/1‑148 (MQ=255)
           aCGATCTGCTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGaaa                                                                                                                                   <  8:241790/149‑1 (MQ=255)
                   cTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCg                                                                                                                                        <  1:140158/136‑1 (MQ=255)
                   cTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGt                                                                                                                           >  3:238656/1‑149 (MQ=255)
                   cTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGATCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCg                                                                                                                                        >  2:140158/1‑136 (MQ=255)
                   atCAGTCGATTGCTGTCAAGAAAAAGAGATTCTAATATCAATCCAAGCGACGATTTCACCCATCTATACGAAATCCTAGAAGCCCTACAAATCAAGTTCGTGAAACCTGCTACATTGATCAACGAGCTTGGTATCGAGAAAGCGAAg                                                                                                                             <  4:100141/146‑1 (MQ=255)
                                tGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCAt                                                                                                              >  6:170583/1‑149 (MQ=255)
                                tGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCCAGAAAGCGAAGGTACGCAAATTCCAt                                                                                                              >  3:6383/1‑149 (MQ=255)
                                 ggCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCAt                                                                                                              >  3:331035/1‑148 (MQ=255)
                                    aagcaaAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAg                                                                                                          >  7:130643/1‑149 (MQ=255)
                                       caaAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCTCCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCt                                                                                                       >  1:234793/1‑149 (MQ=255)
                                       caaAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCt                                                                                                       >  8:206891/1‑149 (MQ=255)
                                        aaaaGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTg                                                                                                      >  1:126168/1‑149 (MQ=255)
                                        aaaaGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTg                                                                                                      <  7:229273/149‑1 (MQ=255)
                                           aGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCaa                                                                                                   <  4:238656/149‑1 (MQ=255)
                                                ttCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGAcc                                                                                              >  2:281220/1‑149 (MQ=255)
                                                  cTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGAccc                                                                                             <  3:49961/148‑1 (MQ=255)
                                                     agagCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCa                                                                                            >  1:169545/1‑146 (MQ=255)
                                                     agagCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCa                                                                                            <  2:169545/146‑1 (MQ=255)
                                                     agagCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATc                                                                                          <  2:2604/148‑1 (MQ=255)
                                                       agCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAAc                                                                                       >  4:147616/1‑149 (MQ=255)
                                                              cAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAgc                                                                                >  5:71686/1‑149 (MQ=255)
                                                                   gacgaGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCCGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCAACAACTGCCAGCCGCCa                                                                           >  1:285116/1‑149 (MQ=255)
                                                                       aGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAATGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGa                                                                       <  1:90745/149‑1 (MQ=255)
                                                                       aGTTCACCAATCTCTACGAAATCGTAGAAGCCCTACAATGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGGAATGACCCATCAACTGCCAGCCGCCAGCGa                                                                       <  2:234793/149‑1 (MQ=255)
                                                                        gTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAg                                                                      >  1:201528/1‑149 (MQ=255)
                                                                         ttCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAgg                                                                     <  2:126168/149‑1 (MQ=255)
                                                                            aCCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAACTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCATTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCt                                                                  >  3:285889/1‑149 (MQ=255)
                                                                            aCCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAACTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCt                                                                  >  5:11604/1‑149 (MQ=255)
                                                                            aCCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAACTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTAAGTCTGCAATGACCAATCAACTGCCAGCCGCCACCGAGGGCt                                                                  >  5:6633/1‑149 (MQ=255)
                                                                              ccATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTgg                                                                <  5:170583/149‑1 (MQ=255)
                                                                                     tACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCaaaa                                                           <  7:206891/147‑1 (MQ=255)
                                                                                     tACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAgg                                                         >  7:165704/1‑149 (MQ=255)
                                                                                       cGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGcc                                                       >  5:91815/1‑149 (MQ=255)
                                                                                           aTCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCAGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCACCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGtt                                                   <  2:285116/149‑1 (MQ=255)
                                                                                              cTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTg                                                                                                      <  4:227678/95‑1 (MQ=255)
                                                                                              cTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTg                                                                                                      >  3:227678/1‑95 (MQ=255)
                                                                                              cTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTg                                                  >  5:225645/1‑147 (MQ=255)
                                                                                              cTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTg                                                  <  6:225645/147‑1 (MQ=255)
                                                                                              cTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGa                                                 >  5:163594/1‑148 (MQ=255)
                                                                                              cTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGa                                                 >  1:345538/1‑148 (MQ=255)
                                                                                                 gAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGaaa                                               <  6:71686/147‑1 (MQ=255)
                                                                                                   aGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAagag                                            <  3:128680/148‑1 (MQ=255)
                                                                                                    gCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGcc                                          <  2:201528/149‑1 (MQ=255)
                                                                                                       cTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGaa                                                                                                                              >  4:196805/1‑62 (MQ=255)
                                                                                                       cTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGaa                                                                                                                              <  3:196805/62‑1 (MQ=255)
                                                                                                        tACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATACAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAgc                                      >  3:46305/1‑149 (MQ=255)
                                                                                                            aaGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTg                                  <  6:163594/149‑1 (MQ=255)
                                                                                                             cccAAGTTCGGGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGa                                 <  8:247141/147‑1 (MQ=255)
                                                                                                              gCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGAt                                <  1:281220/149‑1 (MQ=255)
                                                                                                              gCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGAt                                <  3:147616/149‑1 (MQ=255)
                                                                                                              gCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGAt                                <  8:130643/149‑1 (MQ=255)
                                                                                                              gCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGACAAAGCGAACGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGAt                                <  3:174416/149‑1 (MQ=255)
                                                                                                               cAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGCCCCATCAACTGCCAGCCGCCACCGCGTGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGAt                                >  8:303556/1‑148 (MQ=255)
                                                                                                                        gAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTgg                                                                >  6:153026/1‑107 (MQ=255)
                                                                                                                             cGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGATCAGTACCACCATGCt                 <  4:331035/149‑1 (MQ=255)
                                                                                                                              ggCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGATCAGTACCACCATGCTc                <  6:91815/149‑1 (MQ=255)
                                                                                                                                           aaCGAACTTGGGATCGAGAAAGCGTAGTTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCTCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGATCAGTACCACCATGCTCGAACCTTCGCTTg   <  4:285889/149‑1 (MQ=255)
                                                                                                                                           aaCGAACTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGATCAGTACCACCATGCTCGAACCTTCGCTTg   <  6:11604/149‑1 (MQ=255)
                                                                                                                                           aaCGAACTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACACATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGATCAGTACCACCATGCTCGAACCTTCGCTTg   <  6:6633/149‑1 (MQ=255)
                                                                                                                                            aCGAGCTTGGGATCTAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGATCAGTa                           <  7:66970/124‑1 (MQ=255)
                                                                                                                                            aCGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGATCAGTa                           >  8:66970/1‑124 (MQ=255)
                                                                                                                                            aCGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGATCAGTACCACCATGCTCGAACCTTCGCTTgg  >  3:9930/1‑149 (MQ=255)
                                                                                                                                            aCGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAAATCGATTCAGTCTGCATTGACCCATCAACTGCCAGCGGCCACCGAGGGCTGGCAAAAGGCCCGTTGACAGAGCCCAGCGCTGATCAGTACCACCATGCTCGAACCTTCGCTTgg  <  5:109667/149‑1 (MQ=255)
                                                                                                                                            aCGAGCTTGCGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAgg                                                                     >  7:249283/1‑82 (MQ=255)
                                                                                                                                            aCGAGCTTGCGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAgg                                                                     <  8:249283/82‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCATGCAGATGACGATCTGCTCAGTCGATTGCTGGCAAGCAAAAGCGCTTCTAAGAGCAATCCAAGCGACGAGTTCACCCATCTCTACGAAATCCTAGAAGCCCTACAAAGCAAGTTCGTGAAACCGGCGACATTGATCAACGAGCTTGGGATCGAGAAAGCGAAGGTTCGCAAATTCCATTCAGTCTGCAATGACCCATCAACTGCCAGCCGCCACCGAGGGCTGGCAAAAGGCCCGTTGAAAGAGCCCAGCGCTGATCAGTACCACCATGCTCGAACCTTCGCTTGG  >  NC_002947/4492555‑4492843

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