Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,851,004 G→T 8.3% G148C (GGC→TGC)  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,851,0040GT 8.3% 238.2 / 7.6 81G148C (GGC→TGC) PP_3398curli fiber surface‑exposed nucleator CsgB
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base T (3/4);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGG  >  NC_002947/3850857‑3851147
                                                                                                                                                   |                                                                                                                                               
aGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATgac                                                                                                                                               <  2:1314504/150‑3 (MQ=255)
 gCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAGTTATGACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGcc                                                                                                                                              >  1:1732336/1‑150 (MQ=255)
  ccAACAGGGAAATGGCAAAAGAGCGTCGATTCGCCAGAGCGGCAGCCGCAACAACGCCGCCCTCGAGCAGATCGGCCACGTCAACAGCGCTAGCATCTTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCa                                                                                                                                             <  1:1945542/150‑1 (MQ=255)
    aaCAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAATAACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGc                                                                                                                                           >  1:1342320/1‑150 (MQ=255)
       aGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATg                                                                                                                                        >  1:1397296/1‑150 (MQ=255)
        gggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGCCTgg                                                                                                                                       <  1:2305717/150‑2 (MQ=255)
         ggAAATGGCAACAGTGCGTCGATTCGCTAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACCGCTCCGTGACCCAGGCTGGCAATGGACAGCTTGtt                                                                                                                                      <  2:1274510/150‑1 (MQ=255)
         ggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGtt                                                                                                                                      <  2:101307/150‑1 (MQ=255)
         ggAAATGGCAAAAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGCGCTAAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGtt                                                                                                                                      <  2:777869/150‑1 (MQ=255)
             aTGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGa                                                                                                                                  <  1:1593091/150‑1 (MQ=255)
             aTGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGa                                                                                                                                  <  1:1596120/150‑1 (MQ=255)
                  aaCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAAGAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                             >  2:1373113/1‑150 (MQ=255)
                  aaCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                             <  2:1597470/150‑1 (MQ=255)
                  aaCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                             <  2:1314730/150‑1 (MQ=255)
                   aCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAccc                                                                                                                            <  2:1531567/150‑1 (MQ=255)
                   aCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGATCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAccc                                                                                                                            <  2:1533085/150‑1 (MQ=255)
                        gCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATGGCCAGCATGTTCATATCACCCAATAc                                                                                                                       <  2:1971721/150‑1 (MQ=255)
                         cGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATAcc                                                                                                                      <  2:1936132/150‑1 (MQ=255)
                          gTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCg                                                                                                                     >  2:2254776/1‑150 (MQ=255)
                          gTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCg                                                                                                                     >  1:211799/1‑150 (MQ=255)
                           tCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGa                                                                                                                    <  2:1397296/150‑1 (MQ=255)
                               ttCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGa                                                                                                                <  1:2384196/150‑1 (MQ=255)
                                    cAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTgg                                                                                                           >  1:1655398/1‑150 (MQ=255)
                                     agagCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGa                                                                                                          >  2:2263246/1‑150 (MQ=255)
                                     agagCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTCTACCACGCGTGCAATTCCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGa                                                                                                          <  1:1004424/150‑1 (MQ=255)
                                           gcagcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCa                                                                                                                                    >  1:332177/1‑118 (MQ=255)
                                           gcagcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCa                                                                                                                                    <  2:332177/118‑1 (MQ=255)
                                             agcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGcaccac                                                                                                  <  1:1544655/150‑1 (MQ=255)
                                             agcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGcaccac                                                                                                  >  1:2445066/1‑150 (MQ=255)
                                              gcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATCGACATGGAGGCaccacc                                                                                                 <  2:1732336/150‑1 (MQ=255)
                                              gcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCaccacc                                                                                                 >  1:301656/1‑150 (MQ=255)
                                                agcaACAACGCCGCCATCGCAACGAGCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCAt                                                                                               <  2:1135179/150‑1 (MQ=255)
                                                agcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCAt                                                                                               <  2:2445066/150‑1 (MQ=255)
                                                   agcaacGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTa                                                                                            >  2:2190951/3‑150 (MQ=255)
                                                      aacGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACaa                                                                                         <  1:2263246/150‑1 (MQ=255)
                                                        cgccgcCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGTGGCACCACCATGTACAAgc                                                                                       <  2:1403688/150‑1 (MQ=255)
                                                          ccgccATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgctc                                                                                     >  2:1890238/1‑150 (MQ=255)
                                                           cgccATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGCCTgg                                                                                                                                       <  2:2275642/99‑2 (MQ=255)
                                                           cgccATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGACTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGcctgg                                                                                                                                       >  1:2275642/1‑95 (MQ=255)
                                                              cATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgctcgctc                                                                                 <  2:976062/150‑1 (MQ=255)
                                                               aTCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTcc                                                                                >  1:390627/1‑150 (MQ=255)
                                                                  gAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGACTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGtt                                                                             >  2:1325352/1‑150 (MQ=255)
                                                                           ggCAACGACAACAGCGACAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCACGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGc                                                                    <  2:1231055/150‑1 (MQ=255)
                                                                            gCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCt                                                                   >  1:2406197/1‑150 (MQ=255)
                                                                                  acaacaGCGCCAGCATCGTGCAGTCCGGCTCAGGACTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATAcc                                                                                                                      <  2:736968/93‑1 (MQ=255)
                                                                                  acaacaGCGCCAGCATCGTGCAGTCCGGCTCAGGACTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATAcc                                                                                                                      >  1:736968/1‑93 (MQ=255)
                                                                                    aacaGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCt                                                           >  2:2202081/1‑150 (MQ=255)
                                                                                        gcgcCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                             <  2:1316089/80‑1 (MQ=255)
                                                                                        gcgcCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                             >  1:1316089/1‑80 (MQ=255)
                                                                                         cgcCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGATCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGccc                                                      >  2:1377154/1‑150 (MQ=255)
                                                                                         cgcCAGCATCGTGCAGTCCGACTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGccc                                                      <  1:1325352/150‑1 (MQ=255)
                                                                                            cAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCt                                                                                                             <  1:1075179/92‑1 (MQ=255)
                                                                                            cAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCt                                                                                                             >  2:1075179/1‑92 (MQ=255)
                                                                                            cAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCg                                                   >  2:1993728/1‑150 (MQ=255)
                                                                                              gCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCa                                                 <  1:2190951/150‑1 (MQ=255)
                                                                                               cATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAg                                                <  1:1644656/150‑1 (MQ=255)
                                                                                                        gTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCa                                       >  1:2213983/1‑150 (MQ=255)
                                                                                                         tCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCAt                                      >  1:1508926/1‑150 (MQ=255)
                                                                                                              cTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                             <  2:1832407/58‑1 (MQ=255)
                                                                                                              cTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                             >  1:1832407/1‑58 (MQ=255)
                                                                                                              cTCAGGGATCAACAGCTCGATGAGCCAGCATTTGAATCGCCAAGATCCACAGATCACCCAATACGGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGgctg                                 >  2:1593968/1‑150 (MQ=255)
                                                                                                               tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGca                                                                                                      <  2:380029/80‑1 (MQ=255)
                                                                                                               tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGca                                                                                                      >  1:380029/1‑80 (MQ=255)
                                                                                                               tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGgctgc                                >  1:48120/1‑150 (MQ=255)
                                                                                                                  gggCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACACCCTCGCGCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGa                             <  2:808990/150‑1 (MQ=255)
                                                                                                                   ggCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGAc                            >  2:1875904/1‑150 (MQ=255)
                                                                                                                         aCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGa                                                                                                                                  >  1:1021199/1‑42 (MQ=255)
                                                                                                                         aCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGa                                                                                                                                  <  2:1021199/42‑1 (MQ=255)
                                                                                                                         aCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                      >  2:2037282/1‑150 (MQ=255)
                                                                                                                           aGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTg                                                                        >  2:2333259/1‑98 (MQ=255)
                                                                                                                           aGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTg                                                                        <  1:2333259/98‑1 (MQ=255)
                                                                                                                             cTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGt                  <  1:1373113/150‑1 (MQ=255)
                                                                                                                                 gTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                      >  2:2062065/1‑142 (MQ=255)
                                                                                                                                 gTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                      <  1:2062065/142‑1 (MQ=255)
                                                                                                                                   ggcccATGCTTGCAATGGCGAGCATGTTCAGATCCCCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGa            <  2:390627/148‑1 (MQ=255)
                                                                                                                                   gACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGa            <  2:2068797/150‑1 (MQ=255)
                                                                                                                                   gACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGa            <  1:1875904/150‑1 (MQ=255)
                                                                                                                                        atgCTGGCCATTGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAg       >  2:321503/3‑150 (MQ=255)
                                                                                                                                          gCTGGCCATTGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGa                                                                                                          >  1:1659820/1‑49 (MQ=255)
                                                                                                                                          gCTGGCCATTGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGa                                                                                                          <  2:1659820/49‑1 (MQ=255)
                                                                                                                                            tGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACg   >  1:2341522/1‑150 (MQ=255)
                                                                                                                                             ggCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACACCACCTCGTCCCAGACCCAGGACgg  <  2:211799/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGG  >  NC_002947/3850857‑3851147

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