Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,058 A→G 6.0% T75A (ACG→GCG)  PP_3774 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,0580AG 6.0% 185.7 / 2.4 66T75A (ACG→GCG) PP_3774hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base G (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

GCTGAACGGTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCAC  >  NC_002947/4300910‑4301206
                                                                                                                                                    |                                                                                                                                                    
gCTGAACGGTGCTGAAGCGATTGACCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGa                                                                                                                                                      >  2:77940/1‑149 (MQ=255)
  tGAACGGTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACACAAGACTGGACg                                                                                                                                                    >  7:127590/1‑149 (MQ=255)
    aaCGGTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACg                                                                                                                                                    <  7:2923/147‑1 (MQ=255)
       ggTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGc                                                                                                                                               <  4:345241/149‑1 (MQ=255)
         tGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGc                                                                                                                                               >  2:38692/1‑147 (MQ=255)
            tGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCGGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTg                                                                                                                                                 >  4:276750/1‑142 (MQ=255)
            tGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCGGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTg                                                                                                                                                 <  3:276750/142‑1 (MQ=255)
            tGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGATTGGACGGAAACGCAAGACTGGGCGCTGGCCAgcg                                                                                                                                          >  4:23787/1‑149 (MQ=255)
              aaGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGcc                                                                                                                                        >  3:171344/1‑149 (MQ=255)
               aGCGATTGACCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGccc                                                                                                                                       <  1:77940/149‑1 (MQ=255)
                gCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCa                                                                                                                                      <  1:38692/149‑1 (MQ=255)
                     tGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGAt                                                                                                                                  >  7:88757/1‑148 (MQ=255)
                      gCCCGCTTTCCCAGCATCGTGACGGTGCTCTGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAgcgtc                                                                                                                                        >  1:30916/1‑139 (MQ=255)
                      gCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAgcgtc                                                                                                                                        <  2:30916/141‑3 (MQ=255)
                       cccGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGATTGGACGGAAACGCAAGACTGGGCGCTGGCCAGCGTCCAgtcttgc                                                                                                                               <  3:23787/149‑7 (MQ=255)
                        ccGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACt                                                                                                                              >  4:250484/1‑149 (MQ=255)
                                 cAGTATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGcc                                                                                                                     <  7:85297/149‑1 (MQ=255)
                                      tCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGAc                                                                                                                 >  5:1313/1‑148 (MQ=255)
                                         tGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAAt                                                                                                              >  2:201440/1‑148 (MQ=255)
                                           aCGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGAc                                                                                                           <  8:243601/149‑1 (MQ=255)
                                              gTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACACAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCt                                                                                                        <  8:127590/149‑1 (MQ=255)
                                                    aGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAAt                                                                                                              <  3:17507/137‑1 (MQ=255)
                                                    aGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAAt                                                                                                              >  4:17507/1‑137 (MQ=255)
                                                             tCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTc                                                                                         >  4:149270/1‑149 (MQ=255)
                                                                 tCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCa                                                                                                                    <  5:341191/118‑1 (MQ=255)
                                                                 tCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCa                                                                                                                    >  6:341191/1‑118 (MQ=255)
                                                                    aGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGgcg                                                                                     >  3:106460/1‑146 (MQ=255)
                                                                     gTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGa                                                                                 <  8:88757/149‑1 (MQ=255)
                                                                      ttCTCAGCGCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCCGGATCAGTGGAATCGCCATACAATTACGCTGACCAAGGTAAAATCGGCGCAGaa                                                                                >  1:193508/1‑149 (MQ=255)
                                                                         tcAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggt                                                                             >  6:227655/1‑149 (MQ=255)
                                                                         tcAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggt                                                                             >  6:87585/1‑149 (MQ=255)
                                                                          ccgCCCCACCCCCAGGGGCGTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtg                                                                            <  3:280758/147‑1 (MQ=255)
                                                                          cAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtg                                                                            >  7:21102/1‑149 (MQ=255)
                                                                           aGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtgg                                                                           <  5:87585/149‑1 (MQ=255)
                                                                           aGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtg                                                                            <  3:149270/148‑1 (MQ=255)
                                                                            gCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtggt                                                                          >  2:91759/1‑149 (MQ=255)
                                                                             ccccAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtggt                                                                          >  3:52753/1‑148 (MQ=255)
                                                                               ccaaccaaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCAt                                                                       >  5:331280/1‑149 (MQ=255)
                                                                               ccaaccaaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGCCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCAt                                                                       >  7:320617/1‑149 (MQ=255)
                                                                                caaccaaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATc                                                                      <  1:259120/149‑1 (MQ=255)
                                                                                 aaccaaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCg                                                                     <  4:171344/149‑1 (MQ=255)
                                                                                     aatggggCTTCGAGCAAAAGCACCCACAGGACAGGCTGCGAATGGAAGGAAACGCAATACTGGAAGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCACAGAAGGTGGTCATCGcc                                                                   <  4:277924/144‑1 (MQ=255)
                                                                                     aaaGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCACAGAAGGTGGTCATCGcc                                                                   >  3:277924/1‑147 (MQ=255)
                                                                                     aaaGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATAACTGGAATAGCCAGACAATGACGCTGACCAAGGTCAACTCGGCACAGAAGGTGGTCATCGcc                                                                   >  3:282688/1‑147 (MQ=255)
                                                                                                 aGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTAcg                                                     >  5:142169/1‑149 (MQ=255)
                                                                                                      aaGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTCGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCggg                                                 >  8:244583/1‑148 (MQ=255)
                                                                                                            ccACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGacgacg                                          <  3:250484/149‑1 (MQ=255)
                                                                                                             cacaGGACACGCTGCGAATGGACGGAAACGGAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACgg                                         >  1:29310/1‑149 (MQ=255)
                                                                                                              acaGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCAc                                                          <  5:333704/131‑1 (MQ=255)
                                                                                                              acaGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCAc                                                          >  6:333704/1‑131 (MQ=255)
                                                                                                              acaGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGc                                        >  8:120184/1‑149 (MQ=255)
                                                                                                                aGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGccc                                      <  5:227655/149‑1 (MQ=255)
                                                                                                                  gACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAg                                                                                                 >  6:333107/1‑88 (MQ=255)
                                                                                                                  gACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAg                                                                                                 <  5:333107/88‑1 (MQ=255)
                                                                                                                    cacGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGacgacg                                          <  3:261830/141‑1 (MQ=255)
                                                                                                                    cacGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCccgccg                                  <  1:201440/149‑1 (MQ=255)
                                                                                                                     acGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCgg                                 <  6:142169/149‑1 (MQ=255)
                                                                                                                       gCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGt                                >  8:108756/1‑148 (MQ=255)
                                                                                                                              aTGGACTGAAACGCAAGACTTGACGCTGGCCAGCGCCCAGGATCACTTTAATCGCCAGACAATGACGCTGACCAAGTTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGgtg                         <  2:193508/148‑1 (MQ=255)
                                                                                                                                  aCGGAAACGCAAGACTGGACGCTGGCCAGCGCCCATGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTg                    <  6:1313/149‑1 (MQ=255)
                                                                                                                                    ggAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtg                                                                            >  2:46386/1‑91 (MQ=255)
                                                                                                                                    ggAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtg                                                                            <  1:46386/91‑1 (MQ=255)
                                                                                                                                          gCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAg                            >  5:327484/1‑133 (MQ=255)
                                                                                                                                          gCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAg                            <  6:327484/133‑1 (MQ=255)
                                                                                                                                          gCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAgg            <  7:244583/149‑1 (MQ=255)
                                                                                                                                            aaGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAggg           <  2:330059/148‑1 (MQ=255)
                                                                                                                                            aaGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAggg           >  1:330059/1‑148 (MQ=255)
                                                                                                                                                    aCGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCAc  <  7:120184/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTGAACGGTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCAC  >  NC_002947/4300910‑4301206

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