Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,851,004 G→T 10.8% 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,0040GT10.8% 266.2 / 18.0 93G148C (GGC→TGC) PP_3398curli fiber surface‑exposed nucleator CsgB
Reads supporting (aligned to +/- strand):  ref base G (39/43);  new base T (4/6);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

TCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGCAAG  >  NC_002947/3850855‑3851151
                                                                                                                                                     |                                                                                                                                                   
tcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATg                                                                                                                                                     <  1:102187/150‑1 (MQ=255)
tcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGACTCAACAGCTCCGTGACCCAGGCTGGCAATg                                                                                                                                                     <  2:2258266/150‑1 (MQ=255)
  aGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGGGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGc                                                                                                                                                   >  1:253219/1‑150 (MQ=255)
       acacGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCa                                                                                                                                              <  1:191155/150‑1 (MQ=255)
          gggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGCCTgg                                                                                                                                           <  1:415823/150‑2 (MQ=255)
          gggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGCCTgg                                                                                                                                           <  2:319166/150‑2 (MQ=255)
          gggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGCCTgg                                                                                                                                           <  2:302463/150‑2 (MQ=255)
           ggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACCGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGtt                                                                                                                                          >  1:151109/1‑150 (MQ=255)
            gAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTc                                                                                                                                         <  1:951112/150‑1 (MQ=255)
               aTGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGa                                                                                                                                      <  1:554400/150‑1 (MQ=255)
                  gCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGCCTgg                                                                                                                                           >  1:1616329/1‑141 (MQ=255)
                  gCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGCCTgg                                                                                                                                           <  2:1616329/142‑2 (MQ=255)
                   cAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAc                                                                                                                                  <  2:2297375/150‑1 (MQ=255)
                   cAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAc                                                                                                                                  >  1:1896900/1‑150 (MQ=255)
                     aCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAccc                                                                                                                                <  1:495278/150‑1 (MQ=255)
                      cAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCa                                                                                                                               >  1:1627797/1‑150 (MQ=255)
                         tGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATa                                                                                                                            <  2:151109/150‑1 (MQ=255)
                             tCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGa                                                                                                                        >  1:2255494/1‑150 (MQ=255)
                               gATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATa                                                                                                                      >  2:1944078/1‑150 (MQ=255)
                                   cGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGAcc                                                                                                                  <  2:1223879/150‑1 (MQ=255)
                                    gCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCt                                                                                                                 >  2:325397/1‑150 (MQ=255)
                                    gCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCt                                                                                                                 >  2:169752/1‑150 (MQ=255)
                                     ccAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCCCGGCGGGGCTCAACAGCTCCGTGACCCCGCCACGCGTCCACCAGCATGTTCAGATCACCCAATACCGATAGACCTg                                                                                                                >  1:2340405/1‑150 (MQ=255)
                                      cAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                                 <  2:812589/132‑1 (MQ=255)
                                      cAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                                 >  1:812589/1‑132 (MQ=255)
                                       agagCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGa                                                                                                              <  2:833128/150‑1 (MQ=255)
                                       agagCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGa                                                                                                              <  1:1513612/150‑1 (MQ=255)
                                         agCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAgg                                                                                                            <  2:188647/150‑1 (MQ=255)
                                               agcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGAAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGcaccac                                                                                                      <  1:1340976/150‑1 (MQ=255)
                                                gcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCaccacc                                                                                                     <  1:1194947/150‑1 (MQ=255)
                                                  agcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCAt                                                                                                   <  1:2004273/150‑1 (MQ=255)
                                                     aacaacGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAg                                                                                                                     >  1:1652261/1‑129 (MQ=255)
                                                     aacaacGCCGCCATCGAGCACATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAg                                                                                                                     <  2:1652261/129‑1 (MQ=255)
                                                        aacGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACaa                                                                                             >  1:831462/1‑150 (MQ=255)
                                                            ccgccATCGGGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgctc                                                                                         <  2:253219/150‑1 (MQ=255)
                                                             cgccATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTACAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgctcg                                                                                        >  1:1011722/1‑150 (MQ=255)
                                                                  tCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTGAGGTCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCg                                                                                   >  2:340955/1‑150 (MQ=255)
                                                                    gAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGcctgg                                                                                                                                           >  1:1331570/1‑88 (MQ=255)
                                                                    gAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGCCTgg                                                                                                                                           <  2:1331570/92‑2 (MQ=255)
                                                                    gAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGtt                                                                                 <  2:2340405/150‑1 (MQ=255)
                                                                          aTCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGc                                                                           >  1:708044/1‑150 (MQ=255)
                                                                           tCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGtt                                                                                 <  1:2219116/143‑1 (MQ=255)
                                                                           tCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGtt                                                                                 >  2:2219116/1‑143 (MQ=255)
                                                                              gCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCt                                                                       >  2:1161556/1‑150 (MQ=255)
                                                                              gCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCATTATGCACCACCATGTACAAGCTCGCTCCGTTACGTGCCgat                                                                       >  2:999789/1‑148 (MQ=255)
                                                                               cAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTa                                                                      >  1:552965/1‑150 (MQ=255)
                                                                                    acaacaGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTc                                                                 <  2:2087877/150‑1 (MQ=255)
                                                                                      aacaGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCt                                                               >  2:1394643/1‑150 (MQ=255)
                                                                                       acaGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGCGGAACCACGATGTGCAAGCTCGTTCCGTTGAGTGCCGCTATCGTCCTg                                                              <  2:1835780/150‑1 (MQ=255)
                                                                                        caGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTgg                                                             >  1:1067363/1‑150 (MQ=255)
                                                                                          tcgcCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCCCTCCGTTAAGTGCACCTATCGTCCtatcc                                                           >  2:410925/2‑146 (MQ=255)
                                                                                           cgcCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGccc                                                          >  2:231051/1‑150 (MQ=255)
                                                                                           cgcCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGccc                                                          <  1:416461/150‑1 (MQ=255)
                                                                                              cAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCg                                                       >  1:435575/1‑150 (MQ=255)
                                                                                                 cATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAg                                                    <  2:1896900/150‑1 (MQ=255)
                                                                                                   tCGTGCACTCCGGCTCAGGGCTCAACAGCACCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGt                                                  <  1:169752/150‑1 (MQ=255)
                                                                                                     gTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCa                                                <  1:1944078/150‑1 (MQ=255)
                                                                                                      tGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAg                                               >  2:1574705/1‑150 (MQ=255)
                                                                                                       gCAGTCCGGATCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAgg                                              <  2:925811/150‑1 (MQ=255)
                                                                                                          gTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCa                                           <  1:591521/150‑1 (MQ=255)
                                                                                                           tCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCAt                                          <  2:1780830/150‑1 (MQ=255)
                                                                                                             cGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATgg                                        <  2:708044/150‑1 (MQ=255)
                                                                                                              ggCTCAGGGCTAAACAGCTCCGTGACCCAGGCTGGCAATGGCAAGCATGTTCAGATCACCCAACACCACTACCCCAGTATGCACCACCATGGACAAGCTCGCTCCGTTAAGTGCCGCTAGCGTCCTGGCCCTCGCAGGTCAGGCCATGgc                                       <  1:900304/150‑1 (MQ=255)
                                                                                                                cTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGgctg                                     <  2:1067363/150‑1 (MQ=255)
                                                                                                                 tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTAt                                                                     >  1:1307141/1‑117 (MQ=255)
                                                                                                                 tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTAt                                                                     <  2:1307141/117‑1 (MQ=255)
                                                                                                                    gggCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgct                                                                                          >  2:1138933/1‑93 (MQ=255)
                                                                                                                    gggCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAgct                                                                                          <  1:1138933/93‑1 (MQ=255)
                                                                                                                       cTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAg                              >  2:444726/1‑150 (MQ=255)
                                                                                                                        tCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGc                             >  2:1981073/1‑150 (MQ=255)
                                                                                                                          aaCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCg                                                                   >  1:702411/1‑110 (MQ=255)
                                                                                                                          aaCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCg                                                                   <  2:702411/110‑1 (MQ=255)
                                                                                                                          aaCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAc                           >  2:458993/1‑150 (MQ=255)
                                                                                                                           aCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                          >  2:1743155/1‑150 (MQ=255)
                                                                                                                               cTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAg                                                    <  2:2063700/120‑1 (MQ=255)
                                                                                                                               cTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAg                                                    >  1:2063700/1‑120 (MQ=255)
                                                                                                                                tCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTc                     <  2:1627797/150‑1 (MQ=255)
                                                                                                                                tCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTc                     >  2:1356986/1‑150 (MQ=255)
                                                                                                                                 ccGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTcc                    <  1:1161556/150‑1 (MQ=255)
                                                                                                                                  cGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTccc                   >  1:1810297/1‑150 (MQ=255)
                                                                                                                                   gTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCa                  <  2:1810297/150‑1 (MQ=255)
                                                                                                                                       cccAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAcc              >  2:1579947/1‑150 (MQ=255)
                                                                                                                                        ccAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAccc             >  1:1248364/1‑150 (MQ=255)
                                                                                                                                         catgCTGGCCATTGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCa            >  2:1489307/4‑150 (MQ=255)
                                                                                                                                         cAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCa            <  2:2352328/150‑1 (MQ=255)
                                                                                                                                          atgCTGGCCATTGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCa            >  2:208058/3‑149 (MQ=255)
                                                                                                                                          atgCTGGCCATTGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCa            <  1:208058/147‑1 (MQ=255)
                                                                                                                                            gCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGa         >  1:1559830/1‑150 (MQ=255)
                                                                                                                                                gctctGGCCAGCATGTTCAGATCACCCAACACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGACATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGc     <  1:571437/146‑1 (MQ=255)
                                                                                                                                                 cAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCt                                                               <  1:1165145/91‑1 (MQ=255)
                                                                                                                                                 cAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCt                                                               >  2:1165145/1‑91 (MQ=255)
                                                                                                                                                  aaTGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTaa                                                                               >  2:983130/1‑74 (MQ=255)
                                                                                                                                                  aaTGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTaa                                                                               <  1:983130/74‑1 (MQ=255)
                                                                                                                                                   aTGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGCAAg  <  2:1248364/150‑1 (MQ=255)
                                                                                                                                                   aTGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGCAAg  <  1:458993/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAGGACGGCAAG  >  NC_002947/3850855‑3851151

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