Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,583,420 T→A 12.1% intergenic (+105/‑219) PP_3158 → / → PP_5589 hypothetical protein/rhs repeat‑associated core domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,583,4200TA12.1% 709.9 / 38.8 281intergenic (+105/‑219)PP_3158/PP_5589hypothetical protein/rhs repeat‑associated core domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (118/121);  new base A (36/3);  total (157/124)
Fisher's exact test for biased strand distribution p-value = 1.35e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 303 total aligned reads displayed.

CGCGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT  >  NC_002947/3583273‑3583569
                                                                                                                                                   |                                                                                                                                                     
cgcgCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttt                                                                                                                                                     <  2:6090181/150‑1 (MQ=255)
 gcgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAtttt                                                                                                                                                    <  1:7001919/150‑1 (MQ=255)
 gcgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGGAtttt                                                                                                                                                    >  1:4373817/1‑150 (MQ=255)
  cgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttttt                                                                                                                                                   <  2:5116269/150‑1 (MQ=255)
  cgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttttt                                                                                                                                                   <  1:2359984/150‑1 (MQ=255)
  cgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttttt                                                                                                                                                   <  2:1413845/150‑1 (MQ=255)
         ccTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAAAAATGCCAGtt                                                                                                                                            <  1:2728122/150‑1 (MQ=255)
          cTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAAAAATGCCAGTTg                                                                                                                                           <  2:2693968/150‑1 (MQ=255)
            gTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCa                                                                                                                                         <  2:6463865/150‑1 (MQ=255)
                 tgatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAg                                                                                                                                    >  2:4862317/1‑150 (MQ=255)
                  gatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGGATTTTTTGCCGTTGCAGCGAgt                                                                                                                                   >  2:5228418/1‑150 (MQ=255)
                   atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAgtg                                                                                                                                  >  2:5406934/1‑150 (MQ=255)
                   atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATTGCATTTTTTCGAGTTTCAGCGAgtt                                                                                                                                  >  2:3339312/1‑149 (MQ=255)
                    tgaACTTCCTGCCGGTGCAGTACGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGCGGTTTGTGTTGATGAGTGTTCTTGCGCGAATCAAGGTTTGGTCGGGCTGCAACTGGAATTTTTGCCAGTTGCAGCGAgtgt                                                                                                                                 <  1:1017131/150‑1 (MQ=255)
                    tgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAgtgt                                                                                                                                 >  2:5764149/1‑150 (MQ=255)
                     gaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCCGCGGGTGTa                                                                                                                                >  1:6823313/1‑150 (MQ=255)
                      aaCTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTAt                                                                                                                               <  2:7408537/150‑1 (MQ=255)
                        cTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTACCAGTTGCATCGAGTGTATGa                                                                                                                             <  2:5082097/150‑1 (MQ=255)
                         ttCCTGCCGGTGCAGTAAGCGGCCTCGGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGCGTTGATGAGTGTTCCTGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGAc                                                                                                                            <  2:2134595/150‑1 (MQ=255)
                              gCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCt                                                                                                                       <  1:3693534/150‑1 (MQ=255)
                               ccGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTa                                                                                                                      <  1:6782051/150‑1 (MQ=255)
                               ccGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCAGACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTa                                                                                                                      <  2:1634687/150‑1 (MQ=255)
                                 ggTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCGGCCACAGGCCGTTTTTCCCGGTGCCGCGCGGGTATGGCACCCgaga                                                                                                                    >  1:2007591/1‑150 (MQ=255)
                                  gTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGaa                                                                                                                   <  2:2718990/150‑1 (MQ=255)
                                  gTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCAGGAGCGAATCAACTTTTCGTCGGGCTGCAACTGGCATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGaa                                                                                                                   <  2:6360932/150‑1 (MQ=255)
                                    gCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTTCAGCGGGTGTTTGAAAGCCTACaaat                                                                                                                 >  2:4703534/1‑148 (MQ=255)
                                      aGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGc                                                                                                                                       <  1:3589237/126‑1 (MQ=255)
                                       gTAAGTGGCCTCTGCCGCACGCCCTTTGCGAACTTCCCATACATGAATGAGGATTGATTTGATGAGTGTTCTTGAGCGATTATACGTTTCGTCGGGCTGAAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  2:1893470/150‑1 (MQ=255)
                                       gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGCGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  2:3959234/150‑1 (MQ=255)
                                        tAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAAATGCAATTTTTCCAGATTCAGCGAGTGTATTACAGCCAAGAAATAAAg                                                                                                             >  1:2161088/1‑150 (MQ=255)
                                        tAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATGGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAg                                                                                                             <  1:6129266/150‑1 (MQ=255)
                                        tAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAAGGTGGTGTTTGTTGATGAGTGTTCATGAGCGAAAAAACTTTTCGTCGGGCTGCAACTGGCATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAg                                                                                                             <  1:5650528/150‑1 (MQ=255)
                                        tAAGAGGCCTCGGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGTATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGAAACTGGAATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAg                                                                                                             <  1:4510111/150‑1 (MQ=255)
                                         aaGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCCACAGGCCTTTTTGCCCGGTGCAGCGCGGGTATGACAGCCGAGAATTCAAgg                                                                                                            >  2:1464860/1‑150 (MQ=255)
                                         aaGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAgg                                                                                                            <  2:6282505/150‑1 (MQ=255)
                                            agGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCGAGAATTACAGGGag                                                                                                         >  2:504918/1‑149 (MQ=255)
                                             ggCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGaaa                                                                                                        <  2:5931243/150‑1 (MQ=255)
                                             ggCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCGAGAATTCAAGGGaaa                                                                                                        >  1:5113194/1‑150 (MQ=255)
                                               ccTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGAGTGAGTTGATGAGTGTTCAGGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCCGTTTTGCCCGTTGCAGCGGGTGGATGACAGCCGAGACTTCAAGGGcaaag                                                                                                      >  2:2515395/1‑148 (MQ=255)
                                                ctctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGACATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATgg                                                                                                     >  2:5599569/1‑150 (MQ=255)
                                                 tctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgt                                                                                                    >  2:7349417/1‑150 (MQ=255)
                                                 tctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACATGCATTTTTTGCAGTTGCCGCGGGTGTATGACAGCCTAGAATAcaagggcaagggg                                                                                                    >  1:4556937/1‑149 (MQ=255)
                                                 tctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTGCACTTTCAGCGAGTGTATTAAACCCTAGAATTCAAGGGaaaatgt                                                                                                    >  1:1236463/1‑146 (MQ=255)
                                                  ctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGACTCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgtg                                                                                                   <  2:71449/150‑1 (MQ=255)
                                                  ctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgtg                                                                                                   <  2:6787375/150‑1 (MQ=255)
                                                     ccGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAAAGGCATTTTTGGCAGTTGCACCGAGTGTATTTCACCCGAGAATTCAAGGGTAATGGTGttg                                                                                                >  1:3414014/1‑149 (MQ=255)
                                                      cGCACGCCCTGTGCGAAATGCCCATACTTGAATGAGGATTTAGTTGAGGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGATGCAAGTGGAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTtct                                                                                               <  2:4378063/150‑1 (MQ=255)
                                                           gCCCTTTGCGAAATGCCAATACATTATTGTGTTTTGAGTTGATGTGTTTTATTGAGCGAATCAACGTTTCGTGGGGCTGCAACTGGAATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          <  1:5228418/150‑1 (MQ=255)
                                                             ccTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAATTTTTCCAGTTGCAGCGAGTGTATTACACCCGAGAATTAAAGGGAAATTGTGTTGTGTAAGGt                                                                                        >  1:2897073/1‑150 (MQ=255)
                                                             ccTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCCGTTTTGCCAGTTGCAGCGAGTGTATGACAGCCGAGAATTCAAGGGAAATGGTGTTCTCTAAGGt                                                                                        >  2:3227776/1‑150 (MQ=255)
                                                               tttGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCACTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTTTTCTCTAAGGTCa                                                                                      >  2:231983/1‑150 (MQ=255)
                                                               tttGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTTCCAGTTGCAGCGAGTGTATTACAGCCAAGAATTAAAGGGAAATGGTGGTCTCTAAGGGCa                                                                                      >  1:1625341/1‑150 (MQ=255)
                                                                 tGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCAtt                                                                                    <  2:3431409/150‑1 (MQ=255)
                                                                 tGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGGGTTCTCTAAGGTCAtt                                                                                    >  2:4770840/1‑150 (MQ=255)
                                                                      aaTGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACATGCATTTTTTCCAGTTGAAGCGAGTGTAAGACAGCCGAGAAATCAAGGGAAATGGTGTTAGCTAAGGTGAATaaacc                                                                               >  2:5279903/1‑147 (MQ=255)
                                                                       aTGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGt                                                                              <  2:6569245/150‑1 (MQ=255)
                                                                         gCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACTTTTCGTCGGGCTGCAACTGGAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGa                                                                            <  1:1720827/150‑1 (MQ=255)
                                                                         gCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGa                                                                            <  2:3580065/150‑1 (MQ=255)
                                                                           ccATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCAGTTGCAGCGCGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATAGTGaaa                                                                          >  2:1931492/1‑150 (MQ=255)
                                                                            cATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACATGCAATTTTTCCACTTTCAGCGAGTGTATGTCAGCCAAGAATTAAATGGTAATTGTTTTGTCTAAGGTGATTAATAGTGaaaa                                                                         >  1:4452009/1‑149 (MQ=255)
                                                                                cATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                     <  1:550433/150‑1 (MQ=255)
                                                                                  tgaatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAATTTTTCCACTTTCAGCGAGTGTATTACAGCCGAGAAATAAAGGGAAAAGGTGTTCTCTAAGGTCATTAATAGTGAAACGTTCGc                                                                   >  2:3273208/1‑150 (MQ=255)
                                                                                   gaatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAAAGGCAATTTTTGCAGTTTCAGCGGGTGTATTTAAACCAAGAAATAAACGGCAAATGTTTTGTCTAAGGTGATTAAAAGTTAAACCTTAGCa                                                                  >  1:2800103/1‑150 (MQ=255)
                                                                                     atgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTTGCATTTTTTCCAGTTTCAGCGCGTGTATTACCTCCTAGAAATAAATGGAAATGGTGTTCTCTAATGTCATTAATCGTGAAAAGTTCGCAAg                                                                >  2:6222303/1‑150 (MQ=255)
                                                                                         ggATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTa                                                            <  2:1162146/150‑1 (MQ=255)
                                                                                          gATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTTCCAGTTGCAGCGAGTGTATGACAGCCAAGAAAAAAAGGGAAATGGTGATCTATAAGGTCAATAATAGTGAAAACTTAGCAAGCGTac                                                           >  1:2974722/1‑150 (MQ=255)
                                                                                               aGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              >  2:6134054/1‑94 (MQ=255)
                                                                                               aGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGcc                                                      <  1:7349417/150‑1 (MQ=255)
                                                                                                 ttGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAAAGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCGAGAATTCAAGGGCAAAGGTGTTGTCTAAGGTCATTAATAGTGAAAAGTTCGCAAGCGTACACGccca                                                    >  1:4876050/1‑149 (MQ=255)
                                                                                                         tgtTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTATAAGGTCAATAATAGGGAAAAGTTCGCAAGCGTACACGCCCAAAGGCGCt                                            >  2:7088226/1‑150 (MQ=255)
                                                                                                         tgtTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATGGCATTTTTGCCAGTTGCAGCGAGTTTATGACAGCCTAGAATTCAATGGAAATGGTGGTCTCTAAGGTCATTAATCGTGGAACGTTCGCAAGCGTACACGCCCGACGgcgcg                                            >  2:1818640/1‑149 (MQ=255)
                                                                                                                gAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCTAGCGTACACGCCCTACGGCGCTATTTggg                                     <  2:1065168/150‑1 (MQ=255)
                                                                                                                gAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                     <  1:6630263/150‑1 (MQ=255)
                                                                                                                gAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                     <  2:4271722/150‑1 (MQ=255)
                                                                                                                 aGCGAATCAACGTTTCGTCGGGCTGCAACAGGCCTTTTTGCCAGTTTCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATTGTGTTGTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGGACACGCCCTACGGCGCGAGTTgggg                                    >  2:4576697/1‑150 (MQ=255)
                                                                                                                     aaTCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGAt                                <  1:3961294/150‑1 (MQ=255)
                                                                                                                     aaTCAACGTTTCGTCGGGCTGCAACTGGAATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGAt                                <  1:247448/150‑1 (MQ=255)
                                                                                                                       tCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATgg                              <  2:2924108/150‑1 (MQ=255)
                                                                                                                        cAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAATGTCATTAATAGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGc                             >  2:2382305/1‑150 (MQ=255)
                                                                                                                         aaCGTTTCGTCGGGCTGCAACAGGCAATTTTGCCACTTGCAGCGAGTGTAAGAAAGCCAATAAATAAAGGGAAATTGTGTTCTCTAAGGTCATTAATCGTGTAAAGTTAGCAACCGTAAACGCCCAACCGCGCTAATTTGGGTGTTGGcc                            >  1:4223401/1‑150 (MQ=255)
                                                                                                                             tttCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTg                        <  2:4449337/150‑1 (MQ=255)
                                                                                                                                 gtcgGGCTGCAACTGGCATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTcc                    <  1:6471305/150‑1 (MQ=255)
                                                                                                                                  tcgGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCg                   <  2:3734680/150‑1 (MQ=255)
                                                                                                                                   cgGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCgt                  <  2:1856022/150‑1 (MQ=255)
                                                                                                                                      gCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTg                                                                             >  1:5080200/1‑88 (MQ=255)
                                                                                                                                       cTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTg              >  1:1886065/1‑150 (MQ=255)
                                                                                                                                       cTGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTg              >  2:854962/1‑150 (MQ=255)
                                                                                                                                       cTGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTAGTTGGGGTGATGGCCATTGCTCCGTGTTg              >  2:4042222/1‑150 (MQ=255)
                                                                                                                                          cAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAg           >  1:5444627/1‑150 (MQ=255)
                                                                                                                                          cAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAg           >  1:3199906/1‑150 (MQ=255)
                                                                                                                                            aCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTa         >  2:6928764/1‑150 (MQ=255)
                                                                                                                                            aCTGGATTTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTa         <  1:799534/150‑1 (MQ=255)
                                                                                                                                              tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCgtgt                <  1:7144634/141‑1 (MQ=255)
                                                                                                                                              tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  2:2622135/1‑150 (MQ=255)
                                                                                                                                              tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  1:5534452/1‑150 (MQ=255)
                                                                                                                                               ggCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACaa      <  2:2441509/150‑1 (MQ=255)
                                                                                                                                               ggCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACaa      >  1:5790437/1‑150 (MQ=255)
                                                                                                                                                 cATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATg    <  1:1095240/150‑1 (MQ=255)
                                                                                                                                                  aTTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATgg   >  2:3907781/1‑150 (MQ=255)
                                                                                                                                                   tttttGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGt  <  2:1734757/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT  >  NC_002947/3583273‑3583569

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