Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,851,001 A→C 5.7% 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 5.7% 197.7 / 3.9 70N147H (AAT→CAT) PP_3398curli fiber surface‑exposed nucleator CsgB
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base C (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

CATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAG  >  NC_002947/3850853‑3851142
                                                                                                                                                    |                                                                                                                                             
catcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCaa                                                                                                                                              <  2:1150609/150‑1 (MQ=255)
 atcaGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAAt                                                                                                                                             >  1:963943/1‑150 (MQ=255)
          cAGGGGAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGTCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCAt                                                                                                                                    <  2:45343/150‑1 (MQ=255)
          cAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCAt                                                                                                                                    >  2:469868/1‑150 (MQ=255)
            gggAGATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCATGCTGGCCATTGCCAGCCTgg                                                                                                                                  <  1:165905/150‑2 (MQ=255)
             ggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGCCTgg                                                                                                                                  <  1:536098/149‑2 (MQ=255)
             ggAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCCATTGCCAGCCTgg                                                                                                                                  >  2:536098/1‑148 (MQ=255)
              gAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTCGCAATGGCCAGCATGTTc                                                                                                                                >  2:603530/1‑150 (MQ=255)
              gAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCAAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTc                                                                                                                                >  1:1234056/1‑150 (MQ=255)
                  tGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGACAGCATCGTGCAGTCCGGCTAAGTTGTCAGCATCTCCGTGACCCAGGCTGGCAATGGCCAACATGTTCAGAt                                                                                                                            <  1:349627/150‑1 (MQ=255)
                      aaCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                        <  2:963943/150‑1 (MQ=255)
                      aaCAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCAcc                                                                                                                        >  2:1205518/1‑150 (MQ=255)
                             cGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGACAGCATGTTCAGATCACCCAATAcc                                                                                                                 <  2:79363/150‑1 (MQ=255)
                               tCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGTCAATGGCCAGCATGTTCAGATCACCCAATACCGa                                                                                                               <  1:942659/150‑1 (MQ=255)
                               tCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGa                                                                                                               >  1:459590/1‑150 (MQ=255)
                                    tCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGAc                                                                                                          <  1:789171/150‑1 (MQ=255)
                                      gCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGAc                                                                                                          <  2:666416/148‑1 (MQ=255)
                                      gCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGAc                                                                                                          >  1:666416/1‑148 (MQ=255)
                                       ccAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTg                                                                                                       >  1:624287/1‑150 (MQ=255)
                                        cAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTgg                                                                                                      <  2:846597/150‑1 (MQ=255)
                                        cAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTgg                                                                                                      >  1:1220635/1‑150 (MQ=255)
                                         agagCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGa                                                                                                     >  2:876131/1‑150 (MQ=255)
                                            gCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGc                                                                                                                                          <  2:1145508/110‑1 (MQ=255)
                                            gCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGc                                                                                                                                          >  1:1145508/1‑110 (MQ=255)
                                                 agcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGcaccac                                                                                             <  2:624287/150‑1 (MQ=255)
                                                  gcagcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCaccacc                                                                                            >  2:616841/1‑150 (MQ=255)
                                                    agcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCAt                                                                                          <  1:603530/150‑1 (MQ=255)
                                                     gcaACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATg                                                                                         >  1:1278614/1‑150 (MQ=255)
                                                      caacaaCGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGt                                                                                        >  2:443114/1‑150 (MQ=255)
                                                        acaacGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTAc                                                                                      >  2:39355/1‑150 (MQ=255)
                                                                      gAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGtt                                                                        >  1:52764/1‑150 (MQ=255)
                                                                      gAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGtt                                                                        <  1:39355/150‑1 (MQ=255)
                                                                         cAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCCCCACCATGTACAAGCTCGCTCCGTTAAg                                                                     <  1:876131/150‑1 (MQ=255)
                                                                         cAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAg                                                                     <  2:1234056/150‑1 (MQ=255)
                                                                          aGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGt                                                                    <  1:1205518/150‑1 (MQ=255)
                                                                           gATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGAt                                                                                                              >  1:1229008/1‑107 (MQ=255)
                                                                           gATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGAt                                                                                                              <  2:1229008/107‑1 (MQ=255)
                                                                           gATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAACTCGCTCCGTTAAGTg                                                                   >  2:446584/1‑150 (MQ=255)
                                                                                     gACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGt                                                         >  2:1100585/1‑150 (MQ=255)
                                                                                               ccAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGGCGATATCGTCATGGCCCTc                                               >  2:209077/1‑150 (MQ=255)
                                                                                                cAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATc                                                           <  1:265405/137‑1 (MQ=255)
                                                                                                cAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATc                                                           >  2:265405/1‑137 (MQ=255)
                                                                                                     tCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGt                                         <  1:469868/150‑1 (MQ=255)
                                                                                                      cGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTc                                        <  2:52764/150‑1 (MQ=255)
                                                                                                        tGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCACTGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAg                                      >  1:1283793/1‑150 (MQ=255)
                                                                                                        tGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAg                                                                                                                              >  2:362781/1‑62 (MQ=255)
                                                                                                        tGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAg                                                                                                                              <  1:362781/62‑1 (MQ=255)
                                                                                                             tCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCAt                                 >  1:1222581/1‑150 (MQ=255)
                                                                                                               cGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATgg                               <  2:331236/150‑1 (MQ=255)
                                                                                                                   tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAg                                      >  1:694350/1‑139 (MQ=255)
                                                                                                                   tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAg                                      <  2:694350/139‑1 (MQ=255)
                                                                                                                   tCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGgctgc                           >  1:318192/1‑150 (MQ=255)
                                                                                                                      gggCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGa                        <  2:1220635/150‑1 (MQ=255)
                                                                                                                        gCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACa                      <  2:459590/150‑1 (MQ=255)
                                                                                                                        gCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACa                      >  1:835459/1‑150 (MQ=255)
                                                                                                                          tCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGc                    >  1:990873/1‑150 (MQ=255)
                                                                                                                          tCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGc                    >  1:1220661/1‑150 (MQ=255)
                                                                                                                             aCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                 >  2:534106/1‑150 (MQ=255)
                                                                                                                             aCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                 <  1:443114/150‑1 (MQ=255)
                                                                                                                             aCAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCAACATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                 >  2:475332/1‑150 (MQ=255)
                                                                                                                                 cTCCGTGACCCCGGCAGGCAATGGCCAGCCTGTTCAGATCACCCAATACCGATAGACCTGCAGGCACCACCATGTACAAGCTCGCTCCGTTATGTGCCGCTATCGGCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGTCAGGACCTCGt             >  2:479012/1‑150 (MQ=255)
                                                                                                                                 cTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                 >  1:457417/1‑146 (MQ=255)
                                                                                                                                 cTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCAcc                 <  2:457417/146‑1 (MQ=255)
                                                                                                                                  tCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTc            >  1:454957/1‑150 (MQ=255)
                                                                                                                                  tCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGGGCCGCTTTCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTc            <  2:314335/150‑1 (MQ=255)
                                                                                                                                   ccGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTcc           <  2:990873/150‑1 (MQ=255)
                                                                                                                                       tgaacatgCTGGCCATTGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGa       >  2:963397/8‑150 (MQ=255)
                                                                                                                                         cccAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAcc     <  2:1278614/150‑1 (MQ=255)
                                                                                                                                         cccAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAcc     <  2:1220661/150‑1 (MQ=255)
                                                                                                                                          ccAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGAccc    >  2:1212341/1‑150 (MQ=255)
                                                                                                                                            aGGCTGGCAATGGCCAGCATGTTCAGATCAAACAAAAACGATAAACCTGGAGGCAACACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGACCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAg  >  2:562662/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CATCAGCCAACAGGGAAATGGCAACAGTGCGTCGATTCGCCAGAGCGGCAGCAGCAACAACGCCGCCATCGAGCAGATCGGCAACGACAACAGCGCCAGCATCGTGCAGTCCGGCTCAGGGCTCAACAGCTCCGTGACCCAGGCTGGCAATGGCCAGCATGTTCAGATCACCCAATACCGATAGACCTGGAGGCACCACCATGTACAAGCTCGCTCCGTTAAGTGCCGCTATCGTCCTGGCCCTCGCAGGTCAGGCCATGGCTGCTGACAGCACCTCGTCCCAGACCCAG  >  NC_002947/3850853‑3851142

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: ‑