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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,104,9310GT 6.5% 466.8 / 12.7 167intergenic (+29/‑30)csgA/csgCcurlin subunit, amyloid curli fibers, cryptic/curli assembly protein
Reads supporting (aligned to +/- strand):  ref base G (70/85);  new base T (4/8);  total (74/93)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 170 total aligned reads displayed.

ATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGA  >  NC_000913/1104783‑1105080
                                                                                                                                                    |                                                                                                                                                     
 tGACTGTTAAACAGTGCGGTGGTGGCAACTGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCa                                                                                                                                                     <  2:733823/150‑1 (MQ=255)
 tGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCa                                                                                                                                                     <  2:2057162/150‑1 (MQ=255)
   aCGGTTAAACAGTTCGTTGGTGGAAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGATTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAg                                                                                                                                                   <  1:1970319/150‑1 (MQ=255)
   aCGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAg                                                                                                                                                   >  1:1926733/1‑150 (MQ=255)
    cGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGc                                                                                                                                                  >  1:1241624/1‑150 (MQ=255)
    cGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGc                                                                                                                                                  >  2:545029/1‑150 (MQ=255)
    cGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGc                                                                                                                                                  >  1:1707809/1‑150 (MQ=255)
          aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTTCAGAAACAGGGCGCAAGCCCTGtt                                                                                                                                            >  1:39919/1‑150 (MQ=255)
            cAGTTCGGTGGTGGCAACGGTGCGGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGGTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTGtttc                                                                                                                                          <  1:870233/150‑2 (MQ=255)
            cAGTGCGGTTGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCTTCAACGTGACTCAGGTTGGCTTTTGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTGtttc                                                                                                                                          <  1:1901726/150‑2 (MQ=255)
             aGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCTACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGttttt                                                                                                                                         >  1:80473/1‑150 (MQ=255)
             aGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTGtttct                                                                                                                                         <  1:1596410/150‑3 (MQ=255)
             aGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTGtttct                                                                                                                                         <  2:2019926/150‑3 (MQ=255)
                gcGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGtttttttt                                                                                                                                      <  1:597841/149‑1 (MQ=255)
                   gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCgg                                                                                                                                   <  2:1480075/150‑1 (MQ=255)
                    ttgtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCTTCAACGTTACTCATGTTTGCTTTGGTAACAACTCTACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCggg                                                                                                                                  <  2:270215/148‑1 (MQ=255)
                    tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCggg                                                                                                                                  <  2:1212723/150‑1 (MQ=255)
                    tggtggCAACGGTGATGAATTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGAATCAGGTTGGCTTTGGTAACAACGCGACCGCTCTTCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCggg                                                                                                                                  <  2:787350/150‑1 (MQ=255)
                    ggggggCAACGGTGCTGCAGTTGACCAGACGGCATCTAACTCCTCCGTCAACGTGACTCAGGGTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCggg                                                                                                                                  <  2:1522377/146‑1 (MQ=255)
                     ggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGga                                                                                                                                 <  2:403827/150‑1 (MQ=255)
                     ggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGga                                                                                                                                 <  2:1751050/150‑1 (MQ=255)
                        ggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCCTCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCGCTGTGTTTTTTCGGGagaa                                                                                                                              >  2:1035534/1‑150 (MQ=255)
                           aaCGGTGCTGCAGGTGACCAGACTGCATCTAACTCCTCCGGCAACGGGACTCCGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAgaagaa                                                                                                                           <  2:207869/150‑1 (MQ=255)
                            aCGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAat                                                                                                                          <  1:753900/150‑1 (MQ=255)
                                  ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                    >  1:1103549/1‑150 (MQ=255)
                                    gcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGAAACAACGCGACCGCTCATCAGTACTAATACATCATGTGTATTACAGAAACAGGGCTCAAGCCCTGTTTGTATACTCGAGAACAATATGACTaa                                                                                                                  >  2:959319/1‑149 (MQ=255)
                                     cAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTAGGGAGAAGAATATGAATACg                                                                                                                 >  1:1253061/1‑150 (MQ=255)
                                       gTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGtt                                                                                                               >  2:887855/1‑150 (MQ=255)
                                       gTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGtt                                                                                                               <  1:1802430/150‑1 (MQ=255)
                                          gACCAGACTGCATCTACCTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTACCAACGCGCCCGCTCATCAGTACTAATACATCATTTGTATTTCAGAAACAGGGCGCAAGCCCTGTGTTTTTTCGGGAGAAGAATATGAATACGttatt                                                                                                            >  2:2107045/1‑150 (MQ=255)
                                          gACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGttatt                                                                                                            <  1:1318136/150‑1 (MQ=255)
                                          gACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGttatt                                                                                                            >  2:628484/1‑150 (MQ=255)
                                                aCTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCt                                                                                                      >  1:1583832/1‑150 (MQ=255)
                                                   gCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGc                                                                                                   <  2:1099672/150‑1 (MQ=255)
                                                   gCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGc                                                                                                   >  1:2020111/1‑150 (MQ=255)
                                                     aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                                 <  2:1726829/150‑1 (MQ=255)
                                                      tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                <  2:1289460/150‑1 (MQ=255)
                                                       cTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCa                                                                                               >  2:540143/1‑150 (MQ=255)
                                                        tgaCTCCTCCGTCAACGTGACTCAGGTTGGCTATGGTAACAACGCGACCGCTCATCAGTACTAATACATAATTTGTATTACAGAAACAGGGCTCAAGCCCTGTTATTTTTAGGGATCAGAATATGAATACGTTATTACTCATTGCGGCAc                                                                                              >  1:129128/3‑150 (MQ=255)
                                                        tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                              <  1:79566/150‑1 (MQ=255)
                                                        tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                              >  1:881316/1‑150 (MQ=255)
                                                        tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                              <  1:1525469/150‑1 (MQ=255)
                                                         aaCTCCTCCGTCAACGTGACTCAGGGTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTTTTACAGAAACAGGGCGCAAGCCCTGGTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCaat                                                                                             >  1:1776678/1‑148 (MQ=255)
                                                           ctcctcCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAATCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACttt                                                                                           >  2:1609339/1‑150 (MQ=255)
                                                            acctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTc                                                                                          <  2:125741/149‑1 (MQ=255)
                                                               cccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATCCATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAg                                                                                       <  1:1731189/149‑1 (MQ=255)
                                                                 cGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCAAGTc                                                                                     >  2:2010135/1‑150 (MQ=255)
                                                                   tCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAg                                                                                   <  1:2202136/150‑1 (MQ=255)
                                                                   tCAACGTGACTCAGGTAGGCTTTGGTAACAACGCGACCGCTCATCAGTTATAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAg                                                                                   <  2:459676/150‑1 (MQ=255)
                                                                          gACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCt                                                                            >  2:671828/1‑150 (MQ=255)
                                                                            cTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCttt                                                                          >  1:2112078/1‑150 (MQ=255)
                                                                            cTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCttt                                                                          <  1:1064438/150‑1 (MQ=255)
                                                                            atctggTTGGCTTTGGTAAGAACGCGACCTCTCATCAGTACTTATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTGGGGAGAAGAATATGAATACGTTATTACTCCTTGCTGCACTTTCCAGTCAGATAACCttt                                                                          <  1:1496928/146‑1 (MQ=255)
                                                                             gcAGGTGGGCTTTGGTTACATCGCGACCGCTCTTCAGTACAAATACAACATTTGTATTACAGTAACAGGGCGCAAGCCCTGTTTTTTTTCGTGAGAAGAATATGAATTCGTTATTACTCCTTGCGGCACTTTCGAGTCAGATAACCTTTa                                                                         <  1:212348/149‑1 (MQ=255)
                                                                              cAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTaa                                                                        >  1:2141197/1‑150 (MQ=255)
                                                                              cAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTaa                                                                        <  2:198347/150‑1 (MQ=255)
                                                                              cAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCATGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTaa                                                                        >  2:1409373/1‑150 (MQ=255)
                                                                                ggTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATa                                                                      <  2:939036/150‑1 (MQ=255)
                                                                                  ttGGCTTTGGGAACCACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACACGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACg                                                                    <  2:221633/150‑1 (MQ=255)
                                                                                       tttGGTAACAACGCGACCGCTCATCAGTACTAATACCTCAGTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGTTAGAATAATATGTCTACGTTATTACTCCGTGCGGCACTTTCCAGTCAGATAACCTTTAAGACGAcccc                                                               >  2:1581377/1‑149 (MQ=255)
                                                                                        ttGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAg                                                              <  2:2032725/150‑1 (MQ=255)
                                                                                        ttGGGAACAACTCGACCGGTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCGCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAg                                                              <  2:167584/150‑1 (MQ=255)
                                                                                          ggTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCa                                                            <  2:1944248/150‑1 (MQ=255)
                                                                                          ggTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCa                                                            >  1:377965/1‑150 (MQ=255)
                                                                                              acaacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAggg                                                        <  1:978374/150‑1 (MQ=255)
                                                                                              acaacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAggg                                                        <  2:302322/150‑1 (MQ=255)
                                                                                                aacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                      <  2:1483369/150‑1 (MQ=255)
                                                                                                aacGCGACCGCTCATCAGTACTAATAAATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                      >  1:1399400/1‑150 (MQ=255)
                                                                                                 acGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAt                                                     <  2:1303157/150‑1 (MQ=255)
                                                                                                  cgcgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtg                                                    <  2:968958/150‑1 (MQ=255)
                                                                                                    cgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                                  <  2:569634/150‑1 (MQ=255)
                                                                                                    cgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                                  <  2:1787311/150‑1 (MQ=255)
                                                                                                      cccGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAATAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGta                                                <  1:926607/149‑1 (MQ=255)
                                                                                                      aCCGGTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGACACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGta                                                <  1:1613728/150‑1 (MQ=255)
                                                                                                       ccGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGtat                                               >  2:1774670/1‑150 (MQ=255)
                                                                                                       ccGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGtat                                               >  2:1765784/1‑150 (MQ=255)
                                                                                                        cGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGtata                                              <  2:1926733/150‑1 (MQ=255)
                                                                                                         gCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATAc                                             <  1:1988129/150‑1 (MQ=255)
                                                                                                                   aCTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCt                                   >  1:2010340/1‑150 (MQ=255)
                                                                                                                    cTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTg                                  >  1:1874602/1‑150 (MQ=255)
                                                                                                                      aaTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGaa                                >  1:948170/1‑150 (MQ=255)
                                                                                                                      aaTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACCTTATTACTCCATGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTAAACCATTATTCATGaa                                >  2:1188588/1‑150 (MQ=255)
                                                                                                                       aTACATGATTTGTATTACAGAAACAGGTCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAg                               <  1:706539/150‑1 (MQ=255)
                                                                                                                           atcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCAc                           <  1:1817484/150‑1 (MQ=255)
                                                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATATCCTTTAATACGACCCATCAAGGGGATGTGTATAACATTATTCCTGAAGTCAct                          >  2:94092/1‑150 (MQ=255)
                                                                                                                              atTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCActct                        <  2:514791/150‑1 (MQ=255)
                                                                                                                                   tATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACta                   >  2:1806314/1‑149 (MQ=255)
                                                                                                                                    cgaaaaaaaaCAGGGCTTGCGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTGTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCa                  >  2:933326/8‑150 (MQ=255)
                                                                                                                                    aTTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGTATGTGTATACCATTATTCCTGAAGTCACTCTTAATCa                  >  1:1551199/1‑150 (MQ=255)
                                                                                                                                    aTTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCa                  >  1:628017/1‑150 (MQ=255)
                                                                                                                                       aCAGAAACAGGGCTAAAGCCATGTTTTTTTTCGGGAGAAGAATCTGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCAAGCAAGGGGTTGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc               <  1:959319/150‑1 (MQ=255)
                                                                                                                                        cAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATAAGAATACGTTATTACTCCTTGCGGCACTTTCAAGTCAGATAACCTTTAATAATACCAAGCAAGGGGATGTGTATACCAGTATTACTGAAGTCACTCTTACTCAATCt              >  2:1060908/1‑150 (MQ=255)
                                                                                                                                          gAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCACTCTTg            >  1:987247/1‑150 (MQ=255)
                                                                                                                                            aaCATGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGGCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACAATTATTCCTGAAGTCACTCTTACTCAATCTTGTc          <  2:1885938/150‑1 (MQ=255)
                                                                                                                                            aaCACGGCGCAAGCCCTTTTTTTTTTCGGGAGAAGAATATGAATACGTTAGTACTCCTTGCGGCACTTTCCAGTCAGATACCCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc          <  2:350149/150‑1 (MQ=255)
                                                                                                                                             aCAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCt         <  2:881316/150‑1 (MQ=255)
                                                                                                                                              cAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATGCGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtg        <  2:1646460/150‑1 (MQ=255)
                                                                                                                                              cAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtg        >  1:484255/1‑150 (MQ=255)
                                                                                                                                                gggCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg      >  1:84908/1‑150 (MQ=255)
                                                                                                                                                    gcAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTTTCTGTGCaga  >  1:194295/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGA  >  NC_000913/1104783‑1105080

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