Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,578,474 T→A 7.3% intergenic (+108/‑216) 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,578,4740TA 7.3% 302.3 / 11.8 96intergenic (+108/‑216)PP_3158/PP_5589hypothetical protein/rhs repeat‑associated core domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (46/43);  new base A (4/3);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

CGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT  >  NC_002947/3578326‑3578620
                                                                                                                                                      |                                                                                                                                                  
cgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttttt                                                                                                                                                   <  1:2313280/150‑1 (MQ=255)
         tGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGc                                                                                                                                          >  1:948790/1‑150 (MQ=255)
         tGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGc                                                                                                                                          >  2:60672/1‑150 (MQ=255)
          gTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCa                                                                                                                                         <  1:1440638/150‑1 (MQ=255)
                gatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCgcga                                                                                                                                   >  2:1128060/1‑147 (MQ=255)
                gatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAgt                                                                                                                                   >  2:1582502/1‑150 (MQ=255)
                 atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGa                                                                                                                                     >  1:1074430/1‑147 (MQ=255)
                 atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGa                                                                                                                                     <  2:1074430/147‑1 (MQ=255)
                   gaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCGCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTa                                                                                                                                >  1:2231015/1‑150 (MQ=255)
                       ttCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGAc                                                                                                                            >  1:976020/1‑150 (MQ=255)
                       ttCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGAc                                                                                                                            >  1:2355637/1‑150 (MQ=255)
                        tCCTGCAGGTGCAGTAAGAGGACTCTGGCGCAAGCCTTTTGCGAACTGCCTATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGGGCGAATCAACGTGTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACa                                                                                                                           <  2:266980/150‑1 (MQ=255)
                                  gCAGGAATATGCCTCTGCCGCACGCGCTTTGCGAAATGCCCATACATGAATGAG‑‑TATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAAtt                                                                                                                 <  2:55836/150‑1 (MQ=255)
                                    aGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCa                                                                                                               <  2:142699/150‑1 (MQ=255)
                                     gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATCAAAGCCCGGGGTGAACCAGCGGGAGATGGCAGGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  1:285004/150‑1 (MQ=255)
                                     gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGc                                                                                                                         <  2:2232635/139‑1 (MQ=255)
                                     gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGc                                                                                                                         >  1:2232635/1‑139 (MQ=255)
                                     gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  2:1567934/150‑1 (MQ=255)
                                     gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              >  1:1567934/1‑150 (MQ=255)
                                     gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATCACCGCATCAACGTTTCGTCCATCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  1:296005/150‑1 (MQ=255)
                                     gTAACAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  1:889739/150‑1 (MQ=255)
                                      tAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGtt                                                                                                                                            >  2:1416582/1‑119 (MQ=255)
                                      tAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGtt                                                                                                                                            <  1:1416582/119‑1 (MQ=255)
                                      tAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAg                                                                                                             <  2:178706/150‑1 (MQ=255)
                                         gagGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGa                                                                                                          >  1:1879933/1‑150 (MQ=255)
                                         gagGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTGCATGAGAGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGAGAGGGTATGACAGCCTAGAATTCAAGGGa                                                                                                          >  1:266592/1‑150 (MQ=255)
                                            gCCTCTGCCGCACGCCCTTTGCGAAGTATTTATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAAt                                                                                                       >  1:65159/1‑150 (MQ=255)
                                            gCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GGTAGAGTTGATGAGTGTTCATGAGCCAATCAACGTTTCATCAGGCTGCAGCTGGCATTTTTGCCAGTTGCAGCGAGAGAATGACATCGTAGACGTCAAGGGAAAt                                                                                                       >  2:81895/1‑150 (MQ=255)
                                            gCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTc                                                                                                                <  1:702551/141‑1 (MQ=255)
                                            gCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTc                                                                                                                >  2:702551/1‑141 (MQ=255)
                                               tctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgt                                                                                                    <  2:1879933/150‑1 (MQ=255)
                                                         gCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          <  1:2039984/150‑1 (MQ=255)
                                                         gCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          >  2:2039984/1‑150 (MQ=255)
                                                           ccTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          >  2:216413/1‑148 (MQ=255)
                                                           ccTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          <  1:216413/148‑1 (MQ=255)
                                                           ccTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGt                                                                                        >  1:2181954/1‑150 (MQ=255)
                                                            cTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgt                                                                                                    >  1:237548/1‑137 (MQ=255)
                                                            cTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgt                                                                                                    <  2:237548/137‑1 (MQ=255)
                                                               tGCAAAATGCCCATACATGAATCAG‑‑GATTGAGTTGATGAGTCTTCATGAGCGAATCAACGTTTCATCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGGGAGTGTATGACAGCCGAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCAtt                                                                                    <  1:410035/150‑1 (MQ=255)
                                                                   aaaTGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATc                                                                                <  1:60672/150‑1 (MQ=255)
                                                                     aTGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTTCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGt                                                                              >  2:760969/1‑150 (MQ=255)
                                                                     aTGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGGTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGt                                                                              <  1:1288912/150‑1 (MQ=255)
                                                                       gCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATAATCGTGGCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGa                                                                            <  2:65159/150‑1 (MQ=255)
                                                                         ccATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              >  2:2008050/1‑114 (MQ=255)
                                                                         ccATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  1:2008050/114‑1 (MQ=255)
                                                                             aCATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGa                                                                                                                                     <  1:2088942/87‑1 (MQ=255)
                                                                             aCATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGa                                                                                                                                     >  2:2088942/1‑87 (MQ=255)
                                                                             aCATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGtt                                                                      <  1:1582502/150‑1 (MQ=255)
                                                                              cATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                     <  2:948790/150‑1 (MQ=255)
                                                                              cATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                     <  2:2355637/150‑1 (MQ=255)
                                                                                   atgaG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAg                                                                >  2:54419/1‑150 (MQ=255)
                                                                                       g‑‑gATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAAAAATGCCAGTTGCAg                                                                                                                                        <  1:373532/74‑1 (MQ=25)
                                                                                       g‑‑gATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAAAAATGCCAGTTGCAg                                                                                                                                        >  2:373532/1‑74 (MQ=25)
                                                                                            ttGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTacac                                                         >  1:1599453/1‑150 (MQ=255)
                                                                                               aGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGc                                                                                                                                       <  1:91830/69‑1 (MQ=255)
                                                                                               aGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGc                                                                                                                                       >  2:91830/1‑69 (MQ=255)
                                                                                               aGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGCGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCTCGAAACGTTCGCAATCGTACACGcc                                                      >  2:44967/1‑150 (MQ=255)
                                                                                                    atgaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                     >  2:621330/1‑130 (MQ=255)
                                                                                                    atgaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                     <  1:621330/130‑1 (MQ=255)
                                                                                                        gtgtTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          >  1:671910/1‑105 (MQ=255)
                                                                                                        gtgtTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          <  2:671910/105‑1 (MQ=255)
                                                                                                         tgtTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCt                                            >  1:1848230/1‑150 (MQ=255)
                                                                                                              aTGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTg                                       >  2:66409/1‑150 (MQ=255)
                                                                                                                gAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                     <  2:1599453/150‑1 (MQ=255)
                                                                                                                 aGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGcc                                                      <  2:1852859/132‑1 (MQ=255)
                                                                                                                 aGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGcc                                                      >  1:1852859/1‑132 (MQ=255)
                                                                                                                    gAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAAc                                                                         >  1:683308/1‑110 (MQ=255)
                                                                                                                    gAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAAc                                                                         <  2:683308/110‑1 (MQ=255)
                                                                                                                     aaTCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATg                                                                                                      <  2:771182/80‑1 (MQ=255)
                                                                                                                     aaTCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATg                                                                                                      >  1:771182/1‑80 (MQ=255)
                                                                                                                      aTCAACGTTTCGTCGGGCTGCAAATGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATg                               <  1:81895/150‑1 (MQ=255)
                                                                                                                       tCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATgg                              >  2:1861123/1‑150 (MQ=255)
                                                                                                                         aaCGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGcc                            >  1:1117286/1‑150 (MQ=255)
                                                                                                                          aCGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCa                           >  1:1721053/1‑150 (MQ=255)
                                                                                                                             tttCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTCTTTGGGGTGATGGCCAttc                        >  2:908984/1‑149 (MQ=255)
                                                                                                                               tacactcgCTGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCt                      >  1:705493/8‑150 (MQ=38)
                                                                                                                                agtagggCTGTAACTGGCATTGTTTCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGTACATTTCTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTc                     <  1:760969/146‑1 (MQ=255)
                                                                                                                                  actcgCTGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCg                   >  1:955646/5‑150 (MQ=255)
                                                                                                                                    gggCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                     >  1:1370120/1‑130 (MQ=255)
                                                                                                                                    gggCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                     <  2:1370120/130‑1 (MQ=255)
                                                                                                                                      gCTGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGAt                                <  1:909417/133‑1 (MQ=255)
                                                                                                                                      gCTGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGAt                                >  2:909417/1‑133 (MQ=255)
                                                                                                                                        tGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTg                        <  2:1639774/139‑1 (MQ=255)
                                                                                                                                         gCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa            <  2:2231015/150‑1 (MQ=255)
                                                                                                                                           aaCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGt          >  1:1755557/1‑150 (MQ=255)
                                                                                                                                           aaCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTCGGGGTGATGGCCATTGCTCCGTGTTGAAgg          >  2:1896875/1‑149 (MQ=255)
                                                                                                                                           aaCTGACATTTTTGCCAGTTGCAGCGAGCGTATGACAGCGTAGCATGCACGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGt          >  2:944762/1‑150 (MQ=255)
                                                                                                                                            aCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTa         <  1:1861123/150‑1 (MQ=255)
                                                                                                                                              tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGtt                                                                                                 >  2:40686/1‑60 (MQ=255)
                                                                                                                                              tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGtt                                                                                                 <  1:40686/60‑1 (MQ=255)
                                                                                                                                              tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  2:1247432/1‑150 (MQ=255)
                                                                                                                                                gCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa            <  2:1167179/143‑1 (MQ=255)
                                                                                                                                                gCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa            >  1:1167179/1‑143 (MQ=255)
                                                                                                                                                 cATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATg    <  2:1566330/150‑1 (MQ=255)
                                                                                                                                                 cATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATg    <  1:54419/150‑1 (MQ=255)
                                                                                                                                                  aTTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTgg                                      <  2:2301226/115‑1 (MQ=255)
                                                                                                                                                  aTTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTAAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTgg                                      >  1:2301226/1‑115 (MQ=255)
                                                                                                                                                   tttttGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGt  <  2:1755557/150‑1 (MQ=255)
                                                                                                                                                   tttttGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGt  <  2:976020/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAG‑‑GATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT  >  NC_002947/3578326‑3578620

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