Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,851,001 A→C 7.6% 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,851,0010AC 7.6% 239.2 / 9.7 87N147H (AAT→CAT) PP_3398curli fiber surface‑exposed nucleator CsgB
Reads supporting (aligned to +/- strand):  ref base A (37/43);  new base C (3/4);  total (40/47)
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.39e-01

ATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGG  >  NC_002947/3850854‑3851147
                                                                                                                                                   |                                                                                                                                                  
ttcaGCCAACAGGGAAATGGCAACAGGGCGTCGATTCGCCAGAGCGGCACCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCCTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAAt                                                                                                                                                  <  2:1226879/149‑1 (MQ=255)
atcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAAt                                                                                                                                                  >  1:101307/1‑150 (MQ=255)
atcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAAt                                                                                                                                                  >  1:1936132/1‑150 (MQ=255)
   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:1596120/150‑1 (MQ=255)
                aTGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGa                                                                                                                                  <  1:1593091/150‑1 (MQ=255)
                   gCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAAt                                                                                                                                                  <  1:1653985/131‑1 (MQ=255)
                   gCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAAt                                                                                                                                                  >  2:1653985/1‑131 (MQ=255)
                     aaCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAAGAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                             >  2:1373113/1‑150 (MQ=255)
                     aaCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                             <  2:1314730/150‑1 (MQ=255)
                     aaCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                             <  2:1597470/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                                                                                                                                    <  2:332177/118‑1 (MQ=255)
                                              gcagcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCa                                                                                                                                    >  1:332177/1‑118 (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                                                                                                      >  1:380029/1‑80 (MQ=255)
                                                                                                                  tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGca                                                                                                      <  2:380029/80‑1 (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                                                                        <  1:2333259/98‑1 (MQ=255)
                                                                                                                              aGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTg                                                                        >  2:2333259/1‑98 (MQ=255)
                                                                                                                                cTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGt                  <  1:1373113/150‑1 (MQ=255)
                                                                                                                                    gTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                      <  1:2062065/142‑1 (MQ=255)
                                                                                                                                    gTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                      >  2:2062065/1‑142 (MQ=255)
                                                                                                                                      ggcccATGCTTGCAATGGCGAGCATGTTCAGATCCCCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGa            <  2:390627/148‑1 (MQ=255)
                                                                                                                                      gACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGa            <  1:1875904/150‑1 (MQ=255)
                                                                                                                                      gACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGa            <  2:2068797/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)
                                                                                                                                                   |                                                                                                                                                  
ATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGG  >  NC_002947/3850854‑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: