Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,104,872 T→G F142L (TTT→TTG csgA → curlin subunit, amyloid curli fibers, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,104,8720TG100.0% 252.9 / NA 67F142L (TTT→TTGcsgAcurlin subunit, amyloid curli fibers, cryptic
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (34/33);  total (34/33)

CCAACGTGGCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATAC  >  NC_000913/1104725‑1105013
                                                                                                                                                   |                                                                                                                                             
ccAACGTGGCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGt                                                                                                                                            <  1:390437/151‑1 (MQ=255)
 cAACGTGGCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGt                                                                                                                                            >  2:194182/1‑150 (MQ=255)
  aaCGTGGCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTa                                                                                                                                           >  1:547539/1‑150 (MQ=255)
    cGTGGCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCTAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTaaca                                                                                                                                        <  2:883447/151‑1 (MQ=255)
     gTGGCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTaacaa                                                                                                                                       <  2:653463/151‑1 (MQ=255)
       ggCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTaacaac                                                                                                                                      <  2:465522/150‑1 (MQ=255)
       ggCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTaacaac                                                                                                                                      <  1:528928/150‑1 (MQ=255)
       ggCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAAcg                                                                                                                                     >  1:535579/1‑151 (MQ=255)
        gCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAAcgc                                                                                                                                    <  1:780443/151‑1 (MQ=255)
        gCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAAcgc                                                                                                                                    >  1:33198/1‑151 (MQ=255)
          ttCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAAcgcg                                                                                                                                   >  1:860313/1‑150 (MQ=255)
          ttCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAAcgc                                                                                                                                    >  1:258737/1‑149 (MQ=255)
                aaCAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCtc                                                                                                                            >  2:235672/1‑151 (MQ=255)
                aaCAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCtc                                                                                                                            >  1:832813/1‑151 (MQ=255)
                  cAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCtcat                                                                                                                          <  2:547162/151‑1 (MQ=255)
                   aGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCtcatc                                                                                                                         >  2:24536/1‑151 (MQ=255)
                    gcgcTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCtcatca                                                                                                                        <  2:590887/151‑1 (MQ=255)
                    gcgcTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCtcatc                                                                                                                         >  2:278669/1‑150 (MQ=255)
                     cgcTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCtcatc                                                                                                                         >  2:421647/1‑149 (MQ=255)
                     cgcTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCtcatc                                                                                                                         >  1:109279/1‑149 (MQ=255)
                         aCTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTAc                                                                                                                    <  2:621954/150‑1 (MQ=255)
                              tGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATa                                                                                                               >  1:802636/1‑150 (MQ=255)
                               gATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATAc                                                                                                              <  2:673337/150‑1 (MQ=255)
                                   aGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATAcatcat                                                                                                         <  1:15472/151‑1 (MQ=255)
                                    gTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCAtt                                                                                                        >  1:563180/1‑151 (MQ=255)
                                    gTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCAtt                                                                                                        <  1:53489/151‑1 (MQ=255)
                                    gTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCAtt                                                                                                        >  2:527750/1‑151 (MQ=255)
                                        aaCGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTa                                                                                                    <  1:329710/151‑1 (MQ=255)
                                          cGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTAtt                                                                                                  >  1:666347/1‑151 (MQ=255)
                                            gCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTAc                                                                                                >  1:853678/1‑151 (MQ=255)
                                                  aTTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAAc                                                                                          <  1:471383/151‑1 (MQ=255)
                                                    tCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAAc                                                                                          >  2:769528/1‑149 (MQ=255)
                                                    tCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAAc                                                                                          >  1:825542/1‑149 (MQ=255)
                                                       gAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGGCTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAggg                                                                                      <  2:49525/150‑1 (MQ=255)
                                                       gAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAg                                                                                        <  2:167410/148‑1 (MQ=255)
                                                        aaaTGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGgcg                                                                                    >  1:254798/1‑151 (MQ=255)
                                                        aaaTGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGgcg                                                                                    >  1:826067/1‑151 (MQ=255)
                                                         aaTGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGAAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGgcgc                                                                                   >  1:577635/1‑151 (MQ=255)
                                                                     aCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGtttt                                                                       <  2:440649/151‑1 (MQ=255)
                                                                       aGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGtttttt                                                                     >  2:382729/1‑151 (MQ=255)
                                                                       aGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGtttttt                                                                     >  1:559323/1‑151 (MQ=255)
                                                                          tCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTc                                                                  >  2:300681/1‑151 (MQ=255)
                                                                          tCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTc                                                                  >  2:857076/1‑151 (MQ=255)
                                                                              tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCggg                                                               >  2:49852/1‑150 (MQ=255)
                                                                                gtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCggg                                                               >  2:128980/1‑148 (MQ=255)
                                                                                  ggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGagaag                                                          <  2:860294/151‑1 (MQ=255)
                                                                                    cAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAgaagaa                                                        <  1:218570/151‑1 (MQ=255)
                                                                                     aaCGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAat                                                       >  1:835791/1‑151 (MQ=255)
                                                                                         gTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGa                                                   <  2:333981/151‑1 (MQ=255)
                                                                                         gTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGa                                                   <  1:611209/151‑1 (MQ=255)
                                                                                              gcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATa                                                <  2:56099/149‑1 (MQ=255)
                                                                                              gcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATa                                                <  1:278675/149‑1 (MQ=255)
                                                                                              gcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATAc                                               <  1:250113/150‑1 (MQ=255)
                                                                                              gcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATAc                                               <  1:313607/150‑1 (MQ=255)
                                                                                                aGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGtt                                            >  2:250346/1‑151 (MQ=255)
                                                                                                 gTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGtta                                           <  1:696199/151‑1 (MQ=255)
                                                                                                    gACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGttatt                                         >  2:671028/1‑150 (MQ=255)
                                                                                                             gCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTg                                 <  2:610513/149‑1 (MQ=255)
                                                                                                             gCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGc                                <  2:563165/150‑1 (MQ=255)
                                                                                                                   aaCTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                             <  1:915611/147‑1 (MQ=255)
                                                                                                                   aaCTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCa                            <  2:312258/148‑1 (MQ=255)
                                                                                                                   aaCTCCTCCGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                           <  1:192360/149‑1 (MQ=255)
                                                                                                                       cctccGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTcc                      <  2:832795/150‑1 (MQ=255)
                                                                                                                       cctccGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTcc                      >  1:72202/1‑150 (MQ=255)
                                                                                                                        ctccGTCAACGTGACTCAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTcc                      >  1:838678/1‑149 (MQ=255)
                                                                                                                                        cAGGTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAAt    <  1:358089/151‑1 (MQ=255)
                                                                                                                                          ggTTGGCTTGGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATAc  >  2:500451/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CCAACGTGGCTTCGGTAACAGCGCTACTCTTGATCAGTGGAACGGCAAAAATTCTGAAATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATAC  >  NC_000913/1104725‑1105013

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