Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 870,234 A→C 6.5% A64A (GCT→GCG PP_5456 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947870,2340AC 6.5% 228.6 / 4.1 83A64A (GCT→GCGPP_5456hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base C (2/4);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGCGTT  >  NC_002947/870085‑870383
                                                                                                                                                     |                                                                                                                                                     
tCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGa                                                                                                                                                       >  1:1708303/1‑150 (MQ=255)
  ggCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCTGGGTTTTTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAgc                                                                                                                                                     <  2:928351/150‑1 (MQ=255)
  ggCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGAAGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAgc                                                                                                                                                     <  2:2185216/150‑1 (MQ=255)
              gCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGc                                                                                                                                         <  1:1339638/150‑1 (MQ=255)
                 cGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTc                                                                                                                                      >  2:1831847/1‑150 (MQ=255)
                       ggAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCggg                                                                                                                                >  1:395516/1‑150 (MQ=255)
                       ggAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTATTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGACAACGACGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCggg                                                                                                                                >  2:1377058/1‑150 (MQ=255)
                        gagaCGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGc                                                                                                                               >  2:1092633/1‑150 (MQ=255)
                             gTCCTTATAATCGGGGGGCGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTAGTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCGGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCGTCGGCGGGCAGCTc                                                                                                                          >  2:607824/1‑150 (MQ=255)
                             gTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCGGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCGTCGGCGGGCAGCTc                                                                                                                          >  1:2288487/1‑150 (MQ=255)
                              tCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTcg                                                                                                                         >  1:2102691/1‑150 (MQ=255)
                                 ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCg                                                                                                                                            >  1:320739/1‑128 (MQ=255)
                                 ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCg                                                                                                                                            <  2:320739/128‑1 (MQ=255)
                                      aTCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGAGAACAGAAACGCGGGCTTCGGCGGGCAGCTCGCGCCccag                                                                                                                 >  2:2017996/1‑150 (MQ=255)
                                       tCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagc                                                                                                                >  2:1809754/1‑150 (MQ=255)
                                          gggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAATCCGGGCTTCGGCGGGCAGCTCGCGCCccagccag                                                                                                             <  1:2017996/150‑1 (MQ=255)
                                          gggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTGTGTCAGGCAGCGAGGAGGCGCTGGTGCCGACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagccag                                                                                                             <  2:2427843/150‑1 (MQ=255)
                                          gggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGCGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagccag                                                                                                             >  1:179478/1‑150 (MQ=255)
                                                 aTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCt                                                                                                      <  1:1092633/150‑1 (MQ=255)
                                                               gcgGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTAAGGGCGCACACGCACGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATc                                                                                        >  2:2433452/1‑150 (MQ=255)
                                                                 gggCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTg                                                                                      >  2:2007646/1‑150 (MQ=255)
                                                                  ggCAGAAATGTCGCCATTCGTGGTGTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTAAGACCGACGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGc                                                                                     <  1:147411/150‑1 (MQ=255)
                                                                     aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGaa                                                                                                   >  1:789149/1‑133 (MQ=255)
                                                                     aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGaa                                                                                                   <  2:789149/133‑1 (MQ=255)
                                                                     aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGAt                                                                                         >  2:603837/1‑143 (MQ=255)
                                                                     aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGAt                                                                                         <  1:603837/143‑1 (MQ=255)
                                                                       aaaTGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCCGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCa                                                                                <  1:1840185/150‑1 (MQ=255)
                                                                         aTGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGACCTGCACCCAGt                                                                              <  1:1809754/150‑1 (MQ=255)
                                                                            tCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGa                                                                           <  2:395516/150‑1 (MQ=255)
                                                                            tCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGa                                                                           >  1:2210148/1‑150 (MQ=255)
                                                                            tCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGa                                                                           <  2:1708303/150‑1 (MQ=255)
                                                                             cGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGAt                                                                          <  1:1831847/150‑1 (MQ=255)
                                                                              gCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATg                                                                         >  1:1338583/1‑150 (MQ=255)
                                                                               ccATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGc                                                                        >  1:2338759/1‑150 (MQ=255)
                                                                                cATTCGTGGGTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGcc                                                                       >  2:1852536/1‑150 (MQ=255)
                                                                                      tGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccaggccag                                                                                                        >  2:2108673/1‑111 (MQ=255)
                                                                                      tGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccaggccag                                                                                                        <  1:2108673/111‑1 (MQ=255)
                                                                                      tGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtgg                                                                 >  2:200087/1‑150 (MQ=255)
                                                                                       ggTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGGGGCCGCGCCCACCTGGGCTTCGGCCGGCATCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtggt                                                                <  2:866169/150‑1 (MQ=255)
                                                                                        gTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtggtg                                                               >  2:314698/1‑150 (MQ=255)
                                                                                         tttttATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACACGGGCTTCGGAGGGCAGCTCTACTCGCAGCCATGTTAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGc                                                              <  1:1488086/150‑1 (MQ=255)
                                                                                              aTTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGc                                                         <  2:2338759/150‑1 (MQ=255)
                                                                                               ttGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCt                                                        <  1:2007646/150‑1 (MQ=255)
                                                                                                   aGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCAt                                                    >  2:2235413/1‑150 (MQ=255)
                                                                                                      gTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGc                                                 <  1:2235413/150‑1 (MQ=255)
                                                                                                        ttGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCCAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGATGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCtt                                               <  2:1253089/150‑1 (MQ=255)
                                                                                                           tCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                            <  2:2102691/150‑1 (MQ=255)
                                                                                                            cAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGCGCGCAAGCCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                           <  2:686559/150‑1 (MQ=255)
                                                                                                             aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAACCCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTgg                                                          <  1:1152649/134‑1 (MQ=255)
                                                                                                             aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTgg                                                          >  2:1152649/1‑134 (MQ=255)
                                                                                                              ggCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGATGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGc                                         >  1:1328479/1‑150 (MQ=255)
                                                                                                               gCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGcc                                        >  2:2219909/1‑150 (MQ=255)
                                                                                                               gCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGcc                                        >  1:2333703/1‑150 (MQ=255)
                                                                                                               gCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCACCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGcc                                        >  1:1624330/1‑150 (MQ=255)
                                                                                                                cAGCGAGGAGGCGCTGGTGCCAACGGCGGCTTGAGCGAAAAGAGTGCGGGGCTTGGGCGGGCAGCTCGCACCATCCCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCa                                       <  1:1740672/150‑1 (MQ=255)
                                                                                                                 aGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAg                                      >  1:1122023/1‑150 (MQ=255)
                                                                                                                    gaggagGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGt                                                                              >  2:2357753/1‑107 (MQ=255)
                                                                                                                    gaggagGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGt                                                                              <  1:2357753/107‑1 (MQ=255)
                                                                                                                      ggaggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCAtt                                 >  2:1591966/1‑150 (MQ=255)
                                                                                                                          gcgcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTc                                                                             >  2:2414871/1‑102 (MQ=255)
                                                                                                                          gcgcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTc                                                                             <  1:2414871/102‑1 (MQ=255)
                                                                                                                           cgcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCa                            <  1:1591966/150‑1 (MQ=255)
                                                                                                                            gcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCTGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                                  >  1:2236452/1‑111 (MQ=255)
                                                                                                                            gcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCTGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                                  <  2:2236452/111‑1 (MQ=255)
                                                                                                                            gcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgt                                                                   >  2:95282/1‑110 (MQ=255)
                                                                                                                            gcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgt                                                                   <  1:95282/110‑1 (MQ=255)
                                                                                                                                    cAACGGCGGTTTGCGCGAGCGCAGACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGcc                             <  2:1553830/140‑1 (MQ=255)
                                                                                                                                    cAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGcc                             >  1:1553830/1‑140 (MQ=255)
                                                                                                                                    cAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCa                   >  1:1879001/1‑150 (MQ=255)
                                                                                                                                     aaCGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGAt                                                                                         <  1:2371757/79‑1 (MQ=255)
                                                                                                                                     aaCGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGAt                                                                                         >  2:2371757/1‑79 (MQ=255)
                                                                                                                                      aCGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAg                 <  2:2225794/150‑1 (MQ=255)
                                                                                                                                         cgggTTTGCGCTCGCGCACACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                           <  1:2102103/119‑1 (MQ=255)
                                                                                                                                         cgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                           >  2:2102103/3‑121 (MQ=255)
                                                                                                                                          gggTTTGCGCTCGCGCAAACCGGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg             <  1:1900399/149‑1 (MQ=255)
                                                                                                                                          cggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg             <  1:200087/150‑1 (MQ=255)
                                                                                                                                          cggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg             <  1:1377058/150‑1 (MQ=255)
                                                                                                                                           ggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCggg            <  1:1762770/150‑1 (MQ=255)
                                                                                                                                            gTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGc           >  2:404385/1‑150 (MQ=255)
                                                                                                                                             tttGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCa          >  1:2169331/1‑150 (MQ=255)
                                                                                                                                                gcgcGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAAttt                                                                                                >  2:1581182/1‑61 (MQ=255)
                                                                                                                                                gcgcGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAAttt                                                                                                <  1:1581182/61‑1 (MQ=255)
                                                                                                                                                gcgcGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATaa       <  1:2219909/150‑1 (MQ=255)
                                                                                                                                                gcgcGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATaa       >  1:443714/1‑150 (MQ=255)
                                                                                                                                                     aGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGCGtt  <  1:2356259/150‑1 (MQ=255)
                                                                                                                                                     aGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGCGtt  <  2:179478/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGCGTT  >  NC_002947/870085‑870383

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