Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,475,631 A→T 5.5% intergenic (‑52/+12) PP_3964 ← / ← PP_3965 transposase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,475,6310AT 5.5% 295.6 / 6.5 91intergenic (‑52/+12)PP_3964/PP_3965transposase/transposase
Reads supporting (aligned to +/- strand):  ref base A (49/37);  new base T (2/3);  total (51/40)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

TCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCA  >  NC_002947/4475483‑4475779
                                                                                                                                                    |                                                                                                                                                     
tCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCaaa                                                                                                                                                      <  1:710329/150‑1 (MQ=255)
  aTCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATg                                                                                                                                                    <  1:106578/150‑1 (MQ=255)
  aTCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGCCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATg                                                                                                                                                    >  2:587816/1‑150 (MQ=255)
    cGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGc                                                                                                                                                  <  1:2220521/150‑1 (MQ=255)
        tCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGc                                                                                                                                              <  1:448250/150‑1 (MQ=255)
        tCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGc                                                                                                                                              >  2:1651146/1‑150 (MQ=255)
        tCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGc                                                                                                                                              >  2:448250/1‑150 (MQ=255)
        tCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGc                                                                                                                                              <  1:1651146/150‑1 (MQ=255)
         caaCCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCTATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCATTTGGCCggtc                                                                                                                                             <  1:970071/148‑3 (MQ=255)
           aCCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGtt                                                                                                                                           <  2:1390424/150‑1 (MQ=255)
             cTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTa                                                                                                                                         >  2:486690/1‑150 (MQ=255)
              tGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTAc                                                                                                                                        >  1:1537196/1‑150 (MQ=255)
                ggATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAg                                                                                                                                      <  1:486690/150‑1 (MQ=255)
                  attgcaattgcaACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCa                                                                                                                                    <  2:1537196/150‑1 (MQ=255)
                       aattgcaaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACATCATGCTa                                                                                                                               >  1:1657041/1‑150 (MQ=255)
                       aattgcaaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTa                                                                                                                               >  1:1348779/1‑150 (MQ=255)
                        attgcaaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTaa                                                                                                                              >  1:450170/1‑150 (MQ=255)
                             aaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCATTTGGCCggt                                                                                                                                              <  2:1758121/129‑2 (MQ=255)
                             aaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCATTTGGCCggt                                                                                                                                              >  1:1758121/1‑128 (MQ=255)
                             aaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATcaca                                                                                                                         >  1:315315/1‑150 (MQ=255)
                                 ccTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGc                                                                                                                     <  1:630436/150‑1 (MQ=255)
                                  cTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCATTTGGCCggtctac                                                                                                                                          >  1:2238379/1‑123 (MQ=17)
                                  cTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCATTTGGCCggtctac                                                                                                                                          <  2:2238379/128‑6 (MQ=17)
                                  cTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGc                                                                                                                                                  >  2:1491612/1‑120 (MQ=255)
                                  cTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGc                                                                                                                                                  <  1:1491612/120‑1 (MQ=255)
                                  cTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGcc                                                                                                                    <  1:1635681/150‑1 (MQ=255)
                                    gtgCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTg                                                                                                                  >  2:1817361/1‑150 (MQ=255)
                                    gtgCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTg                                                                                                                  <  2:222359/150‑1 (MQ=255)
                                       tcggcAAGGCACCAAGCTGTGCCGGCGTCATCTGGTCGAGGTTGGGCGGGAAAGTCGCGCCGCCGATTGTTCCAGAGCCGGCGCGAGGCGACGAAAAGGAGACCGGCCAAATGGTCTGCTTTTACAGCATGCTAATCACACCGCCTGCAc                                                                                                               <  1:92160/146‑1 (MQ=255)
                                        cgccAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCAcc                                                                                                              <  2:156985/150‑1 (MQ=255)
                                         gccACGGCACGAAGCTGTTCCGGGGTCATCTGGTCGGGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCAcca                                                                                                             >  2:92160/1‑149 (MQ=255)
                                            aaGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGAcg                                                                                                          >  2:301112/1‑150 (MQ=255)
                                                   gAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAg                                                                                                                                      <  2:1890809/115‑1 (MQ=255)
                                                   gAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAg                                                                                                                                      >  1:1890809/1‑115 (MQ=255)
                                                   gAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGa                                                                                                            >  1:2397615/1‑141 (MQ=255)
                                                   gAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGa                                                                                                            <  2:2397615/141‑1 (MQ=255)
                                                   gAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGcc                                                                                                   >  1:1915709/1‑150 (MQ=255)
                                                          ttCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAg                                                                                            >  2:120658/1‑150 (MQ=255)
                                                                ggTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAg                                                                                      >  1:1807315/1‑150 (MQ=255)
                                                                  tCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGtt                                                                                                                                           <  1:2355019/95‑1 (MQ=255)
                                                                  tCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGtt                                                                                                                                           >  2:2355019/1‑95 (MQ=255)
                                                                    aTCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCAcc                                                                                  <  2:1461367/150‑1 (MQ=255)
                                                                     tCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGTCCAGCACCa                                                                                 >  1:517494/1‑150 (MQ=255)
                                                                     tCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCa                                                                                 >  1:1630626/1‑150 (MQ=255)
                                                                     tCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCa                                                                                 >  2:1457543/1‑150 (MQ=255)
                                                                       tGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAAc                                                                               >  2:16409/1‑150 (MQ=255)
                                                                         gTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAAt                                                                                                                             <  1:1992948/102‑1 (MQ=255)
                                                                         gTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAAt                                                                                                                             >  2:1992948/1‑102 (MQ=255)
                                                                         gTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGc                                                                             >  2:210517/1‑150 (MQ=255)
                                                                         gTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGc                                                                             <  1:1333633/150‑1 (MQ=255)
                                                                           cGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCt                                                                           <  1:1594711/150‑1 (MQ=255)
                                                                            gAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAATTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCa                                                                                       <  1:1496260/137‑1 (MQ=255)
                                                                            gAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCa                                                                                       >  2:1496260/1‑137 (MQ=255)
                                                                            gAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTg                                                                          >  1:2249629/1‑150 (MQ=255)
                                                                            gAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTg                                                                          >  1:1351306/1‑150 (MQ=255)
                                                                              ggTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGa                                                                        <  1:1457543/150‑1 (MQ=255)
                                                                                    gCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGcc                                                                                         >  1:2295372/1‑127 (MQ=255)
                                                                                    gCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGcc                                                                                         <  2:2295372/127‑1 (MQ=255)
                                                                                    gCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAgttgtt                                                                 >  1:2070488/1‑150 (MQ=255)
                                                                                     cGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTc                                                                <  2:1630626/150‑1 (MQ=255)
                                                                                       agaAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGc                                                              <  1:210517/150‑1 (MQ=255)
                                                                                         aaaGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTg                                                            >  1:643014/1‑150 (MQ=255)
                                                                                         aaaGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTg                                                            >  2:2195314/1‑150 (MQ=255)
                                                                                               aTGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCa                                                                                                                                    <  2:2004568/73‑1 (MQ=255)
                                                                                               aTGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCa                                                                                                                                    >  1:2004568/1‑73 (MQ=255)
                                                                                                tGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTAGCTGTCGAGtt                                                     >  1:715565/1‑150 (MQ=255)
                                                                                                  ccgccgATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCg                                                   >  1:140201/1‑150 (MQ=255)
                                                                                                    gccgATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGAc                                                 <  2:315315/150‑1 (MQ=255)
                                                                                                    gccgATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGAc                                                 <  2:1807315/150‑1 (MQ=255)
                                                                                                       gATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCATCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTc                                              >  2:1805726/1‑150 (MQ=255)
                                                                                                        aTTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCg                                             >  2:459762/1‑150 (MQ=255)
                                                                                                        aTTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCg                                             >  2:2106867/1‑150 (MQ=255)
                                                                                                           gtgCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCATCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAg                                          <  1:1805726/150‑1 (MQ=255)
                                                                                                             ggcAGAGCTGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGcc                                        <  2:2067430/148‑1 (MQ=255)
                                                                                                              ccAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCg                                       >  2:684089/1‑150 (MQ=255)
                                                                                                               cAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCt                                                                           >  1:929382/1‑114 (MQ=255)
                                                                                                               cAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCt                                                                           <  2:929382/114‑1 (MQ=255)
                                                                                                               cAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGt                                      >  2:204202/1‑150 (MQ=255)
                                                                                                                       gcgcGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATg                              >  1:1849314/1‑150 (MQ=255)
                                                                                                                       gcgcGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATg                              >  1:1652236/1‑150 (MQ=255)
                                                                                                                          cgGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCg                           <  2:1348779/150‑1 (MQ=255)
                                                                                                                          cgGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGAACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCg                           <  1:1793283/150‑1 (MQ=255)
                                                                                                                              gCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAgtt                                                                    <  2:2228472/105‑1 (MQ=255)
                                                                                                                              gCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAgtt                                                                    >  1:2228472/1‑105 (MQ=255)
                                                                                                                              gCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGcc                       >  1:66471/1‑150 (MQ=255)
                                                                                                                                gacgacAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAg                     <  1:16409/150‑1 (MQ=255)
                                                                                                                                  agaacaGTAGACGCGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTg                   <  2:140201/145‑1 (MQ=255)
                                                                                                                                   gacAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTgg                  >  1:991066/1‑150 (MQ=255)
                                                                                                                                           gACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGcc          >  1:1814858/1‑150 (MQ=255)
                                                                                                                                            aCCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCt         >  2:2334310/1‑150 (MQ=255)
                                                                                                                                            aCCGGACAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCATGCACCGAAGGGCTGCCACGACATGGCCGGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCt         <  1:285865/150‑1 (MQ=17)
                                                                                                                                                 ccAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGAATGGACAATGCCGGAAGATCAGACAACCGCGCTGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCGCGGtt    >  1:284608/1‑150 (MQ=12)
                                                                                                                                                   aaaTGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCa  <  2:1915709/150‑1 (MQ=255)
                                                                                                                                                   aaaTGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCa  <  2:450170/150‑1 (MQ=255)
                                                                                                                                                    aaTGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGa                                                                        <  2:2364232/80‑1 (MQ=255)
                                                                                                                                                    aaTGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGa                                                                        >  1:2364232/1‑80 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCA  >  NC_002947/4475483‑4475779

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