Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,104,932 C→T 5.9% intergenic (+30/‑29) 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,9320CT 5.9% 424.9 / 8.9 134intergenic (+30/‑29)csgA/csgCcurlin subunit, amyloid curli fibers, cryptic/curli assembly protein
Reads supporting (aligned to +/- strand):  ref base C (58/68);  new base T (3/5);  total (61/73)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 140 total aligned reads displayed.

ATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAG  >  NC_000913/1104783‑1105081
                                                                                                                                                     |                                                                                                                                                      
aTGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGgcgc                                                                                                                                                        >  1:1528429/1‑150 (MQ=255)
    cGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGc                                                                                                                                                    <  1:916084/150‑1 (MQ=255)
          aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGCCCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑tt                                                                                                                                             >  1:1386034/1‑150 (MQ=255)
          aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑tt                                                                                                                                             >  2:169837/1‑150 (MQ=255)
          aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑tt                                                                                                                                             >  2:56974/1‑150 (MQ=255)
           aCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑ttt                                                                                                                                            >  1:321648/1‑150 (MQ=255)
           aCAGTTCGGGGGGGGCAACGGGGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑ttt                                                                                                                                            <  2:1269050/150‑1 (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)
              gtgcGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑tttttt                                                                                                                                         <  2:1319727/147‑1 (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)
               ttCGGTGGCGGCAACGGGGCGGCAGGTGACCAGACTGCATCTAACTCCTCCGTAAACGTGTCTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCCGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttctgt                                                                                                                                        <  2:275247/150‑5 (MQ=255)
                 gggtgggGGCAACGGTGCTGCATTTTACCAGACTGAATCTAACTCCTCAGTCAAAGTGACGCAGGTTGGCGTTGGTAACAACGCGACCGCTCATCAGTACTAATACAGCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTc                                                                                                                                      <  2:1003612/149‑1 (MQ=255)
                 cggcggTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTc                                                                                                                                      <  2:1175579/150‑1 (MQ=255)
                  ggtggtGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCg                                                                                                                                     <  1:1496647/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)
                    tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCggg                                                                                                                                   <  2:1488725/150‑1 (MQ=255)
                      gtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGgag                                                                                                                                 <  2:405638/150‑1 (MQ=255)
                      gtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGgag                                                                                                                                 >  1:874788/1‑150 (MQ=255)
                      gtggCAACGGTGCTGCAGGGGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACCGGGCGCAAGCCCTG‑TTTTTTTTCGGgag                                                                                                                                 <  1:122047/150‑1 (MQ=255)
                          cAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGagaaga                                                                                                                             <  2:1220223/150‑1 (MQ=255)
                               gTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATg                                                                                                                        <  2:1016514/150‑1 (MQ=255)
                                  ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                     >  2:709457/1‑150 (MQ=255)
                                  ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                     >  2:892825/1‑150 (MQ=255)
                                  ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                     >  1:1395053/1‑150 (MQ=255)
                                   tgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATa                                                                                                                    >  1:1725827/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:1451279/150‑1 (MQ=255)
                                       gTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGtt                                                                                                                >  2:1785677/1‑150 (MQ=255)
                                           aCCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGttatta                                                                                                            >  2:1183086/1‑150 (MQ=255)
                                             cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACt                                                                                                          <  2:1725827/150‑1 (MQ=255)
                                              aGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTc                                                                                                         <  1:1222408/150‑1 (MQ=255)
                                               gACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTATTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTcc                                                                                                        <  1:717490/150‑1 (MQ=255)
                                                    cATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCg                                                                                                   <  2:321648/150‑1 (MQ=255)
                                                    cATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCg                                                                                                   >  2:1359183/1‑150 (MQ=255)
                                                     aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                                  <  2:1550737/150‑1 (MQ=255)
                                                     aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                                  <  1:1230007/150‑1 (MQ=255)
                                                      tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                 <  2:775805/150‑1 (MQ=255)
                                                      tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                 <  2:1763041/150‑1 (MQ=255)
                                                      tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                 <  2:1477277/150‑1 (MQ=255)
                                                      tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                 >  2:520559/1‑150 (MQ=255)
                                                        tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                               >  1:1468621/1‑150 (MQ=255)
                                                        tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                               >  1:1584981/1‑150 (MQ=255)
                                                             cctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTcc                                                                                          <  1:1160780/150‑1 (MQ=255)
                                                              ctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCa                                                                                         <  2:1525218/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)
                                                                 cGTCAACGTGACTCAGGTTGTCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTc                                                                                      <  1:542844/150‑1 (MQ=255)
                                                                 cGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTc                                                                                      <  2:1386034/150‑1 (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)
                                                                               aGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAAt                                                                        >  2:254975/1‑150 (MQ=255)
                                                                               aGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACACGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAAt                                                                        <  2:1156562/150‑1 (MQ=255)
                                                                                 gTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATAc                                                                      <  1:555601/150‑1 (MQ=255)
                                                                                     gCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                                  <  1:310719/150‑1 (MQ=255)
                                                                                      cTTTGGTAACAACGCGACCGCTCATCAGTAGTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAccc                                                                 >  1:720948/1‑150 (MQ=255)
                                                                                      cTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAccc                                                                 >  1:1030762/1‑150 (MQ=255)
                                                                                      cTTTGGTAACAACGCGAACGCTCATCAGTACTAATACATCATTAGGATTACAGAAACAGGGCGCAAGCCCTG‑GTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGAGTCACTTTCCCTTCAGATAACCTTTAATAAGAccc                                                                 >  1:290951/1‑150 (MQ=255)
                                                                                           gTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCa                                                             >  1:1356115/1‑150 (MQ=255)
                                                                                           gTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCaa                                                            >  1:967106/1‑150 (MQ=255)
                                                                                             aacaacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgg                                                          <  1:709457/150‑1 (MQ=255)
                                                                                               caacGCTACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACCGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAcggg                                                        >  2:695787/1‑150 (MQ=255)
                                                                                                  cgcgACCGCTCATCATTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtg                                                     <  2:1013127/150‑1 (MQ=255)
                                                                                                    cgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                                   <  1:1073670/150‑1 (MQ=255)
                                                                                                    cgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                                   <  1:1785677/150‑1 (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)
                                                                                                              tcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TATTTTTTCGGGAGAAGAATATGAATACGTTATAACTCCTTGCGGCACTTTCCAGTCAGATAACCTTAAAAACGACCCAGCAAGGGGATGTGTATACCatta                                         >  2:1052933/1‑150 (MQ=255)
                                                                                                                aGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCattatt                                       >  1:642926/1‑150 (MQ=255)
                                                                                                                     tAATACATCATTTGTATTACAGAAACAGTGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGa                                  <  1:1143400/150‑1 (MQ=255)
                                                                                                                       aTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAg                                <  1:1298061/150‑1 (MQ=255)
                                                                                                                          catcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCa                             >  2:155321/1‑150 (MQ=255)
                                                                                                                             catTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCActc                          <  1:229214/150‑1 (MQ=255)
                                                                                                                               tttGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCtt                        <  2:967106/150‑1 (MQ=255)
                                                                                                                                 tGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTAc                      >  1:452764/1‑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)
                                                                                                                                         aGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCtt              >  2:954590/1‑150 (MQ=255)
                                                                                                                                          gAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTg             <  2:1433315/150‑1 (MQ=255)
                                                                                                                                            aaCAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc           >  2:123260/1‑150 (MQ=255)
                                                                                                                                             aCAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCt          <  2:874788/150‑1 (MQ=255)
                                                                                                                                                gggCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg       <  2:546368/150‑1 (MQ=255)
                                                                                                                                                gggCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg       <  1:892825/150‑1 (MQ=255)
                                                                                                                                                 ggCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc      <  1:764872/150‑1 (MQ=255)
                                                                                                                                                    gcAAGCCCTG‑TTTTTTTTCGCGCGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCACGTGGTTGTGTATACCATTATTCCTGAAGTCACTCTTACTCAAACTTGTCTGTGCaga   <  2:90004/150‑1 (MQ=255)
                                                                                                                                                     cAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATATGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCagag  <  1:978435/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAG  >  NC_000913/1104783‑1105081

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