Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,104,932 C→T 6.0% intergenic (+30/‑29) csgA → / → csgC curlin, major subunit/inhibitor of CsgA amyloid formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,104,9320CT 6.0% 520.9 / 14.0 165intergenic (+30/‑29)csgA/csgCcurlin, major subunit/inhibitor of CsgA amyloid formation
Reads supporting (aligned to +/- strand):  ref base C (67/88);  new base T (4/6);  total (71/94)
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.48e-01

Only 100 of 168 total aligned reads displayed.

TGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGA  >  NC_000913/1104784‑1105080
                                                                                                                                                    |                                                                                                                                                     
tGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGTGCa                                                                                                                                                      <  2:37383/150‑1 (MQ=255)
   cGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGc                                                                                                                                                   >  1:147169/1‑150 (MQ=255)
      ttAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCt                                                                                                                                                >  2:315834/1‑150 (MQ=255)
        aaaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGt                                                                                                                                             <  1:1196553/150‑1 (MQ=255)
        aaaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGt                                                                                                                                             <  2:1896316/150‑1 (MQ=255)
              ttCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAGACAGGGCGCAAGCCCTGttttttt                                                                                                                                       <  1:1928155/150‑1 (MQ=255)
              tggtgtggtgGCTAGTGTGCTGCAGTTGACCAGACTGCATCTACCTCCTCCGTCAACGTGAGTGAGGTTGGCTTTGGTAACAACGCGACCTCTCATGAGTACTAATACATGATATGTATTACATAAACATGGCGCTATCCCTGttttttt                                                                                                                                       <  2:1014107/146‑1 (MQ=255)
                 ggtggtGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCg                                                                                                                                    >  1:1081161/1‑150 (MQ=255)
                  gtggtgGCAACGGTGGTGCATTTGACCATACTTCATCTAACTCCTTCGTCAACGTGACTCAGGTTGGCTTTTGTAAAAACGCGACCGGTCATTATTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCgg                                                                                                                                   <  2:1859539/150‑1 (MQ=255)
                  gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCgg                                                                                                                                   >  2:2035215/1‑150 (MQ=255)
                  gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCgg                                                                                                                                   <  1:595731/150‑1 (MQ=255)
                  gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCgg                                                                                                                                   <  1:1455629/150‑1 (MQ=255)
                  gtgggcgCAACTGTTCTGCAGTTGACCAGACTGCATCTTACTTATCTGTCAACGTGACTAAGGTTTTCTTTTGTAACATCGAGACCTATCATCATTAATAATATATCATTTTTATTAAATAAACAGGGCGCAAGCCCTGTTTTTTTTCgg                                                                                                                                   <  2:494261/144‑1 (MQ=255)
                  ggggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCgg                                                                                                                                   <  1:746321/148‑1 (MQ=255)
                     gtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTATTTTTCGGgag                                                                                                                                <  2:1923247/150‑1 (MQ=255)
                      tggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGgaga                                                                                                                               >  1:2207036/1‑150 (MQ=255)
                              gTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATg                                                                                                                       <  2:433676/150‑1 (MQ=255)
                               tgctgcAGTTGTCCAGTCTGCTTCTAAGTCCTCCTTCAACGTTTCTCAGGTTGGATTTGGTATTAACACGACCGCTTATCAGTAATAATACATCATTTTTATTTGAGTAATAGGGCGCAAGCGCTTTTTTTTTTCGGGAGAAGAATATgt                                                                                                                      <  2:235035/150‑2 (MQ=255)
                               tgctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGa                                                                                                                      <  1:1546919/150‑1 (MQ=255)
                                 ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCAATTCTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTAGGGAGAAGAATATGAAt                                                                                                                    >  1:905387/1‑150 (MQ=255)
                                  tgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATa                                                                                                                   <  1:2186636/150‑1 (MQ=255)
                                  tgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATa                                                                                                                   <  1:2161922/150‑1 (MQ=255)
                                    cAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACg                                                                                                                 <  1:1354071/150‑1 (MQ=255)
                                    cAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACg                                                                                                                 <  2:1121251/150‑1 (MQ=255)
                                      gTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGtt                                                                                                               >  1:1991157/1‑150 (MQ=255)
                                         gACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGttatt                                                                                                            <  2:1730775/150‑1 (MQ=255)
                                          aCCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGttatta                                                                                                           <  2:141653/150‑1 (MQ=255)
                                           ccAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTAc                                                                                                          >  2:1249882/1‑150 (MQ=255)
                                           ccAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTAc                                                                                                          <  1:431126/150‑1 (MQ=255)
                                              gACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTcc                                                                                                       >  1:1842038/1‑150 (MQ=255)
                                               aCTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCt                                                                                                      <  2:1081161/150‑1 (MQ=255)
                                               aCTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCt                                                                                                      >  2:694773/1‑150 (MQ=255)
                                                cTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCtt                                                                                                     >  2:649556/1‑150 (MQ=255)
                                                  gCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCg                                                                                                  <  2:988467/150‑1 (MQ=255)
                                                   cATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCg                                                                                                  >  1:2274787/1‑150 (MQ=255)
                                                    aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                                 <  1:2179623/150‑1 (MQ=255)
                                                    aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                                 >  2:788130/1‑150 (MQ=255)
                                                     tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                <  1:2411272/150‑1 (MQ=255)
                                                      cTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATTAATACGTTATTACTCCTTGCggga                                                                                               <  1:100294/150‑3 (MQ=255)
                                                       tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                              <  1:1876050/150‑1 (MQ=255)
                                                         aCTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACtt                                                                                            <  1:794933/150‑1 (MQ=255)
                                                            cctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTcc                                                                                         <  1:2035215/150‑1 (MQ=255)
                                                              tccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAg                                                                                       <  1:1965164/150‑1 (MQ=255)
                                                              tccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAg                                                                                       <  2:564383/150‑1 (MQ=255)
                                                               ccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGt                                                                                      <  1:694773/150‑1 (MQ=255)
                                                                cGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGt                                                                                      >  2:1664569/1‑150 (MQ=255)
                                                                 gTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCa                                                                                    >  2:592632/1‑150 (MQ=255)
                                                                     aCGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATa                                                                                <  1:872893/150‑1 (MQ=255)
                                                                            tCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTa                                                                         <  1:1068374/150‑1 (MQ=255)
                                                                             cAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTaa                                                                        <  2:2445857/150‑1 (MQ=255)
                                                                                 ttGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACg                                                                    >  1:1333694/1‑150 (MQ=255)
                                                                                   ggCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAc                                                                  <  1:2401713/150‑1 (MQ=255)
                                                                                   ggCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAc                                                                  <  2:2207036/150‑1 (MQ=255)
                                                                                   ggCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATAAGAc                                                                  >  2:2214711/1‑150 (MQ=255)
                                                                                      ttttGTTACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCa                                                               <  1:2074150/150‑1 (MQ=255)
                                                                                      tttGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCa                                                               <  2:1568162/150‑1 (MQ=255)
                                                                                       ttGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAg                                                              <  2:510992/150‑1 (MQ=255)
                                                                                         ggTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCa                                                            <  1:611950/150‑1 (MQ=255)
                                                                                         ggTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATAAGACCCAGCa                                                            >  1:1948304/1‑150 (MQ=255)
                                                                                         ggTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTACTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCa                                                            >  2:1163919/1‑150 (MQ=255)
                                                                                               aacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                      <  2:651179/150‑1 (MQ=255)
                                                                                                 cgcgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAc                                                                  <  2:2256267/136‑1 (MQ=255)
                                                                                                 cgcgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAc                                                                  >  1:2256267/1‑136 (MQ=255)
                                                                                                 cgcgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtg                                                    <  2:1991157/150‑1 (MQ=255)
                                                                                                 cgcgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtg                                                    <  2:1681756/150‑1 (MQ=255)
                                                                                                     aCCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGTGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGta                                                <  1:1213332/150‑1 (MQ=255)
                                                                                                     aCCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGta                                                <  2:530274/150‑1 (MQ=255)
                                                                                                       cGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGtata                                              <  1:997559/150‑1 (MQ=255)
                                                                                                       cGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGtata                                              <  2:959326/150‑1 (MQ=255)
                                                                                                         cTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATAcc                                            <  1:1452981/150‑1 (MQ=255)
                                                                                                           catcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCat                                          >  2:18133/1‑150 (MQ=255)
                                                                                                              caGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCattat                                       >  2:1228548/1‑150 (MQ=255)
                                                                                                               aGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCattatt                                      >  1:1679465/1‑150 (MQ=255)
                                                                                                                gTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                     >  1:343086/1‑150 (MQ=255)
                                                                                                                   cTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTg                                  >  1:1677988/1‑150 (MQ=255)
                                                                                                                   cTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTg                                  >  2:581687/1‑150 (MQ=255)
                                                                                                                   cTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTg                                  <  2:2287018/150‑1 (MQ=255)
                                                                                                                   cTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTg                                  >  1:2067507/1‑150 (MQ=255)
                                                                                                                          atcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCAc                           >  1:1076869/1‑150 (MQ=255)
                                                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCAct                          <  2:2274787/150‑1 (MQ=255)
                                                                                                                            catTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCActc                         >  1:563869/1‑150 (MQ=255)
                                                                                                                                  tATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTc                   <  1:581328/150‑1 (MQ=255)
                                                                                                                                  tATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTc                   >  2:1268371/1‑150 (MQ=255)
                                                                                                                                   cgaaaaaaaaCAGGGCTTGCGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCa                  >  1:2316773/8‑150 (MQ=255)
                                                                                                                                   aTTACATAAACACGGCGCAAGCCCTGTTTTTTTTCGGGAGAATAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCa                  <  2:2472347/150‑1 (MQ=255)
                                                                                                                                   aTTACAGAAACAGGGAGCAAGCCCTGGTTTTTTTCGGGAGAAGAATATGAATACGGTATTACGCCTTGCGGCACTTTCCAGTCAGAGAACCTTTAATACTCCCCAGCAAGGGGATGTGTCTACCATTATTCCTGAAGTCACTCGTaagca                  >  2:1405629/1‑146 (MQ=255)
                                                                                                                                    ttACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAGGTCACTCTTACTCaa                 >  2:722304/1‑150 (MQ=255)
                                                                                                                                     tACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGa                                 >  1:463367/1‑133 (MQ=255)
                                                                                                                                     tACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAAt                >  2:1791150/1‑150 (MQ=255)
                                                                                                                                      aCAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc               <  2:1598476/150‑1 (MQ=255)
                                                                                                                                       cAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCt              >  1:1268170/1‑150 (MQ=255)
                                                                                                                                       cAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCt              <  2:2380565/150‑1 (MQ=255)
                                                                                                                                        aaaaaCAGGGCTTGCGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCtt             >  1:883432/3‑150 (MQ=255)
                                                                                                                                         gAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTg            <  2:1719270/150‑1 (MQ=255)
                                                                                                                                           aaCAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc          <  1:1228548/150‑1 (MQ=255)
                                                                                                                                               gggCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg      >  1:1403883/1‑150 (MQ=255)
                                                                                                                                                ggCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc     >  1:69951/1‑150 (MQ=255)
                                                                                                                                                  cgcAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCag   <  2:147169/150‑1 (MQ=255)
                                                                                                                                                  cgcAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCag   <  1:169538/150‑1 (MQ=255)
                                                                                                                                                   gcAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCaga  >  2:1195537/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGA  >  NC_000913/1104784‑1105080

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