Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,852,923 A→C 6.3% N147H (AAT→CAT)  PP_3398 → curli fiber surface‑exposed nucleator CsgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,852,9230AC 6.3% 464.7 / 16.6 142N147H (AAT→CAT) PP_3398curli fiber surface‑exposed nucleator CsgB
Reads supporting (aligned to +/- strand):  ref base A (61/72);  new base C (7/2);  total (68/74)
Fisher's exact test for biased strand distribution p-value = 8.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 148 total aligned reads displayed.

TCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGCAA  >  NC_002947/3852774‑3853072
                                                                                                                                                     |                                                                                                                                                     
  atcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAAt                                                                                                                                                     <  1:3395227/150‑1 (MQ=255)
   tcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATg                                                                                                                                                    >  1:4079270/1‑150 (MQ=255)
   tcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATg                                                                                                                                                    <  1:4282407/150‑1 (MQ=255)
   tcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATg                                                                                                                                                    <  1:4843521/150‑1 (MQ=255)
   tcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATg                                                                                                                                                    <  1:3440814/150‑1 (MQ=255)
         aaCAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGc                                                                                                                                              >  1:679139/1‑150 (MQ=255)
         aaCAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGc                                                                                                                                              >  2:1540264/1‑150 (MQ=255)
         aaCAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGc                                                                                                                                              >  2:1540301/1‑150 (MQ=255)
         aaCAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGc                                                                                                                                              >  2:4985976/1‑150 (MQ=255)
          aCAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCa                                                                                                                                             >  2:783112/1‑150 (MQ=255)
            aGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATg                                                                                                                                           >  1:1812511/1‑150 (MQ=255)
              ggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGtt                                                                                                                                         >  2:3547387/1‑150 (MQ=255)
              ggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGtt                                                                                                                                         <  2:2920482/150‑1 (MQ=255)
              ggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGtt                                                                                                                                         <  2:4074750/150‑1 (MQ=255)
               gAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTc                                                                                                                                        >  1:2407550/1‑150 (MQ=255)
               gAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTc                                                                                                                                        <  1:2882446/150‑1 (MQ=255)
                  aTGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACACCGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGa                                                                                                                                     <  2:3928412/150‑1 (MQ=255)
                  aTGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGa                                                                                                                                     <  1:2109247/150‑1 (MQ=255)
                   tGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGAt                                                                                                                                    <  2:3025349/150‑1 (MQ=255)
                   tGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGAt                                                                                                                                    <  2:1185018/150‑1 (MQ=255)
                      cAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAc                                                                                                                                 <  1:1879937/150‑1 (MQ=255)
                       aaCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                                >  1:308311/1‑150 (MQ=255)
                        aCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAccc                                                                                                                               <  1:487690/150‑1 (MQ=255)
                               gTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCg                                                                                                                        >  2:502986/1‑150 (MQ=255)
                                  gATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATa                                                                                                                     <  1:2957608/150‑1 (MQ=255)
                                  gATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATa                                                                                                                     <  1:3278137/150‑1 (MQ=255)
                                   aTTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAg                                                                                                                    <  2:462019/150‑1 (MQ=255)
                                    ttCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGa                                                                                                                   >  1:1116671/1‑150 (MQ=255)
                                       gCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCt                                                                                                                <  2:986524/150‑1 (MQ=255)
                                       gCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCt                                                                                                                <  2:986409/150‑1 (MQ=255)
                                        ccAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTg                                                                                                               >  1:972076/1‑150 (MQ=255)
                                          agagCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCGACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGa                                                                                                             <  2:549376/150‑1 (MQ=255)
                                          agagCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGa                                                                                                             <  1:1804638/150‑1 (MQ=255)
                                           gagCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAg                                                                                                            <  2:1628476/150‑1 (MQ=255)
                                             gCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATCCCGATAGACCTGGAGGc                                                                                                          <  2:5114976/150‑1 (MQ=255)
                                                  agcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGcaccac                                                                                                     <  2:3790825/150‑1 (MQ=255)
                                                   gcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGcaccagc                                                                                                    <  1:2533642/150‑3 (MQ=255)
                                                   gcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCaccacc                                                                                                    >  2:847175/1‑150 (MQ=255)
                                                     agcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCAt                                                                                                  >  1:546314/1‑150 (MQ=255)
                                                      gcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATg                                                                                                 <  1:3515544/150‑1 (MQ=255)
                                                        aacaacGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTa                                                                                               >  2:1406424/1‑150 (MQ=255)
                                                         acaacGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTAc                                                                                              <  2:2532745/150‑1 (MQ=255)
                                                           aacGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACaa                                                                                            <  1:228792/150‑1 (MQ=255)
                                                           aacGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACaa                                                                                            >  1:3338247/1‑150 (MQ=255)
                                                             cgccgcCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgc                                                                                          >  1:2445422/1‑150 (MQ=255)
                                                              gccgccATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgct                                                                                         <  1:4170378/150‑1 (MQ=255)
                                                               ccgccATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgctc                                                                                        <  1:141362/150‑1 (MQ=255)
                                                               ccgccATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgctc                                                                                        <  1:4640456/150‑1 (MQ=255)
                                                                    aTCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGGCCTGGAGGCACCACCATGTACAAGCTCGCTcc                                                                                   <  2:4342421/150‑1 (MQ=255)
                                                                      cGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGt                                                                                 <  1:314452/150‑1 (MQ=255)
                                                                       gAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGtt                                                                                >  2:2247276/1‑150 (MQ=255)
                                                                         gCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTaa                                                                              >  1:5088583/1‑150 (MQ=255)
                                                                         gCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTaa                                                                              >  1:899401/1‑150 (MQ=255)
                                                                         gCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTaa                                                                              <  2:1087716/150‑1 (MQ=255)
                                                                           aGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGt                                                                            <  1:1139867/150‑1 (MQ=255)
                                                                                 gCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCt                                                                      >  2:1911878/1‑150 (MQ=255)
                                                                                 gCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCt                                                                      >  2:5115917/1‑150 (MQ=255)
                                                                                    aCGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATc                                                                   >  2:4509181/1‑150 (MQ=255)
                                                                                    aCGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATc                                                                   <  2:1564555/150‑1 (MQ=255)
                                                                                         aacaGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCt                                                              <  1:2276502/150‑1 (MQ=255)
                                                                                         aacaGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCt                                                              >  2:357217/1‑150 (MQ=255)
                                                                                         aacaGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCt                                                              <  1:4754596/150‑1 (MQ=255)
                                                                                           caGCGCCAGCATCGTGCAGTCCGGCTCATGGCTCAACAGCTACGTGACCCATTCTGGCAATGGCCAGAATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACACGCTCGCTCCGTTAAGTGCCGCTATCGTCCTgg                                                            <  1:2058478/150‑1 (MQ=255)
                                                                                           caGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTgg                                                            <  1:4886773/150‑1 (MQ=255)
                                                                                                 cAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCg                                                      >  2:1266677/1‑150 (MQ=255)
                                                                                                    cATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAg                                                   <  1:1102276/150‑1 (MQ=255)
                                                                                                    cATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAg                                                   >  2:2078812/1‑150 (MQ=255)
                                                                                                        gTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCa                                               >  1:1900770/1‑150 (MQ=255)
                                                                                                         tGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAg                                              >  1:1631393/1‑150 (MQ=255)
                                                                                                         tGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAg                                              >  1:1631453/1‑150 (MQ=255)
                                                                                                          gCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAgg                                             <  2:624573/150‑1 (MQ=255)
                                                                                                                    tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGgctgc                                   >  2:1917280/1‑150 (MQ=255)
                                                                                                                    tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGgctgc                                   <  1:1479968/150‑1 (MQ=255)
                                                                                                                    tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGgctgc                                   <  1:355697/150‑1 (MQ=255)
                                                                                                                       gggCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGa                                <  1:3745411/150‑1 (MQ=255)
                                                                                                                       gggCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGa                                <  1:3192509/150‑1 (MQ=255)
                                                                                                                         gCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACa                              <  2:1812511/150‑1 (MQ=255)
                                                                                                                         gCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACa                              >  2:3169929/1‑150 (MQ=255)
                                                                                                                           tCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGt                                                 >  1:363705/1‑129 (MQ=255)
                                                                                                                           tCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGc                            >  2:4081225/1‑150 (MQ=255)
                                                                                                                             aaCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGTAGGCACCACCATGTACAAGCTCGCTCCGTGAAGTGCCGCTAACGTCCAGGCCCTCGCAGGTCAGGCCATGGCTCCTGACAGCAc                          >  2:3828061/1‑150 (MQ=255)
                                                                                                                              aCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                         >  1:1023273/1‑150 (MQ=255)
                                                                                                                              aCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                         <  1:683320/150‑1 (MQ=255)
                                                                                                                               cAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCt                        >  1:2690448/1‑150 (MQ=255)
                                                                                                                                  cTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGt                     >  2:3366021/1‑150 (MQ=255)
                                                                                                                                  cTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGt                     >  2:1570946/1‑150 (MQ=255)
                                                                                                                                       tGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAg                >  2:4224382/1‑150 (MQ=255)
                                                                                                                                        tgaacatgCTGGCCATTGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGa               >  2:1657864/8‑150 (MQ=255)
                                                                                                                                          aacatgCTGGCCATTGCCAGCCTGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAcc             >  1:3056726/6‑150 (MQ=255)
                                                                                                                                          aacatgCTGGCCATTGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAcc             >  1:2850007/6‑150 (MQ=255)
                                                                                                                                           ccAGGCTGGCACTGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAccc            >  2:2562895/1‑150 (MQ=255)
                                                                                                                                           ccAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAccc            >  1:2993282/1‑150 (MQ=255)
                                                                                                                                           ccAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAccc            >  1:436103/1‑150 (MQ=255)
                                                                                                                                           acatgCTGGCCATTGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAccc            >  1:2525360/5‑150 (MQ=255)
                                                                                                                                                 tGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACg      <  1:1540264/150‑1 (MQ=255)
                                                                                                                                                  ggCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACgg     <  2:3300243/150‑1 (MQ=255)
                                                                                                                                                    cAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGCa   >  2:182356/1‑150 (MQ=255)
                                                                                                                                                     aaTGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGCaa  <  1:3826894/150‑1 (MQ=255)
                                                                                                                                                     aaTGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGCaa  <  2:4613662/150‑1 (MQ=255)
                                                                                                                                                     aaTGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGCaa  <  1:3043257/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGCAA  >  NC_002947/3852774‑3853072

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