Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,520,127 T→A 5.7% intergenic (+703/‑177) PP_5580 → / → PP_5582 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,520,1270TA 5.7% 286.3 / 5.6 87intergenic (+703/‑177)PP_5580/PP_5582hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/47);  new base A (3/2);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTTGACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTGTGTCTCAG  >  NC_002947/3519978‑3520275
                                                                                                                                                     |                                                                                                                                                    
tttGACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTAt                                                                                                                                                      <  1:146487/150‑1 (MQ=255)
     cGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTAGATAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATgcaag                                                                                                                                                 <  1:1230345/150‑1 (MQ=255)
        tGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATgcaagcaa                                                                                                                                              <  2:100276/150‑1 (MQ=255)
         gAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAAt                                                                                                                                             >  1:592989/1‑150 (MQ=255)
          aaGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATa                                                                                                                                            <  1:2140448/150‑1 (MQ=255)
            ggAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGt                                                                                                                                          >  1:282063/1‑150 (MQ=255)
            ggAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGt                                                                                                                                          <  2:1081699/150‑1 (MQ=255)
                   cTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCCATAGTGCTCCTg                                                                                                                                   <  2:1842899/150‑1 (MQ=255)
                   cTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTg                                                                                                                                   >  1:1822581/1‑150 (MQ=255)
                     aTTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGt                                                                                                                                 <  1:1031421/150‑1 (MQ=255)
                      ttAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTa                                                                                                                                <  2:1822581/150‑1 (MQ=255)
                       tAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTAc                                                                                                                               >  2:2087954/1‑150 (MQ=255)
                                tGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATAtt                                                                                                                      <  2:968806/150‑1 (MQ=255)
                                  gggTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCa                                                                                                                    <  1:287102/150‑1 (MQ=255)
                                     tCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTg                                                                                                                 <  2:2382181/150‑1 (MQ=255)
                                        tatTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGc                                                                                                              <  1:2160966/150‑1 (MQ=255)
                                         atTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATgcaag                                                                                                                                                 >  2:2258917/1‑114 (MQ=255)
                                         atTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATgcaag                                                                                                                                                 <  1:2258917/114‑1 (MQ=255)
                                          tttGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCga                                                                                                            >  2:1515370/1‑150 (MQ=255)
                                                aaTCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGc                                                                                                      >  2:1571139/1‑150 (MQ=255)
                                                aaTCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGc                                                                                                      <  2:282063/150‑1 (MQ=255)
                                                aaTCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGCGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGc                                                                                                      >  2:2053318/1‑150 (MQ=255)
                                                 aTCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATa                                                                                                                                            >  2:1937084/1‑111 (MQ=255)
                                                 aTCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATa                                                                                                                                            <  1:1937084/111‑1 (MQ=255)
                                                     aGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATa                                                                                                                                            <  1:806963/107‑1 (MQ=255)
                                                     aGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATa                                                                                                                                            >  2:806963/1‑107 (MQ=255)
                                                           aGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTa                                                                                                                                <  2:1179870/113‑1 (MQ=255)
                                                           aGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTa                                                                                                                                >  1:1179870/1‑113 (MQ=255)
                                                           aGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGt                                                                                           <  1:88429/150‑1 (MQ=255)
                                                           aGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGGGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGt                                                                                           >  1:2291152/1‑150 (MQ=255)
                                                             atatatCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTg                                                                                         <  2:2119521/150‑1 (MQ=255)
                                                              tatatCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAg                                                                                                                                           >  2:1415658/1‑99 (MQ=255)
                                                              tatatCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAg                                                                                                                                           <  1:1415658/99‑1 (MQ=255)
                                                                 atCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTAc                                                                                                                               >  2:1827859/1‑108 (MQ=255)
                                                                 atCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTAc                                                                                                                               <  1:1827859/108‑1 (MQ=255)
                                                                 atCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCtt                                                                                     <  1:1693202/150‑1 (MQ=255)
                                                                    aTGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTc                                                                                                                                      <  1:570856/98‑1 (MQ=255)
                                                                    aTGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTc                                                                                                                                      >  2:570856/1‑98 (MQ=255)
                                                                          ccTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTTGAGATATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTa                                                                            <  2:592989/150‑1 (MQ=255)
                                                                           cTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTAc                                                                           >  1:499909/1‑150 (MQ=255)
                                                                           cTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTAc                                                                           <  1:918650/150‑1 (MQ=255)
                                                                           cTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTAc                                                                           <  1:2087954/150‑1 (MQ=255)
                                                                             ggaggaCGATGGTGAATATTGGTGCTATGAATATAGGCGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGa                                                                         <  1:2053318/150‑1 (MQ=255)
                                                                                        gTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTAt                                                                                 >  2:37152/1‑131 (MQ=255)
                                                                                        gTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTAt                                                                                 <  1:37152/131‑1 (MQ=255)
                                                                                          gAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACa                                                                                                                              <  2:404290/84‑1 (MQ=255)
                                                                                          gAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACa                                                                                                                              >  1:404290/1‑84 (MQ=255)
                                                                                          gAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCaaaa                                                            <  1:1258091/150‑1 (MQ=255)
                                                                                              atTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGc                                                        >  2:1255089/1‑150 (MQ=255)
                                                                                                  gTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCt                                                    >  1:2071490/1‑150 (MQ=255)
                                                                                                  gTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCt                                                    >  1:2363666/1‑150 (MQ=255)
                                                                                                  gTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCt                                                    >  2:1281934/1‑150 (MQ=255)
                                                                                                  gTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCt                                                    >  1:919957/1‑150 (MQ=255)
                                                                                                         gAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACaaa                                             <  1:1515370/150‑1 (MQ=255)
                                                                                                          aaTATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTAt                                                                                 <  1:17728/113‑1 (MQ=255)
                                                                                                          aaTATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTAt                                                                                 >  2:17728/1‑113 (MQ=255)
                                                                                                           atatAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCaaaa                                                            >  2:1291734/1‑133 (MQ=255)
                                                                                                           atatAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCaaaa                                                            <  1:1291734/133‑1 (MQ=255)
                                                                                                            tataGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTc                                          >  1:1189131/1‑150 (MQ=255)
                                                                                                            tataGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTc                                          >  1:705520/1‑150 (MQ=255)
                                                                                                               aGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAg                                       <  1:1571139/150‑1 (MQ=255)
                                                                                                                   gtgAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTAc                                                                                                                               <  1:750565/58‑1 (MQ=255)
                                                                                                                   gtgAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTAc                                                                                                                               >  2:750565/1‑58 (MQ=255)
                                                                                                                      aaGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAg                                       >  1:1567393/1‑143 (MQ=255)
                                                                                                                      aaGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAg                                       <  2:1567393/143‑1 (MQ=255)
                                                                                                                       agcagcACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGaa                               <  2:2071490/150‑1 (MQ=255)
                                                                                                                           gcaCGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTa                           >  2:309746/1‑150 (MQ=255)
                                                                                                                              cGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCa                                                       <  1:2266057/119‑1 (MQ=255)
                                                                                                                              cGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCa                                                       >  2:2266057/1‑119 (MQ=255)
                                                                                                                                taTAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCTTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGt                      <  1:1362985/149‑1 (MQ=255)
                                                                                                                                  tAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGa                    <  1:1255089/150‑1 (MQ=255)
                                                                                                                                     tATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttg                 >  2:1793427/1‑150 (MQ=255)
                                                                                                                                      aTCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAATTATTCATTGAGCGAGATGGCGaaa                                                                                                  <  1:1863539/68‑1 (MQ=255)
                                                                                                                                      aTCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGaaa                                                                                                  >  2:1863539/1‑68 (MQ=255)
                                                                                                                                      aTCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGt                                                          >  1:2028284/1‑108 (MQ=255)
                                                                                                                                      aTCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGt                                                          <  2:2028284/108‑1 (MQ=255)
                                                                                                                                      aTCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttgt                >  1:2092176/1‑150 (MQ=255)
                                                                                                                                       ctattgcttgcaTAAGCAAGCATTAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAg                                       >  1:68488/5‑126 (MQ=255)
                                                                                                                                       ctattgcttgcaTAAGCAAGCATTAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAg                                       <  2:68488/122‑1 (MQ=255)
                                                                                                                                        tattgcttgcaTAAGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttgttt              >  2:1573742/4‑150 (MQ=255)
                                                                                                                                          aTGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttgtttgt            <  2:1189131/150‑1 (MQ=255)
                                                                                                                                            gcttgcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTttgttg          <  2:499909/150‑1 (MQ=255)
                                                                                                                                             cttgcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTtgttgt         <  2:919957/150‑1 (MQ=255)
                                                                                                                                             cttgcatAAGCAAGCATTAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTAc                                                                           >  2:1820100/1‑84 (MQ=255)
                                                                                                                                             cttgcatAAGCAAGCATTAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTAc                                                                           <  1:1820100/84‑1 (MQ=255)
                                                                                                                                               tgcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTgtgt       <  2:2363666/150‑1 (MQ=255)
                                                                                                                                                   tATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTGTGTCTCa   <  2:2291152/150‑1 (MQ=255)
                                                                                                                                                    aTGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTGTGTCTCAg  <  1:1281934/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTGACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTGTGTCTCAG  >  NC_002947/3519978‑3520275

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