Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,104,933 A→G 5.2% intergenic (+31/‑28) 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,9330AG 5.2% 1799.2 / 52.9 533intergenic (+31/‑28)csgA/csgCcurlin subunit, amyloid curli fibers, cryptic/curli assembly protein
Reads supporting (aligned to +/- strand):  ref base A (246/258);  new base G (22/6);  total (268/265)
Fisher's exact test for biased strand distribution p-value = 2.93e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 540 total aligned reads displayed.

TGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAG  >  NC_000913/1104784‑1105081
                                                                                                                                                     |                                                                                                                                                      
    ggTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGcc                                                                                                                                                    <  2:4212215/150‑1 (MQ=255)
    ggTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTCGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGcc                                                                                                                                                    <  1:6844566/150‑1 (MQ=255)
     gTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGccc                                                                                                                                                   >  1:4296131/1‑150 (MQ=255)
      ttAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑t                                                                                                                                                 >  1:6679424/1‑150 (MQ=255)
      ttAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑t                                                                                                                                                 >  1:8431170/1‑150 (MQ=255)
         aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑tt                                                                                                                                             >  2:3470497/1‑150 (MQ=255)
          aCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑ttt                                                                                                                                            >  2:5178038/1‑150 (MQ=255)
          aCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑ttt                                                                                                                                            >  2:2134224/1‑150 (MQ=255)
           cAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑tttt                                                                                                                                           <  2:7046978/150‑1 (MQ=255)
            aGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACACGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCC‑TG‑tttct                                                                                                                                          <  2:7408167/150‑3 (MQ=255)
              ttCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑ttttttt                                                                                                                                        >  2:1320118/1‑150 (MQ=255)
               tCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑tttttttt                                                                                                                                       <  1:3793631/150‑1 (MQ=255)
                cGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTTGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTc                                                                                                                                      <  2:6571301/150‑1 (MQ=255)
                cGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTc                                                                                                                                      >  2:1020855/1‑150 (MQ=255)
                cGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTc                                                                                                                                      <  1:2368112/150‑1 (MQ=255)
                 ggtggtGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCg                                                                                                                                     <  2:5790123/150‑1 (MQ=255)
                  gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCgg                                                                                                                                    <  2:8417514/150‑1 (MQ=255)
                  gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCgg                                                                                                                                    <  2:2299119/150‑1 (MQ=255)
                  gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCgg                                                                                                                                    <  1:583530/150‑1 (MQ=255)
                    ggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGga                                                                                                                                  <  1:562656/150‑1 (MQ=255)
                      tggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGgaga                                                                                                                                >  1:2691526/1‑150 (MQ=255)
                        gCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGagaag                                                                                                                              >  1:7195193/1‑150 (MQ=255)
                             ggTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAatat                                                                                                                         <  1:6371381/150‑1 (MQ=255)
                             ggTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAatat                                                                                                                         >  2:3266425/1‑150 (MQ=255)
                                gctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGaa                                                                                                                      <  1:465364/150‑1 (MQ=255)
                                 ctgcAGGTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                     <  2:2820622/150‑1 (MQ=255)
                                       ttGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACCGGGCGCAAGCCC‑TG‑ATTTTTTTCGGGAGAAGAATATGAATACGtta                                                                                                               <  1:1635886/150‑1 (MQ=255)
                                       ttGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGtta                                                                                                               >  2:7133402/1‑150 (MQ=255)
                                        tGACCAGACTGCATCTAACTCCTCTGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGttat                                                                                                              <  1:3592104/150‑1 (MQ=255)
                                          aCCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGttatta                                                                                                            >  1:7653855/1‑150 (MQ=255)
                                          aCCAGACTGCATCTAACGCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAAGACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGATTATGAATACGttatta                                                                                                            <  2:4794865/150‑1 (MQ=255)
                                            cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACt                                                                                                          <  1:8108816/150‑1 (MQ=255)
                                            cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACt                                                                                                          >  2:7315074/1‑150 (MQ=255)
                                                 tGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTg                                                                                                     >  1:7139549/1‑150 (MQ=255)
                                                  gCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGc                                                                                                    <  2:6688443/150‑1 (MQ=255)
                                                    aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                                  <  2:7923769/150‑1 (MQ=255)
                                                      cTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCa                                                                                                >  1:7079830/1‑150 (MQ=255)
                                                      cTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCa                                                                                                >  1:6132869/1‑150 (MQ=255)
                                                       tAACTCCTCCGTCAACTTTAATCATGTTTGCTTTGTTAAAAACGCGTACTCTAATCAGTACTTATACATCATTTTTATTACATAAACAGGGCGCAAGCCCTTT‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCa                                                                                                <  1:4732272/150‑1 (MQ=255)
                                                       tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                               <  2:4031341/150‑1 (MQ=255)
                                                       tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                               >  1:4985039/1‑150 (MQ=255)
                                                        aaCTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACt                                                                                              >  2:3013639/1‑150 (MQ=255)
                                                        aaCTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCACTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACt                                                                                              <  2:2599473/150‑1 (MQ=255)
                                                         aCTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACtt                                                                                             >  2:7614812/1‑150 (MQ=255)
                                                             ctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCa                                                                                         >  1:7949145/1‑150 (MQ=255)
                                                             ctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCa                                                                                         <  1:2139057/150‑1 (MQ=255)
                                                                  tCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAg                                                                                    >  1:2544200/1‑150 (MQ=255)
                                                                       gTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAAc                                                                               <  2:5852239/150‑1 (MQ=255)
                                                                           cTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCttt                                                                           <  1:548433/150‑1 (MQ=255)
                                                                             cAGGTTGGCTTTGGTAACAACGCGACCGCTCATCATTACTAATACATAATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTaa                                                                         <  2:4489419/150‑1 (MQ=255)
                                                                             cAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTaa                                                                         <  2:1182030/150‑1 (MQ=255)
                                                                             cAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTaa                                                                         <  1:1020855/150‑1 (MQ=255)
                                                                               ggTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATa                                                                       >  2:8234249/1‑150 (MQ=255)
                                                                                gTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATAc                                                                      >  1:3560446/1‑150 (MQ=255)
                                                                                gTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCATTGCGGCACTTTCCAGTCAGATAACCTTTAATaa                                                                      >  2:1647992/1‑149 (MQ=255)
                                                                                  tGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGa                                                                    <  2:1248861/150‑1 (MQ=255)
                                                                                   ggCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAc                                                                   <  2:2792430/150‑1 (MQ=255)
                                                                                       ttGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAg                                                               <  2:7810158/150‑1 (MQ=255)
                                                                                        tGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGc                                                              >  1:1291383/1‑150 (MQ=255)
                                                                                         ggTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCa                                                             >  2:5096003/1‑150 (MQ=255)
                                                                                         ggTAACAACGCGACCGCTCATAAGTACTAATACATCATTTGGATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCa                                                             >  2:4942502/1‑150 (MQ=255)
                                                                                          gTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCaa                                                            <  2:621279/150‑1 (MQ=255)
                                                                                           tAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTc                                                                                      <  2:1836542/123‑1 (MQ=255)
                                                                                           tAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTc                                                                                      >  1:1836542/1‑123 (MQ=255)
                                                                                             acaacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAggg                                                         >  2:4860679/1‑150 (MQ=255)
                                                                                                acGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAt                                                      <  1:5415912/150‑1 (MQ=255)
                                                                                                    gACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATgtgt                                                  <  1:1163761/150‑1 (MQ=255)
                                                                                                    gACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATgtgt                                                  <  1:4869412/150‑1 (MQ=255)
                                                                                                      ccGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGtat                                                <  1:2231064/150‑1 (MQ=255)
                                                                                                       cGCTCATCAGTACTAATACATCATTTGTATCACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGtata                                               >  2:1975021/1‑150 (MQ=255)
                                                                                                       cGCTCATCAGTACTAATACATCATTTGTAGTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTACGGACAAGGATAAGAATAAGTTATTACTCCGTTCGGCACATTCCAGGACGATTACCGTTAATGCGCCCCAGCCAGGGGGTGTGtatt                                               >  2:549991/1‑149 (MQ=255)
                                                                                                          tcatcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCa                                            <  1:8282508/150‑1 (MQ=255)
                                                                                                             tcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCatta                                         >  2:3115876/1‑150 (MQ=255)
                                                                                                                gTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                      <  1:3407181/150‑1 (MQ=255)
                                                                                                                gTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                      <  1:3266425/150‑1 (MQ=255)
                                                                                                                gTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                      >  1:5070962/1‑150 (MQ=255)
                                                                                                                    tAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGa                                  >  2:1609776/1‑150 (MQ=255)
                                                                                                                    tAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGa                                  <  1:5220933/150‑1 (MQ=255)
                                                                                                                    tAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGa                                  >  2:6023878/1‑150 (MQ=255)
                                                                                                                     aaTACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGaa                                 <  2:3779877/150‑1 (MQ=255)
                                                                                                                        aCATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTc                              >  2:7470911/1‑150 (MQ=255)
                                                                                                                         catcatTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCa                             <  1:4455350/150‑1 (MQ=255)
                                                                                                                          atcatTTGTATTACAGAAACAGGGCGCAAGCCC‑TGTTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCa                             >  2:5285298/1‑150 (MQ=255)
                                                                                                                          atcatTTGTATTACAGAAACAGGGCGCAAGCCC‑AG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCAc                            >  1:7355870/1‑150 (MQ=255)
                                                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCAct                           <  1:5178039/150‑1 (MQ=255)
                                                                                                                            catTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCActc                          <  1:6584367/150‑1 (MQ=255)
                                                                                                                             atTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCActct                         >  2:1408131/1‑150 (MQ=255)
                                                                                                                             atTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCActct                         >  1:3951677/1‑150 (MQ=255)
                                                                                                                              tttGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCtt                        >  1:5957777/1‑150 (MQ=255)
                                                                                                                                   aTTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCa                   <  2:6132869/150‑1 (MQ=255)
                                                                                                                                    gaaaaaaaaCAGGGCTTGCGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGGGTATACCATTATTCCTGAAGTCACTCTTACTCaa                  >  2:4313/7‑150 (MQ=255)
                                                                                                                                    gaaaaaaaaCAGGGCTTGCGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGGTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCaa                  >  2:4280949/7‑150 (MQ=255)
                                                                                                                                      aCAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc                <  2:6999580/150‑1 (MQ=255)
                                                                                                                                       cAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCt               >  2:3694529/1‑150 (MQ=255)
                                                                                                                                        aGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCtt              >  2:8011532/1‑150 (MQ=255)
                                                                                                                                          aaaCAGGGCTTGCGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCCAGGGGCTGTGTATACCCTTATTCCTGAAGTCACTCTTACTCAAACTGGt            >  2:6884214/1‑150 (MQ=255)
                                                                                                                                            aCAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCt          <  1:3108654/150‑1 (MQ=255)
                                                                                                                                               gggCGCAAGCCC‑TGTTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgt        <  1:412741/150‑1 (MQ=255)
                                                                                                                                                ggCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc      >  1:7447819/1‑150 (MQ=255)
                                                                                                                                                    cAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCagag  <  1:8234469/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCC‑TG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAG  >  NC_000913/1104784‑1105081

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑