Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,104,931 G→T 5.7% intergenic (+29/‑30) csgA → / → csgC curlin subunit, amyloid curli fibers, cryptic/curli assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,104,9310GT 5.7% 395.8 / 6.1 134intergenic (+29/‑30)csgA/csgCcurlin subunit, amyloid curli fibers, cryptic/curli assembly protein
Reads supporting (aligned to +/- strand):  ref base G (58/67);  new base T (3/6);  total (61/73)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

ATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGA  >  NC_000913/1104783‑1105080
                                                                                                                                                    |                                                                                                                                                      
aTGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGgcgc                                                                                                                                                       >  1:1528429/1‑150 (MQ=255)
    cGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGc                                                                                                                                                   <  1:916084/150‑1 (MQ=255)
      gTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGccc                                                                                                                                                 <  1:1823969/150‑1 (MQ=255)
          aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑tt                                                                                                                                            >  2:169837/1‑150 (MQ=255)
          aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑tt                                                                                                                                            >  2:56974/1‑150 (MQ=255)
          aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGACCTG‑tt                                                                                                                                            >  1:658219/1‑150 (MQ=255)
           aCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑ttt                                                                                                                                           >  1:321648/1‑150 (MQ=255)
            cAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACTGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttc                                                                                                                                          <  1:391189/150‑2 (MQ=255)
            cAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttc                                                                                                                                          <  2:975380/150‑2 (MQ=255)
            cAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCAGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttc                                                                                                                                          <  1:459792/150‑2 (MQ=255)
              gTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttctg                                                                                                                                        <  2:1342194/150‑4 (MQ=255)
               ttCGGTGGTTGCAAAGTTGCTTCAGTTTAACAGACTTCATATAAATACTCCTTAAACGTGACTCATTTTGGCTTTGGTAACAAATCGACATATCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTCAAGCCCTG‑ttttttt                                                                                                                                       <  2:1288252/150‑1 (MQ=255)
               ttCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑ttttttt                                                                                                                                       <  1:64461/150‑1 (MQ=255)
                 gggtgggGGCAACGGTGCTGCATTTTACCAGACTGAATCTAACTCCTCAGTCAAAGTGACGCAGGTTGGCGTTGGTAACAACGCGACCGCTCATCAGTACTAATACAGCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTc                                                                                                                                     <  2:1003612/149‑1 (MQ=255)
                 cggcggTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTc                                                                                                                                     <  2:1175579/150‑1 (MQ=255)
                  ggtggtGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTGTGGTAACAACGCGACCGCTCATCAGTACTTATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCg                                                                                                                                    <  1:797884/150‑1 (MQ=255)
                  gggggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCg                                                                                                                                    <  1:618987/147‑1 (MQ=255)
                   gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGGTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAACCAGGGCGCAAGCCCTG‑TTTTTTTTCgg                                                                                                                                   <  2:249228/150‑1 (MQ=255)
                   gtggtgGCAAAGGTGCTGCAGGTGTCCAGCCTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTTGTAACAACGCGACCGCGCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCgg                                                                                                                                   <  2:49675/150‑1 (MQ=255)
                    tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACTACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCggg                                                                                                                                  <  2:1548391/150‑1 (MQ=255)
                    tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCggg                                                                                                                                  <  2:1488725/150‑1 (MQ=255)
                    tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCggg                                                                                                                                  <  2:1189460/150‑1 (MQ=255)
                      gtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGgag                                                                                                                                <  2:405638/150‑1 (MQ=255)
                      gtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGgag                                                                                                                                >  1:874788/1‑150 (MQ=255)
                          cAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGagaaga                                                                                                                            <  2:1220223/150‑1 (MQ=255)
                                  ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                    >  2:709457/1‑150 (MQ=255)
                                  ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                    >  1:1395053/1‑150 (MQ=255)
                                  ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                    >  2:892825/1‑150 (MQ=255)
                                   tgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGGGACTCAGGTTGGCTTTGGTAACAACGCGCCCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGAGCTAGCCCTG‑TTTTTTTTCGGGAGATGAATATGACTa                                                                                                                   >  2:1684255/1‑150 (MQ=255)
                                    gcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATAc                                                                                                                  >  1:1649863/1‑150 (MQ=255)
                                     cAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACg                                                                                                                 >  1:1403000/1‑150 (MQ=255)
                                      aGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGt                                                                                                                >  1:409999/1‑150 (MQ=255)
                                      aGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGt                                                                                                                <  1:1451279/150‑1 (MQ=255)
                                       gTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGtt                                                                                                               >  2:1785677/1‑150 (MQ=255)
                                           aCCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGttatta                                                                                                           >  2:1078441/1‑150 (MQ=255)
                                             cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGACCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACt                                                                                                         >  2:1334893/1‑150 (MQ=255)
                                              aGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTc                                                                                                        <  1:1222408/150‑1 (MQ=255)
                                               gACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTATTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTcc                                                                                                       <  1:717490/150‑1 (MQ=255)
                                                    cATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCg                                                                                                  >  2:1359183/1‑150 (MQ=255)
                                                    cATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCg                                                                                                  <  2:321648/150‑1 (MQ=255)
                                                     aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                                 <  1:1230007/150‑1 (MQ=255)
                                                     aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                                 <  2:1550737/150‑1 (MQ=255)
                                                      tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                <  2:775805/150‑1 (MQ=255)
                                                      tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                <  2:1477277/150‑1 (MQ=255)
                                                      tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                <  2:1763041/150‑1 (MQ=255)
                                                      tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                >  2:520559/1‑150 (MQ=255)
                                                        tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                              >  1:1584981/1‑150 (MQ=255)
                                                        tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                              >  1:1468621/1‑150 (MQ=255)
                                                             cctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTcc                                                                                         <  1:1160780/150‑1 (MQ=255)
                                                               tccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAg                                                                                       >  2:649871/1‑150 (MQ=255)
                                                               cccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTAATAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAg                                                                                       <  2:562423/149‑1 (MQ=255)
                                                                ccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGAGCAAGCGCTG‑TTTTTTTTCTGGACAAGAATATGAATACGTAATTACTCCTTGCGGCACTTTCCAGt                                                                                      >  2:339189/1‑150 (MQ=255)
                                                                  gTCCACGTGACTCAGGGTGGCTTTGGTAACAACGCGACCGCTCATCATGACTAATACATCATTTTTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCa                                                                                    <  2:142059/150‑1 (MQ=255)
                                                                   tccaCGGGACTCAGGTTGGCTTTGGTAACAACGCTACCGCTCATCAGTAATAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAg                                                                                   <  2:1113439/147‑1 (MQ=255)
                                                                   tCAAAGTGACTCAGGTTGGATTTGGTAACAACGCGTCCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAg                                                                                   <  2:938391/150‑1 (MQ=255)
                                                                          gACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCt                                                                            <  2:511758/150‑1 (MQ=255)
                                                                           aCTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCtt                                                                           >  1:1488949/1‑150 (MQ=255)
                                                                               aGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAAt                                                                       >  2:254975/1‑150 (MQ=255)
                                                                                ggTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGGTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATa                                                                      >  2:1168044/1‑150 (MQ=255)
                                                                                 gTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATAc                                                                     <  1:555601/150‑1 (MQ=255)
                                                                                      cTTTGGTAACAACGCGACCGCTCATCAGTAGTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAccc                                                                >  1:720948/1‑150 (MQ=255)
                                                                                      cTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAccc                                                                >  1:1030762/1‑150 (MQ=255)
                                                                                      cTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAccc                                                                >  1:1708468/1‑150 (MQ=255)
                                                                                      cTTTGGTAACAACGCGAACGCTCATCAGTACTAATACATCATTAGGATTACAGAAACAGGGCGCAAGCCCTG‑GTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGAGTCACTTTCCCTTCAGATAACCTTTAATAAGAccc                                                                >  1:290951/1‑150 (MQ=255)
                                                                                        ttGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAg                                                              <  2:1754682/150‑1 (MQ=255)
                                                                                           gTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCa                                                            >  1:1356115/1‑150 (MQ=255)
                                                                                               caacGCTACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACCGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAcggg                                                       >  2:695787/1‑150 (MQ=255)
                                                                                                  cgcgACCGCTCATCATTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtg                                                    <  2:1013127/150‑1 (MQ=255)
                                                                                                  cgcgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtg                                                    <  2:1020594/150‑1 (MQ=255)
                                                                                                    cgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                                  <  1:1785677/150‑1 (MQ=255)
                                                                                                    cgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                                  <  1:1073670/150‑1 (MQ=255)
                                                                                                       ccGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACGTTCCAGTCAGATAACCTTTACTACGACCCAGCAAGGGGATGTGtat                                               >  2:278347/1‑150 (MQ=255)
                                                                                                         gCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATAc                                             >  2:477282/1‑150 (MQ=255)
                                                                                                          ggaatcaCTACGAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATAcc                                            <  2:1044883/147‑1 (MQ=255)
                                                                                                            catcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCATAAAACATTTAATACGACCCAGCAAGGGGATGTGTATACCat                                          >  1:107939/1‑150 (MQ=255)
                                                                                                                     tAATACATCATTTGTATTACAGAAACAGTGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGa                                 <  1:1143400/150‑1 (MQ=255)
                                                                                                                       aTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAg                               <  1:1298061/150‑1 (MQ=255)
                                                                                                                          catcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGCCCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCa                            >  2:1286438/1‑150 (MQ=255)
                                                                                                                          catcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCa                            >  1:486427/1‑150 (MQ=255)
                                                                                                                          catcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCa                            >  2:155321/1‑150 (MQ=255)
                                                                                                                               tttGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCtt                       <  2:967106/150‑1 (MQ=255)
                                                                                                                               tttGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCtt                       >  1:118879/1‑150 (MQ=255)
                                                                                                                                 tGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGAACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTAc                     >  1:292682/1‑150 (MQ=255)
                                                                                                                                     gaaaaaaaaCAGGGCTTGCGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTTTTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCaa                 >  1:746620/7‑150 (MQ=255)
                                                                                                                                       aCAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGTATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc               <  2:720948/150‑1 (MQ=255)
                                                                                                                                       aCAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc               >  1:176989/1‑150 (MQ=255)
                                                                                                                                       aCAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTGGCTGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGCTGTGGATACCATTATTCCTGAAGTCACTCTTACTCAAta               >  2:1116184/1‑149 (MQ=255)
                                                                                                                                        cAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCt              >  1:1623233/1‑150 (MQ=255)
                                                                                                                                         aaaaaCAGGGCTTGCGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGAGTATACCATTATTCCTGAAGTCACTCTTACTCAATCtt             >  1:9932/3‑150 (MQ=255)
                                                                                                                                         aGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCtt             >  2:954590/1‑150 (MQ=255)
                                                                                                                                          gAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTg            <  2:1433315/150‑1 (MQ=255)
                                                                                                                                            aaCAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc          <  1:1359183/150‑1 (MQ=255)
                                                                                                                                            aaCAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc          >  2:123260/1‑150 (MQ=255)
                                                                                                                                            aaCAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc          >  1:800111/1‑150 (MQ=255)
                                                                                                                                             aCAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCt         <  2:874788/150‑1 (MQ=255)
                                                                                                                                                gggCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg      <  2:546368/150‑1 (MQ=255)
                                                                                                                                                 ggCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc     >  1:928571/1‑150 (MQ=255)
                                                                                                                                                  gcgcAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCCGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCa    <  2:1403000/150‑1 (MQ=255)
                                                                                                                                                   cgcAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCag   >  1:636377/1‑150 (MQ=255)
                                                                                                                                                    gcAAGCCCTG‑TTTTTTTTCGCGCGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCACGTGGTTGTGTATACCATTATTCCTGAAGTCACTCTTACTCAAACTTGTCTGTGCaga  <  2:90004/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGA  >  NC_000913/1104783‑1105080

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