Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,104,934 A→C 5.2% intergenic (+32/‑27) 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,9340AC 5.2% 422.7 / 8.4 139intergenic (+32/‑27)csgA/csgCcurlin subunit, amyloid curli fibers, cryptic/curli assembly protein
Reads supporting (aligned to +/- strand):  ref base A (59/71);  new base C (3/5);  total (63/76)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 139 total aligned reads displayed.

CGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAG  >  NC_000913/1104787‑1105081
                                                                                                                                                   |                                                                                                                                                    
  gTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGccc                                                                                                                                                  <  1:1823969/150‑1 (MQ=255)
      aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGtt                                                                                                                                             >  2:169837/1‑150 (MQ=255)
      aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGtt                                                                                                                                             >  2:56974/1‑150 (MQ=255)
      aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGACCTGtt                                                                                                                                             >  1:658219/1‑150 (MQ=255)
        cAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACTGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTGtttc                                                                                                                                           <  1:391189/150‑2 (MQ=255)
        cAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTGtttc                                                                                                                                           <  2:975380/150‑2 (MQ=255)
          gTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTGtttctg                                                                                                                                         <  2:1342194/150‑4 (MQ=255)
           ttCGGTGGTTGCAAAGTTGCTTCAGTTTAACAGACTTCATATAAATACTCCTTAAACGTGACTCATTTTGGCTTTGGTAACAAATCGACATATCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTCAAGCCCTGttttttt                                                                                                                                        <  2:1288252/150‑1 (MQ=255)
             gggtgggGGCAACGGTGCTGCATTTTACCAGACTGAATCTAACTCCTCAGTCAAAGTGACGCAGGTTGGCGTTGGTAACAACGCGACCGCTCATCAGTACTAATACAGCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTc                                                                                                                                      <  2:1003612/149‑1 (MQ=255)
             cggcggTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTc                                                                                                                                      <  2:1175579/150‑1 (MQ=255)
              gggggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCg                                                                                                                                     <  1:618987/147‑1 (MQ=255)
               gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGGTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAACCAGGGCGCAAGCCCTGTTTTTTTTCgg                                                                                                                                    <  2:249228/150‑1 (MQ=255)
               gtggtgGCAAAGGTGCTGCAGGTGTCCAGCCTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTTGTAACAACGCGACCGCGCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCgg                                                                                                                                    <  2:49675/150‑1 (MQ=255)
                tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACTACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCggg                                                                                                                                   <  2:1548391/150‑1 (MQ=255)
                tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCggg                                                                                                                                   <  2:1488725/150‑1 (MQ=255)
                tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCggg                                                                                                                                   <  2:1189460/150‑1 (MQ=255)
                tggtggCAACGGTGCGGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCggg                                                                                                                                   <  1:1591924/150‑1 (MQ=255)
                  gtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGgag                                                                                                                                 >  1:874788/1‑150 (MQ=255)
                  gtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGgag                                                                                                                                 <  2:405638/150‑1 (MQ=255)
                      cAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGagaaga                                                                                                                             <  2:1220223/150‑1 (MQ=255)
                              ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                     >  2:892825/1‑150 (MQ=255)
                              ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                     >  1:1395053/1‑150 (MQ=255)
                              ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                     >  2:709457/1‑150 (MQ=255)
                               tgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATa                                                                                                                    >  1:1725827/1‑150 (MQ=255)
                               tgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGGGACTCAGGTTGGCTTTGGTAACAACGCGCCCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGAGCTAGCCCTGTTTTTTTTCGGGAGATGAATATGACTa                                                                                                                    >  2:1684255/1‑150 (MQ=255)
                                gcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATAc                                                                                                                   >  1:1649863/1‑150 (MQ=255)
                                 cAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACg                                                                                                                  >  1:1403000/1‑150 (MQ=255)
                                  aGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGt                                                                                                                 >  1:409999/1‑150 (MQ=255)
                                  aGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGt                                                                                                                 <  1:1451279/150‑1 (MQ=255)
                                   gTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGtt                                                                                                                >  2:1785677/1‑150 (MQ=255)
                                       aCCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGttatta                                                                                                            >  2:1078441/1‑150 (MQ=255)
                                         cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGACCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACt                                                                                                          >  2:1334893/1‑150 (MQ=255)
                                          aGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTc                                                                                                         <  1:1222408/150‑1 (MQ=255)
                                                cATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCg                                                                                                   >  2:1359183/1‑150 (MQ=255)
                                                cATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCg                                                                                                   <  2:321648/150‑1 (MQ=255)
                                                 aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                                  <  1:1230007/150‑1 (MQ=255)
                                                 aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                                  <  2:1550737/150‑1 (MQ=255)
                                                  tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                 >  2:520559/1‑150 (MQ=255)
                                                  tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                 <  2:1763041/150‑1 (MQ=255)
                                                  tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                 <  2:1477277/150‑1 (MQ=255)
                                                  tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                 <  2:775805/150‑1 (MQ=255)
                                                    tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                               >  1:1468621/1‑150 (MQ=255)
                                                    tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                               >  1:1584981/1‑150 (MQ=255)
                                                         cctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTcc                                                                                          <  1:1160780/150‑1 (MQ=255)
                                                           tccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAg                                                                                        >  2:649871/1‑150 (MQ=255)
                                                           cccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTAATAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAg                                                                                        <  2:562423/149‑1 (MQ=255)
                                                            ccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGAGCAAGCGCTGTTTTTTTTCTGGACAAGAATATGAATACGTAATTACTCCTTGCGGCACTTTCCAGt                                                                                       >  2:339189/1‑150 (MQ=255)
                                                             cGTCAACGTGACTCAGGTTGTCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTc                                                                                      <  1:542844/150‑1 (MQ=255)
                                                              gTCCACGTGACTCAGGGTGGCTTTGGTAACAACGCGACCGCTCATCATGACTAATACATCATTTTTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCa                                                                                     <  2:142059/150‑1 (MQ=255)
                                                               tccaCGGGACTCAGGTTGGCTTTGGTAACAACGCTACCGCTCATCAGTAATAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAg                                                                                    <  2:1113439/147‑1 (MQ=255)
                                                               tCAAAGTGACTCAGGTTGGATTTGGTAACAACGCGTCCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAg                                                                                    <  2:938391/150‑1 (MQ=255)
                                                                      gACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCt                                                                             <  2:511758/150‑1 (MQ=255)
                                                                       aCTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCtt                                                                            >  1:1488949/1‑150 (MQ=255)
                                                                           aGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAAt                                                                        >  2:254975/1‑150 (MQ=255)
                                                                            ggTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGGTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATa                                                                       >  2:1168044/1‑150 (MQ=255)
                                                                                 gCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                                  <  1:310719/150‑1 (MQ=255)
                                                                                  cTTTGGTAACAACGCGACCGCTCATCAGTAGTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAccc                                                                 >  1:720948/1‑150 (MQ=255)
                                                                                  cTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAccc                                                                 >  1:1030762/1‑150 (MQ=255)
                                                                                  cTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAccc                                                                 >  1:1708468/1‑150 (MQ=255)
                                                                                    ttGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAg                                                               <  2:1754682/150‑1 (MQ=255)
                                                                                       gTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCa                                                             >  1:1356115/1‑150 (MQ=255)
                                                                                       gTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCaa                                                            >  1:967106/1‑150 (MQ=255)
                                                                                         aacaacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgg                                                          <  1:709457/150‑1 (MQ=255)
                                                                                           caacGCTACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACCGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAcggg                                                        >  2:695787/1‑150 (MQ=255)
                                                                                              cgcgACCGCTCATCATTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtg                                                     <  2:1013127/150‑1 (MQ=255)
                                                                                              cgcgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtg                                                     <  2:1020594/150‑1 (MQ=255)
                                                                                                cgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                                   <  1:1785677/150‑1 (MQ=255)
                                                                                                cgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                                   <  1:1073670/150‑1 (MQ=255)
                                                                                                   ccGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACGTTCCAGTCAGATAACCTTTACTACGACCCAGCAAGGGGATGTGtat                                                >  2:278347/1‑150 (MQ=255)
                                                                                                     gCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATAc                                              >  2:477282/1‑150 (MQ=255)
                                                                                                      ggaatcaCTACGAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATAcc                                             <  2:1044883/147‑1 (MQ=255)
                                                                                                        catcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCATAAAACATTTAATACGACCCAGCAAGGGGATGTGTATACCat                                           >  1:107939/1‑150 (MQ=255)
                                                                                                            aGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCattatt                                       >  1:642926/1‑150 (MQ=255)
                                                                                                                 tAATACATCATTTGTATTACAGAAACAGTGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGa                                  <  1:1143400/150‑1 (MQ=255)
                                                                                                                   aTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAg                                <  1:1298061/150‑1 (MQ=255)
                                                                                                                      catcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGCCCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCa                             >  2:1286438/1‑150 (MQ=255)
                                                                                                                      catcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCa                             >  1:486427/1‑150 (MQ=255)
                                                                                                                      catcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCa                             >  2:155321/1‑150 (MQ=255)
                                                                                                                         catTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCActc                          <  1:229214/150‑1 (MQ=255)
                                                                                                                           tttGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCtt                        <  2:967106/150‑1 (MQ=255)
                                                                                                                           tttGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCtt                        >  1:118879/1‑150 (MQ=255)
                                                                                                                             tGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTAc                      >  1:452764/1‑150 (MQ=255)
                                                                                                                             tGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGAACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTAc                      >  1:292682/1‑150 (MQ=255)
                                                                                                                                 gaaaaaaaaCAGGGCTTGCGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTTTTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCaa                  >  1:746620/7‑150 (MQ=255)
                                                                                                                                   aCAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGTATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc                <  2:720948/150‑1 (MQ=255)
                                                                                                                                   aCAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTGGCTGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGCTGTGGATACCATTATTCCTGAAGTCACTCTTACTCAAta                >  2:1116184/1‑149 (MQ=255)
                                                                                                                                     aaaaaCAGGGCTTGCGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGAGTATACCATTATTCCTGAAGTCACTCTTACTCAATCtt              >  1:9932/3‑150 (MQ=255)
                                                                                                                                     aGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCtt              >  2:954590/1‑150 (MQ=255)
                                                                                                                                      gAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTg             <  2:1433315/150‑1 (MQ=255)
                                                                                                                                        aaCAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc           <  1:1359183/150‑1 (MQ=255)
                                                                                                                                        aaCAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc           >  2:123260/1‑150 (MQ=255)
                                                                                                                                        aaCAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc           >  1:800111/1‑150 (MQ=255)
                                                                                                                                         aCAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCt          <  2:874788/150‑1 (MQ=255)
                                                                                                                                            gggCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg       <  2:546368/150‑1 (MQ=255)
                                                                                                                                            gggCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg       <  1:892825/150‑1 (MQ=255)
                                                                                                                                             ggCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc      >  1:928571/1‑150 (MQ=255)
                                                                                                                                             ggCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc      <  1:764872/150‑1 (MQ=255)
                                                                                                                                              gcgcAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCCGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCa     <  2:1403000/150‑1 (MQ=255)
                                                                                                                                               cgcAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCag    >  1:636377/1‑150 (MQ=255)
                                                                                                                                                gcAAGCCCTGTTTTTTTTCGCGCGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCACGTGGTTGTGTATACCATTATTCCTGAAGTCACTCTTACTCAAACTTGTCTGTGCaga   <  2:90004/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAG  >  NC_000913/1104787‑1105081

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