Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 27,584 A→G 5.7% intergenic (‑431/+15) PP_0022 ← / ← PP_0023 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294727,5840AG 5.7% 799.7 / 21.3 244intergenic (‑431/+15)PP_0022/PP_0023hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (121/109);  new base G (9/5);  total (130/114)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 259 total aligned reads displayed.

CAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCAAAA  >  NC_002947/27435‑27733
                                                                                                                                                     |                                                                                                                                                     
cAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCga                                                                                                                                                       <  1:2696851/150‑1 (MQ=255)
  gTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCgaga                                                                                                                                                     >  1:1055924/1‑150 (MQ=255)
   tctcGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGagag                                                                                                                                                    <  1:4071380/150‑1 (MQ=255)
   tctcGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGagag                                                                                                                                                    >  2:3976232/1‑150 (MQ=255)
   tctcGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGagag                                                                                                                                                    >  2:5462237/1‑150 (MQ=255)
         cgtcCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCtt                                                                                                                                              <  1:1533869/150‑1 (MQ=255)
           tcCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCtttt                                                                                                                                            >  2:1929890/1‑150 (MQ=255)
           tcCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCtttt                                                                                                                                            >  2:6430007/1‑150 (MQ=255)
             cGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCtttttt                                                                                                                                          <  1:2193328/150‑1 (MQ=255)
             cGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCtttttt                                                                                                                                          >  2:6996967/1‑150 (MQ=255)
               cACTGAAGGGATACCCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCt                                                                                                                                        <  1:4209048/150‑1 (MQ=255)
               cACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCt                                                                                                                                        <  2:1196243/150‑1 (MQ=255)
                 cTGAAGGGATCACCGATCCGGAGGTAGCTTCAGAATCCCTCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCttg                                                                                                                                      <  2:3914938/150‑2 (MQ=255)
                  tGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCga                                                                                                                                                       >  1:6951980/1‑132 (MQ=255)
                  tGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTaa                                                                                                                                     >  2:5873616/1‑149 (MQ=255)
                  tGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTAc                                                                                                                                     >  2:5137204/1‑150 (MQ=255)
                  tGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTAc                                                                                                                                     >  1:5234234/1‑150 (MQ=255)
                   gAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACt                                                                                                                                    >  1:6646025/1‑150 (MQ=255)
                    aaGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTc                                                                                                                                   >  2:2433323/1‑150 (MQ=255)
                    aaGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTc                                                                                                                                   <  1:3841034/150‑1 (MQ=255)
                      gggATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTcac                                                                                                                                 >  2:7254338/1‑148 (MQ=255)
                         aTAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCt                                                                                                                              >  1:2323394/1‑150 (MQ=255)
                          tAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGTAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCtt                                                                                                                             >  1:752001/1‑150 (MQ=255)
                           aaCCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCCGCGAGACGGGCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGGGTGCTGCTTTTTGCTTACTCCCCCtgg                                                                                                                            >  2:6657563/1‑148 (MQ=255)
                            aCCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCtttt                                                                                                                           >  1:1329608/1‑150 (MQ=255)
                             ccGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCttttt                                                                                                                          <  2:1802711/150‑1 (MQ=255)
                               gATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTtgt                                                                                                                        >  1:4025441/1‑150 (MQ=255)
                               gATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTtgt                                                                                                                        <  1:1404915/150‑1 (MQ=255)
                                      aGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTtgtg                                                                                                                 <  1:6887399/150‑1 (MQ=255)
                                       ggTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTgg                                                                                                                <  2:4690866/150‑1 (MQ=255)
                                       ggTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTgg                                                                                                                <  1:6693404/150‑1 (MQ=255)
                                        gTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTTTGCTTGTggg                                                                                                               >  2:6364818/1‑150 (MQ=255)
                                          aGCTTCAGAATCCATCACGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGt                                                                                                             <  1:4986088/150‑1 (MQ=255)
                                            cTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGttt                                                                                                           <  2:3586451/150‑1 (MQ=255)
                                             ttCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTc                                                                                                          >  2:5046534/1‑150 (MQ=255)
                                              tCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCa                                                                                                         <  2:647096/150‑1 (MQ=255)
                                               cAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAg                                                                                                        <  1:4306477/150‑1 (MQ=255)
                                               cAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAg                                                                                                        <  1:3540143/150‑1 (MQ=255)
                                                 gAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGt                                                                                                      >  2:4042744/1‑150 (MQ=255)
                                                   aTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTcc                                                                                                    <  1:3673518/150‑1 (MQ=255)
                                                      cATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGa                                                                                                 <  2:5537833/150‑1 (MQ=255)
                                                       aTCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGAt                                                                                                >  2:654921/1‑150 (MQ=255)
                                                          aagccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGt                                                                                             <  2:816179/150‑1 (MQ=255)
                                                           agccacgcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTc                                                                                            <  2:1015337/150‑1 (MQ=255)
                                                           agccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTc                                                                                            >  1:4003432/1‑150 (MQ=255)
                                                               acgcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTc                                                                                        <  1:2717972/148‑1 (MQ=255)
                                                               aagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTc                                                                                        <  1:7023257/150‑1 (MQ=255)
                                                                agcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCTGGTCCGGATGGTCCTTCt                                                                                       >  2:5703781/1‑150 (MQ=255)
                                                                 gcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTg                                                                                      <  1:86785/150‑1 (MQ=255)
                                                                 gcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTg                                                                                      <  1:7200945/150‑1 (MQ=255)
                                                                    gAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCtt                                                                                   <  1:2308833/150‑1 (MQ=255)
                                                                     aaGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCttt                                                                                  >  1:6376778/1‑150 (MQ=255)
                                                                        aCTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCa                                                                               >  1:1736276/1‑150 (MQ=255)
                                                                         cTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAgg                                                                             <  1:2877411/150‑1 (MQ=255)
                                                                          tCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAgg                                                                             >  2:4922169/1‑150 (MQ=255)
                                                                           cACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAggg                                                                            <  1:1675907/150‑1 (MQ=255)
                                                                             cgacgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGt                                                                          <  2:7227071/150‑1 (MQ=255)
                                                                                   ggTCCGATTGATCAGGGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCg                                                                    <  2:5834287/150‑1 (MQ=255)
                                                                                   ggTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCg                                                                    >  2:2329201/1‑150 (MQ=255)
                                                                                       cGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGc                                                                <  1:3931859/150‑1 (MQ=255)
                                                                                            gATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCtttt                                                           <  2:3731945/150‑1 (MQ=255)
                                                                                              tCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGt                                                         >  2:5563691/1‑150 (MQ=255)
                                                                                              tCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGt                                                         >  1:941361/1‑150 (MQ=255)
                                                                                              tCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGt                                                         >  2:4434628/1‑150 (MQ=255)
                                                                                              tCAGCGAGACGGTCGTGAGGGATTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCGTTTTGt                                                         >  1:1228604/1‑150 (MQ=255)
                                                                                                aGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTcc                                                       <  1:460749/150‑1 (MQ=255)
                                                                                                aGCGAGACGGTCGAGAGGGTTTGCACACAGCCCCCCAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTcc                                                       <  2:2146014/150‑1 (MQ=255)
                                                                                                 gCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCg                                                      >  2:658498/1‑150 (MQ=255)
                                                                                                   gagaCGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGa                                                    <  1:3442530/150‑1 (MQ=255)
                                                                                                    agaCGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAg                                                   <  2:407311/150‑1 (MQ=255)
                                                                                                             tGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCtt                                          >  1:3220603/1‑150 (MQ=255)
                                                                                                             tGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCGTTTGTCCGGAGTGCCATCtt                                          >  1:579562/1‑150 (MQ=255)
                                                                                                              gAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCGTc                                         >  1:5954457/1‑150 (MQ=255)
                                                                                                                gggTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCAt                                       >  2:3384328/1‑150 (MQ=255)
                                                                                                                     tGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTg                                  >  1:2469234/1‑150 (MQ=255)
                                                                                                                     tGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCAGCTGCATTGGTg                                  >  2:486987/1‑150 (MQ=255)
                                                                                                                      gCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTgc                                 <  1:2018430/150‑1 (MQ=255)
                                                                                                                        acacaGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCGTCTTCATTGGTgcgc                               >  2:1383225/1‑150 (MQ=255)
                                                                                                                          acaGCCCCAAAAAAAAAGCAGCACCCGAGAGCTGCTATTTTCTTACTCCCCCTCTTTGTGCTAGTGGGGTTTCAGGTCCGGATCGTCCTTCTGCTTTTCAGGCGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCtt                             <  2:2654943/150‑1 (MQ=255)
                                                                                                                           caGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCttt                            <  2:6087609/150‑1 (MQ=255)
                                                                                                                            aGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTg                           <  2:3067860/150‑1 (MQ=255)
                                                                                                                             gCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCTCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTgg                          <  1:5648444/150‑1 (MQ=255)
                                                                                                                               cccAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGc                        >  2:2200611/1‑150 (MQ=255)
                                                                                                                                 taagaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTTTTACTCCCCCTGTTTGTGCGTGTGGGGTTTCAGGTCCGGATGGTCCGTCTGCTTGTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCCTTGGTGCGCTTTGGGCCt                      >  2:4415239/5‑150 (MQ=38)
                                                                                                                                 ccaaaaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCt                      <  2:7137853/148‑1 (MQ=255)
                                                                                                                                 cAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCCTCTTCATTGGTGCGCTTTGGGCCt                      >  2:1392355/1‑150 (MQ=255)
                                                                                                                                 cAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCt                      >  1:4633195/1‑150 (MQ=255)
                                                                                                                                 cAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCt                      <  1:3976232/150‑1 (MQ=255)
                                                                                                                                 cAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCt                      >  1:2996298/1‑150 (MQ=255)
                                                                                                                                  aaaaaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTg                     >  1:5557354/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCt                  >  1:6762674/1‑150 (MQ=255)
                                                                                                                                       aaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTg                >  1:6176627/1‑150 (MQ=255)
                                                                                                                                        aaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGa                                                                                                 >  2:2117026/1‑68 (MQ=255)
                                                                                                                                         aaGCAGCTCTCGGGAGCTGCTTTTTTTTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATgg                                                                                              >  1:5554119/1‑70 (MQ=37)
                                                                                                                                          accagcTCCCGAGAGCTGCAATTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGTGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATg             <  2:1329608/148‑1 (MQ=255)
                                                                                                                                                tCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTGTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCATCATTGGTGAGCTTTGGGCCTCGCTTGATGCTCATg       >  2:2699339/1‑150 (MQ=255)
                                                                                                                                                  ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTc                                                                                                          >  1:6502034/1‑49 (MQ=255)
                                                                                                                                                  ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCa     >  2:1155673/1‑150 (MQ=255)
                                                                                                                                                   cGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCaa    <  1:5498577/150‑1 (MQ=255)
                                                                                                                                                     agagCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCaaaa  <  1:4156690/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCAAAA  >  NC_002947/27435‑27733

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