Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,488,991 A→G 8.2% intergenic (+344/+227) PP_3981 → / ← PP_3982 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,488,9910AG 8.2% 201.5 / 6.7 73intergenic (+344/+227)PP_3981/PP_3982transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base G (3/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCACCACCGC  >  NC_002947/4488848‑4489137
                                                                                                                                               |                                                                                                                                                  
tGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTcc                                                                                                                                               <  1:225643/149‑1 (MQ=255)
         cTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCatat                                                                                                                                      <  2:79592/149‑1 (MQ=255)
         atCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAATATCCACCGATCGCCAAGATCAATCATCGCGCCCCGCGTGGAGCTGAAGTATAAAGGGCTTATCTGAGGCCTTTCTTGCACACGGGCGGGAACCTCCCGCCCatat                                                                                                                                      <  8:227358/148‑1 (MQ=255)
          tCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACATGGGCGGGAGGTTCCCGccc                                                                                                                                          <  4:121419/144‑1 (MQ=255)
          tCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACATGGGCGGGAGGTTCCCGccc                                                                                                                                          >  3:121419/1‑144 (MQ=255)
          tCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATg                                                                                                                                     >  8:2943/1‑149 (MQ=255)
          tCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATg                                                                                                                                     >  3:33346/1‑149 (MQ=255)
            gACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATgcg                                                                                                                                   >  2:295008/1‑149 (MQ=255)
            gACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATgc                                                                                                                                    <  5:236585/148‑1 (MQ=255)
            gACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATgc                                                                                                                                    <  7:2943/148‑1 (MQ=255)
              cATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGcgcg                                                                                                                                 >  8:360109/1‑149 (MQ=255)
                ttGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTa                                                                                                                               >  1:113715/1‑149 (MQ=255)
                    gcCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACtt                                                                                                                            <  3:19587/148‑1 (MQ=255)
                     ccATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCa                                                                                                                          <  8:202836/149‑1 (MQ=255)
                            gacTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCg                                                                                                                                             <  1:62886/123‑1 (MQ=255)
                            gacTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCg                                                                                                                                             >  2:62886/1‑123 (MQ=255)
                              cTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGGTCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGa                                                                                                                  >  7:104571/1‑148 (MQ=255)
                               tGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAGGTTCCCGcccg                                                                                                                                         <  2:102846/124‑2 (MQ=255)
                               tGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAGGTTCCCGcccg                                                                                                                                         >  1:102846/1‑123 (MQ=255)
                               tGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGAc                                                                                                                 >  1:361670/1‑148 (MQ=255)
                                     aTGCAGGTTGTCGAAAATCGACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCAc                                                                                                          <  4:33346/149‑1 (MQ=255)
                                       gCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTc                                                                                                        >  6:138663/1‑149 (MQ=255)
                                         aGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTcc                                                                                                       >  5:3533/1‑148 (MQ=255)
                                          ggTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCaa                                                                                                     <  2:345606/149‑1 (MQ=255)
                                                  aaaaTCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTca                                                                                              >  3:127734/1‑148 (MQ=255)
                                                     aTCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGTGCGTACTTCATAGCGCGACAATc                                                                                                             >  7:201208/1‑130 (MQ=255)
                                                     aTCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGTGCGTACTTCATAGCGCGACAATc                                                                                                             <  8:201208/130‑1 (MQ=255)
                                                     aTCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACAtc                                                                                          >  6:255083/1‑149 (MQ=255)
                                                          ccGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACt                                                                                     <  1:295008/149‑1 (MQ=255)
                                                                  cAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGaa                                                                             >  4:13959/1‑149 (MQ=255)
                                                                  cAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGaa                                                                             >  5:348188/1‑149 (MQ=255)
                                                                  cAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGaa                                                                             >  6:57547/1‑149 (MQ=255)
                                                                      aTCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTgg                                                                         >  2:286165/1‑149 (MQ=255)
                                                                      aTCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTgg                                                                         >  8:124554/1‑149 (MQ=255)
                                                                       tCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCACGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACAACACACTTGGATGAAcacgg                                                                        >  8:58360/1‑145 (MQ=255)
                                                                         aaGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGc                                                                       <  7:124554/148‑1 (MQ=255)
                                                                           gCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCt                                                                     >  7:317978/1‑148 (MQ=255)
                                                                              tCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTaaaa                                                                 >  2:287546/1‑149 (MQ=255)
                                                                              tCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAAACCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTaaaa                                                                 <  1:286165/149‑1 (MQ=255)
                                                                                gcgcCCAGCGTGGAGCTGAAGGATAAAGGGCTTAGCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTaaaa                                                                 <  6:348188/147‑1 (MQ=255)
                                                                                 cgcCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCATTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAActc                                                              <  1:287546/149‑1 (MQ=255)
                                                                                 cgcCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAActc                                                              <  5:57547/149‑1 (MQ=255)
                                                                                       gCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATa                                                        <  7:324429/149‑1 (MQ=255)
                                                                                        cGTGGAGCTGAAGGATAAAGGGCTTATCTGAGTCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATATCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAg                                                       <  7:360109/149‑1 (MQ=255)
                                                                                                tGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTcac                                                                                             >  1:65423/1‑103 (MQ=255)
                                                                                                tGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTcac                                                                                             <  2:65423/103‑1 (MQ=255)
                                                                                                  aaGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCt                                              <  3:13959/148‑1 (MQ=255)
                                                                                                     gATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACACTCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACg                                          >  4:285070/1‑149 (MQ=255)
                                                                                                     gATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACg                                          >  3:9089/1‑149 (MQ=255)
                                                                                                        aaaGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCt                                       <  5:138663/149‑1 (MQ=255)
                                                                                                         aaGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACCCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCt                                       <  8:104571/148‑1 (MQ=255)
                                                                                                          aGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGc                                     <  2:113715/149‑1 (MQ=255)
                                                                                                          aGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGc                                     >  7:185802/1‑149 (MQ=255)
                                                                                                               ttATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGAGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCcaca                                 >  7:204977/1‑148 (MQ=255)
                                                                                                                 aTCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACAt                                <  5:255083/147‑1 (MQ=255)
                                                                                                                       ggCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATAcc                         >  3:67648/1‑148 (MQ=255)
                                                                                                                         ccTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAg                      <  7:58360/149‑1 (MQ=255)
                                                                                                                          cttgcttgCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAAt                            <  3:359477/142‑1 (MQ=255)
                                                                                                                          cttgcttgCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAAt                            >  4:359477/1‑142 (MQ=255)
                                                                                                                                   atatgggcgggAGGTTCCCGCCCGTATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGtt            >  1:134570/5‑149 (MQ=255)
                                                                                                                                   acacGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGtt            >  6:97559/1‑149 (MQ=255)
                                                                                                                                     acGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGttt           >  6:53883/1‑148 (MQ=255)
                                                                                                                                     acGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTc          <  4:67648/149‑1 (MQ=255)
                                                                                                                                      tgggcgggAGGTTCCCGCCCGTATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTc          <  2:134570/147‑1 (MQ=255)
                                                                                                                                      cgggcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTca         <  1:74663/149‑1 (MQ=255)
                                                                                                                                       gggcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTcac        <  5:97559/149‑1 (MQ=255)
                                                                                                                                         gcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTcacca      >  7:272958/1‑149 (MQ=255)
                                                                                                                                         gagggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTc          <  3:44154/143‑1 (MQ=255)
                                                                                                                                         gagggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTc          >  4:44154/3‑145 (MQ=255)
                                                                                                                                          cgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTcaccac     <  5:53883/149‑1 (MQ=255)
                                                                                                                                           gggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGttt           <  2:401240/142‑1 (MQ=255)
                                                                                                                                           gggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGttt           >  1:401240/1‑142 (MQ=255)
                                                                                                                                             gAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCACCACCGc  >  2:337709/1‑149 (MQ=255)
                                                                                                                                              aaCCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCACCACCGc  >  8:50039/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCACCACCGC  >  NC_002947/4488848‑4489137

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: