Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,104,949 C→T intergenic (+47/‑12) 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,9490CT98.8% 284.5 / ‑3.8 83intergenic (+47/‑12)csgA/csgCcurlin subunit, amyloid curli fibers, cryptic/curli assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (38/44);  minor base . (0/1);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAG  >  NC_000913/1104815‑1105079
                                                                                                                                      |                                                                                                                                  
tgctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGga                                                                                                                                >  2:303797/1‑139 (MQ=255)
 gctgtAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgag                                                                                                                               <  1:351179/139‑1 (MQ=255)
 gctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgag                                                                                                                               <  2:1162306/139‑1 (MQ=255)
 gctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgag                                                                                                                               <  1:303797/139‑1 (MQ=255)
  ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGgaga                                                                                                                              <  1:266643/139‑1 (MQ=255)
  ctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTGGGagaa                                                                                                                             <  2:607419/139‑1 (MQ=255)
     cAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGagaaga                                                                                                                           <  2:971469/139‑1 (MQ=255)
         tGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAatat                                                                                                                       >  2:1150478/1‑139 (MQ=255)
              aGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATa                                                                                                                  <  1:473292/139‑1 (MQ=255)
                         aaCTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTc                                                                                                       >  2:1080676/1‑139 (MQ=255)
                              ctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGc                                                                                                  >  1:880874/1‑139 (MQ=255)
                                 cGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGc                                                                                               <  2:832010/139‑1 (MQ=255)
                                   tCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCAc                                                                                             >  1:55777/1‑139 (MQ=255)
                                       cGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTc                                                                                         >  2:1025167/1‑139 (MQ=255)
                                            cTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTc                                                                                    <  2:55777/139‑1 (MQ=255)
                                                 gTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATa                                                                               <  1:1025167/139‑1 (MQ=255)
                                                           gTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATa                                                                     >  2:302930/1‑139 (MQ=255)
                                                                     gACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                                >  2:508799/1‑134 (MQ=255)
                                                                     gACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                                <  1:508799/134‑1 (MQ=255)
                                                                           tcatcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                     <  2:1162707/139‑1 (MQ=255)
                                                                           tcatcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGa                                                     <  1:1080676/139‑1 (MQ=255)
                                                                            catcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAt                                                    <  1:302930/139‑1 (MQ=255)
                                                                            catcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAt                                                    <  2:880874/139‑1 (MQ=255)
                                                                                 gTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGta                                               <  2:390600/139‑1 (MQ=255)
                                                                                    cTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGagaaga                                                                                                                           <  1:204065/60‑1 (MQ=255)
                                                                                    cTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGagaaga                                                                                                                           >  2:204065/1‑60 (MQ=255)
                                                                                       aTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGagaa                                                                                                                             <  2:584846/55‑1 (MQ=255)
                                                                                       aTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGagaa                                                                                                                             >  1:584846/1‑55 (MQ=255)
                                                                                        tACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAGGGGGATGTGTATACCatt                                        <  1:530485/139‑1 (MQ=255)
                                                                                         aCATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGaa                                                                                                                    <  1:566046/62‑1 (MQ=255)
                                                                                         aCATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGaa                                                                                                                    >  2:566046/1‑62 (MQ=255)
                                                                                          catcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCa                                                                                       >  2:601025/1‑90 (MQ=255)
                                                                                          catcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCa                                                                                       <  1:601025/90‑1 (MQ=255)
                                                                                            tcatTTGTATTACAGCAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGAt                                                                                >  2:827557/1‑95 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTTGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAggg                                                       <  2:887280/120‑1 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTGTTACTCCTTGCGGCACTTTCCAGTCAGAt                                                                                <  1:827557/95‑1 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACt                                                                                                        <  1:777786/71‑1 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACt                                                                                                        >  2:777786/1‑71 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGAt                                                                                >  2:731858/1‑95 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGAt                                                                                <  1:731858/95‑1 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAAcc                                                                            <  2:1219605/99‑1 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAAcc                                                                            >  1:1219605/1‑99 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCtt                                                                          <  1:812809/101‑1 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCtt                                                                          >  2:812809/1‑101 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                                <  2:1036674/111‑1 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                                >  1:1036674/1‑111 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAggg                                                       >  1:887280/1‑120 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                    >  2:1051565/1‑139 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                    >  1:14661/1‑139 (MQ=255)
                                                                                            tcatTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTc                                    >  1:547781/1‑139 (MQ=255)
                                                                                             catTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGa                                                                  <  2:508332/108‑1 (MQ=255)
                                                                                             catTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGa                                                                  >  1:508332/1‑108 (MQ=255)
                                                                                               tttGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAt                                                    <  1:449050/120‑1 (MQ=255)
                                                                                               tttGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAt                                                    >  2:449050/1‑120 (MQ=255)
                                                                                               tttGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTg                                 >  2:1021794/1‑139 (MQ=255)
                                                                                                 tGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTcc                                   >  1:1043842/1‑135 (MQ=255)
                                                                                                 tGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTcc                                   <  2:1043842/135‑1 (MQ=255)
                                                                                                  gTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGGGAAGAATATGAATACGTTATTACCCCTTGCGGCACTTTCCAGTCAGGTAACCTTTAATAAGACCCAGCAGGGGGGTGTGTATACCATTATTCCTGaaa                              >  2:530485/1‑138 (MQ=255)
                                                                                                  gTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                                <  1:882628/105‑1 (MQ=255)
                                                                                                  gTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                                >  2:882628/1‑105 (MQ=255)
                                                                                                  gTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCat                                         <  1:1241380/128‑1 (MQ=255)
                                                                                                  gTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCat                                         >  2:1241380/1‑128 (MQ=255)
                                                                                                   tATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGt                             <  2:14661/139‑1 (MQ=255)
                                                                                                   tATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGt                             <  2:775542/139‑1 (MQ=255)
                                                                                                   tATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGt                             <  2:547781/139‑1 (MQ=255)
                                                                                                   tATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGt                             <  1:1051565/139‑1 (MQ=255)
                                                                                                     ttACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGaa                                                                                                                    <  2:841304/50‑1 (MQ=255)
                                                                                                     ttACAGAAACAGGGCGCAAGACCTGTTTTTTTTTGGGAGAAGAATATGaa                                                                                                                    >  1:841304/1‑50 (MQ=255)
                                                                                                         aGAAACAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCActct                       >  2:1231245/1‑139 (MQ=255)
                                                                                                           aaaCAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTa                     >  2:463932/1‑139 (MQ=255)
                                                                                                             aCAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCt                                  <  1:1067411/124‑1 (MQ=255)
                                                                                                             aCAGGGCGCAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCt                                  >  2:1067411/1‑124 (MQ=255)
                                                                                                                  gcgcAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgggg                                                      <  1:1110545/99‑1 (MQ=255)
                                                                                                                  gcgcAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAgggg                                                      >  2:1110545/1‑99 (MQ=255)
                                                                                                                  gcgcAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc              >  1:181843/1‑139 (MQ=255)
                                                                                                                  gcgcAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc              >  1:74409/1‑139 (MQ=255)
                                                                                                                  gcgcAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc              >  2:865239/1‑139 (MQ=255)
                                                                                                                   cgcAAGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCt             <  2:74409/139‑1 (MQ=255)
                                                                                                                       aGCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTc         >  1:37383/1‑139 (MQ=255)
                                                                                                                        gCCCTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCt        <  1:921038/139‑1 (MQ=255)
                                                                                                                           cTGTTTTTTTTTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCtgtg     <  1:1021794/139‑1 (MQ=255)
                                                                                                                              ttttttttTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCag  <  2:37383/139‑1 (MQ=255)
                                                                                                                              ttttttttTGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCag  <  1:1231245/139‑1 (MQ=255)
                                                                                                                                   ttttGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCAct                         <  1:537656/111‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
TGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAG  >  NC_000913/1104815‑1105079

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