Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,104,933 A→G 5.3% 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.3% 517.3 / 13.3 168intergenic (+31/‑28)csgA/csgCcurlin subunit, amyloid curli fibers, cryptic/curli assembly protein
Reads supporting (aligned to +/- strand):  ref base A (74/85);  new base G (2/7);  total (76/92)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 170 total aligned reads displayed.

TGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAG  >  NC_000913/1104784‑1105081
                                                                                                                                                     |                                                                                                                                                      
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)
     gTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGccc                                                                                                                                                   <  1:940879/150‑1 (MQ=255)
         aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTTCAGAAACAGGGCGCAAGCCCTG‑‑tt                                                                                                                                             >  1:39919/1‑150 (MQ=255)
         aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑‑tt                                                                                                                                             <  2:1236040/150‑1 (MQ=255)
           cAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACTTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑tttt                                                                                                                                           <  2:1943764/150‑1 (MQ=255)
           cAGTGCGGTTGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCTTCAACGTGACTCAGGTTGGCTTTTGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑‑tttc                                                                                                                                           <  1:1901726/150‑2 (MQ=255)
            aGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑‑tttct                                                                                                                                          <  2:2019926/150‑3 (MQ=255)
            aGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑‑tttct                                                                                                                                          <  1:1596410/150‑3 (MQ=255)
                 ggtggtGGCATCGGTGCTGGAGTTGATCAGTCTGTTGCATAGTCGTCCGTCAACGTGACTCTGGTTGGCTTTGGTAACTACTCGACCGCGCATCTGTTCTATTACTTCATTTGTTTTACAGAATCAGGGCGCAAGCCCTG‑‑TTTTTTTTCg                                                                                                                                     <  1:317774/150‑1 (MQ=255)
                  gttgaggCCACGGGGCTGCAGTTGACCAGACTGGATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCGCATCAGTACTACTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCgg                                                                                                                                    <  2:170277/145‑1 (MQ=255)
                   ttgtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCTTCAACGTTACTCATGTTTGCTTTGGTAACAACTCTACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCggg                                                                                                                                   <  2:270215/148‑1 (MQ=255)
                   tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCggg                                                                                                                                   <  2:1212723/150‑1 (MQ=255)
                   tggtggCAACGGTGATGAATTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGAATCAGGTTGGCTTTGGTAACAACGCGACCGCTCTTCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCggg                                                                                                                                   <  2:787350/150‑1 (MQ=255)
                    ggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGga                                                                                                                                  <  2:1751050/150‑1 (MQ=255)
                    ggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGga                                                                                                                                  >  2:1988129/1‑150 (MQ=255)
                    ggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGga                                                                                                                                  <  2:403827/150‑1 (MQ=255)
                       ggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCCTCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCGCTG‑‑TGTTTTTTCGGGagaa                                                                                                                               >  2:1035534/1‑150 (MQ=255)
                          aaCGGTGCTGCAGGTGACCAGACTGCATCTAACTCCTCCGGCAACGGGACTCCGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAgaagaa                                                                                                                            <  2:207869/150‑1 (MQ=255)
                                 ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAAt                                                                                                                     >  1:1103549/1‑150 (MQ=255)
                                   gcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGAAACAACGCGACCGCTCATCAGTACTAATACATCATGTGTATTACAGAAACAGGGCTCAAGCCCTG‑‑TTTGTATACTCGAGAACAATATGACTaa                                                                                                                   >  2:959319/1‑149 (MQ=255)
                                    cAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTAGGGAGAAGAATATGAATACg                                                                                                                  >  1:1253061/1‑150 (MQ=255)
                                    cAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATATGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACg                                                                                                                  <  1:856139/150‑1 (MQ=255)
                                      gTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGtt                                                                                                                >  2:887855/1‑150 (MQ=255)
                                      gTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGtt                                                                                                                <  1:1802430/150‑1 (MQ=255)
                                       ttGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAAAACGtta                                                                                                               >  2:1950872/1‑150 (MQ=255)
                                         gACCAGACTGCATCTACCTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTACCAACGCGCCCGCTCATCAGTACTAATACATCATTTGTATTTCAGAAACAGGGCGCAAGCCCTG‑‑TGTTTTTTCGGGAGAAGAATATGAATACGttatt                                                                                                             >  2:2107045/1‑150 (MQ=255)
                                         gACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGttatt                                                                                                             >  2:628484/1‑150 (MQ=255)
                                         gACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGttatt                                                                                                             <  1:1318136/150‑1 (MQ=255)
                                         gACCAGACTGCATCTAACTCCTCCGTCAAAGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCATG‑‑TTTTTTTTCGTGAGAAGAATATGAATACGttatt                                                                                                             >  1:916544/1‑150 (MQ=255)
                                               aCTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCt                                                                                                       >  1:1583832/1‑150 (MQ=255)
                                                     tCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                                 <  2:1064917/150‑1 (MQ=255)
                                                      cTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCa                                                                                                >  2:540143/1‑150 (MQ=255)
                                                       tgaCTCCTCCGTCAACGTGACTCAGGTTGGCTATGGTAACAACGCGACCGCTCATCAGTACTAATACATAATTTGTATTACAGAAACAGGGCTCAAGCCCTG‑‑TTATTTTTAGGGATCAGAATATGAATACGTTATTACTCATTGCGGCAc                                                                                               >  1:129128/3‑150 (MQ=255)
                                                       tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                               >  1:881316/1‑150 (MQ=255)
                                                       tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                               <  1:1525469/150‑1 (MQ=255)
                                                       tAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                               <  1:79566/150‑1 (MQ=255)
                                                        aaCTCCTCCGTCAACGTGACTCAGGGTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTTTTACAGAAACAGGGCGCAAGCCCTG‑‑GTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCaat                                                                                              >  1:1776678/1‑148 (MQ=255)
                                                          ctcctcCGTCAACGTGACTCAGGTTGGCTTTCGTAACAAGGCGACCGCTCATCAGTACTTATACATCATTTGTATTACAGAAACAGGGCGCAAGACCTG‑‑TTTTTTTTCGTGAGAACAATATGAAGACGTTATTACACCTTGCgcccattt                                                                                            >  2:1131922/1‑143 (MQ=255)
                                                              cccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATCCATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAg                                                                                        <  1:1731189/149‑1 (MQ=255)
                                                              accGTCCACGGGACCCCGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGACACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAg                                                                                        <  1:969687/149‑1 (MQ=255)
                                                               ccGTCAACGTGACTCAGGTTGGCTTTGGAAACAACGCGACCGCTCATCAGTACAAATACATCATTTGTATAACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGt                                                                                       <  1:825375/150‑1 (MQ=255)
                                                                  tCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAg                                                                                    <  2:1764132/150‑1 (MQ=255)
                                                                  tCAACGTGACTCAGGTAGGCTTTGGTAACAACGCGACCGCTCATCAGTTATAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAg                                                                                    <  2:459676/150‑1 (MQ=255)
                                                                   cAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGa                                                                                   <  2:1856983/150‑1 (MQ=255)
                                                                    aaCGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACATTATTAATCCTTGCGGCACTTTCCAGTCAGAt                                                                                  >  2:1774037/1‑150 (MQ=255)
                                                                         gACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCt                                                                             >  2:671828/1‑150 (MQ=255)
                                                                           cTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCttt                                                                           <  1:1064438/150‑1 (MQ=255)
                                                                           atctggTTGGCTTTGGTAAGAACGCGACCTCTCATCAGTACTTATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTGGGGAGAAGAATATGAATACGTTATTACTCCTTGCTGCACTTTCCAGTCAGATAACCttt                                                                           <  1:1496928/146‑1 (MQ=255)
                                                                            gcAGGTGGGCTTTGGTTACATCGCGACCGCTCTTCAGTACAAATACAACATTTGTATTACAGTAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGTGAGAAGAATATGAATTCGTTATTACTCCTTGCGGCACTTTCGAGTCAGATAACCTTTa                                                                          <  1:212348/149‑1 (MQ=255)
                                                                              tggTTGGCTTTGTTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACATGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAAt                                                                        <  1:725032/149‑1 (MQ=255)
                                                                               ggTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATa                                                                       <  2:939036/150‑1 (MQ=255)
                                                                               ggTTGGCTTTGGTAACAACGCGACAGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATa                                                                       >  2:1587273/1‑150 (MQ=255)
                                                                                 ttGGCTTTGGGAACCACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACACGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACg                                                                     <  2:221633/150‑1 (MQ=255)
                                                                                       ttGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAg                                                               <  2:2032725/150‑1 (MQ=255)
                                                                                        tGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGc                                                              <  2:1671639/150‑1 (MQ=255)
                                                                                           tAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAg                                                           >  2:1544854/1‑150 (MQ=255)
                                                                                             acaacGCGTACTCTCATCATTACTAATACTTCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAg                                                           <  2:180072/150‑1 (MQ=255)
                                                                                             acaacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAggg                                                         <  2:302322/150‑1 (MQ=255)
                                                                                             acaacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAggg                                                         >  1:983654/1‑150 (MQ=255)
                                                                                               aacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                       <  2:1483369/150‑1 (MQ=255)
                                                                                               aacGCGACCGCTCATCAGTACTAATAAATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                       >  1:1399400/1‑150 (MQ=255)
                                                                                                acGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAt                                                      <  2:1303157/150‑1 (MQ=255)
                                                                                                 cgcgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtg                                                     <  2:968958/150‑1 (MQ=255)
                                                                                                   cgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                                   <  2:569634/150‑1 (MQ=255)
                                                                                                   cgACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                                   <  2:1787311/150‑1 (MQ=255)
                                                                                                     aCCGGTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGACACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGta                                                 <  1:1613728/150‑1 (MQ=255)
                                                                                                       cGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGtata                                               <  2:1926733/150‑1 (MQ=255)
                                                                                                        gCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATAc                                              <  1:1988129/150‑1 (MQ=255)
                                                                                                             tcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCatta                                         <  1:2154218/150‑1 (MQ=255)
                                                                                                              caGTACTAATACATCATATGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCattat                                        <  1:463073/150‑1 (MQ=255)
                                                                                                                  aCTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCt                                    >  1:2010340/1‑150 (MQ=255)
                                                                                                                   cTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTg                                   >  1:514791/1‑150 (MQ=255)
                                                                                                                   cTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTg                                   >  1:1874602/1‑150 (MQ=255)
                                                                                                                    tAATACATCATTTGTATTACAGAAACAGGGAGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTGCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGa                                  <  2:1037743/150‑1 (MQ=255)
                                                                                                                     aaTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACACAGCAAGGGGATGTGTATACGATTATTCCTGaa                                 >  2:1324624/1‑150 (MQ=255)
                                                                                                                      aTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAg                                <  1:800818/150‑1 (MQ=255)
                                                                                                                          atcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCAc                            <  1:1817484/150‑1 (MQ=255)
                                                                                                                          atcatTTGTATTACAGAAACAGGGCGAAAGCCCTT‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCATTCAGATAACCTTTTATACGACCCAGCCAGGGGATGTGTATACCATTATTCCTGAAGTCAc                            <  2:1487815/150‑1 (MQ=255)
                                                                                                                           tcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATATCCTTTAATACGACCCATCAAGGGGATGTGTATAACATTATTCCTGAAGTCAct                           >  2:94092/1‑150 (MQ=255)
                                                                                                                             atTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCActct                         <  2:514791/150‑1 (MQ=255)
                                                                                                                                   cgaaaaaaaaCAGGGCTTGCGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTGTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCa                   >  2:933326/8‑150 (MQ=255)
                                                                                                                                   aTTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGTATGTGTATACCATTATTCCTGAAGTCACTCTTAATCa                   >  1:1551199/1‑150 (MQ=255)
                                                                                                                                   aTTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCa                   >  1:628017/1‑150 (MQ=255)
                                                                                                                                    ttACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCaa                  >  1:2136941/1‑150 (MQ=255)
                                                                                                                                      aaaaaaaCAGGGCTTGCGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc                >  2:147006/5‑150 (MQ=255)
                                                                                                                                      aCAGAAACAGGGCTAAAGCCATG‑‑TTTTTTTTCGGGAGAAGAATCTGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCAAGCAAGGGGTTGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc                <  1:959319/150‑1 (MQ=255)
                                                                                                                                         gAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTg             <  2:1283675/150‑1 (MQ=255)
                                                                                                                                           aaCACGGCGCAAGCCCTT‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTAGTACTCCTTGCGGCACTTTCCAGTCAGATACCCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc           <  2:350149/150‑1 (MQ=255)
                                                                                                                                            aCAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCt          <  2:881316/150‑1 (MQ=255)
                                                                                                                                             cAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATGCGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtg         <  2:1646460/150‑1 (MQ=255)
                                                                                                                                             cAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtg         >  1:631141/1‑150 (MQ=255)
                                                                                                                                                ggCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc      >  1:2073854/1‑150 (MQ=255)
                                                                                                                                                   gcAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTTTCTGTGCaga   >  1:194295/1‑150 (MQ=255)
                                                                                                                                                    cAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCagag  <  2:2020111/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAG  >  NC_000913/1104784‑1105081

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