Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,102,451 A→G intergenic (‑32/+435) csgD ← / ← b4518 CsgD/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,102,4510AG99.6% 814.5 / NA 247intergenic (‑32/+435)csgD/b4518CsgD/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (121/125);  total (121/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

GAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTA  >  NC_000913/1102361‑1102540
                                                                                          |                                                                                         
gaagaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGcc                                                                                         >  1:446812/1‑93 (MQ=255)
gaagaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGcc                                                                                         <  2:2330792/93‑1 (MQ=255)
gaagaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGcc                                                                                         <  2:3051073/93‑1 (MQ=255)
gaagaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGcc                                                                                         >  1:1702170/1‑93 (MQ=255)
  agaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTg                                                                                       <  2:5765591/93‑1 (MQ=255)
   gaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGa                                                                                      >  2:3922788/1‑93 (MQ=255)
   gaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGa                                                                                      >  2:4106818/1‑93 (MQ=255)
     tttAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACa                                                                                    >  1:2863489/1‑93 (MQ=255)
     tttAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACa                                                                                    >  1:993914/1‑93 (MQ=255)
       tAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGc                                                                                  <  1:1371290/93‑1 (MQ=255)
       tAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGc                                                                                  <  2:5967065/93‑1 (MQ=255)
        aGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCt                                                                                 >  1:150628/1‑93 (MQ=255)
           gATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGcc                                                                              <  2:1798378/93‑1 (MQ=255)
               aacaaTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTa                                                                          >  1:4735031/1‑93 (MQ=255)
                   ataatGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAAt                                                                      <  1:933118/93‑1 (MQ=255)
                    taatGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATa                                                                     <  2:1441264/93‑1 (MQ=255)
                     aatGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAg                                                                    >  1:2625996/1‑93 (MQ=255)
                      atGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGa                                                                   >  2:376340/1‑93 (MQ=255)
                      atGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGa                                                                   >  2:34155/1‑93 (MQ=255)
                      atGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGa                                                                   <  2:5612366/93‑1 (MQ=255)
                       tGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGaa                                                                  >  1:651506/1‑93 (MQ=255)
                          aTGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCa                                                               >  1:4655489/1‑93 (MQ=255)
                           tGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCAc                                                              >  1:276696/1‑93 (MQ=255)
                           tGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCAc                                                              <  2:337565/93‑1 (MQ=255)
                             aCCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCa                                                            <  2:4210030/93‑1 (MQ=255)
                               cATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGa                                                          >  2:222480/1‑93 (MQ=255)
                                aTGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTTACAGCTGCCTCTAAAATAGAAGCACCAGaa                                                         <  2:2827266/93‑1 (MQ=255)
                                   aaTACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTa                                                      <  2:1053961/93‑1 (MQ=255)
                                    aTACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTAc                                                     <  1:5227188/93‑1 (MQ=255)
                                    aTACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTAc                                                     <  2:6181350/93‑1 (MQ=255)
                                     tACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACt                                                    <  1:633837/93‑1 (MQ=255)
                                     tACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACt                                                    <  1:4132816/93‑1 (MQ=255)
                                     tACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACt                                                    <  2:3446558/93‑1 (MQ=255)
                                      aCTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTg                                                   <  1:5451319/93‑1 (MQ=255)
                                       cTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGa                                                  <  2:4436600/93‑1 (MQ=255)
                                        tATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTggc                                                 <  1:5309088/93‑3 (MQ=255)
                                        tATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGAc                                                 >  1:2952384/1‑93 (MQ=255)
                                        tATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGAc                                                 >  2:4565098/1‑93 (MQ=255)
                                         aTGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTAATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACa                                                >  2:5390995/1‑93 (MQ=255)
                                          tGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAg                                               >  1:3452330/1‑93 (MQ=255)
                                           ggACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGa                                              <  2:1594692/93‑1 (MQ=255)
                                            gACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGAt                                             >  2:3981680/1‑93 (MQ=255)
                                            gACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGAt                                             >  2:3567151/1‑93 (MQ=255)
                                             aCTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATg                                            <  1:1008694/93‑1 (MQ=255)
                                             aCTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATg                                            <  2:509161/93‑1 (MQ=255)
                                              cTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGt                                           <  2:4527479/93‑1 (MQ=255)
                                              cTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGt                                           <  1:5595240/93‑1 (MQ=255)
                                              cTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGt                                           >  2:6733424/1‑93 (MQ=255)
                                              cTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGt                                           >  1:3414627/1‑93 (MQ=255)
                                              cTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGt                                           <  2:5934917/93‑1 (MQ=255)
                                               ttCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGtt                                          >  1:5863130/1‑93 (MQ=255)
                                               ttCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGtt                                          >  2:4932830/1‑93 (MQ=255)
                                               ttCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGtt                                          >  2:2169324/1‑93 (MQ=255)
                                                tCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTg                                         >  2:4067276/1‑93 (MQ=255)
                                                 cATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGc                                        <  1:2533423/93‑1 (MQ=255)
                                                 cATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGc                                        >  1:2404734/1‑93 (MQ=255)
                                                  aTTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCa                                       >  1:6254287/1‑93 (MQ=255)
                                                   ttAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCAc                                      >  2:5807942/1‑93 (MQ=255)
                                                      aaCATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActgc                                   <  1:1620150/93‑1 (MQ=255)
                                                       aCATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActgct                                  >  2:706019/1‑93 (MQ=255)
                                                          tgatgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtg                               >  1:3221832/1‑93 (MQ=255)
                                                          tgatgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtg                               <  1:4002685/93‑1 (MQ=255)
                                                          tgatgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtg                               >  2:4482092/1‑93 (MQ=255)
                                                           gatgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtgt                              <  2:2081808/93‑1 (MQ=255)
                                                           gatgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtgt                              >  1:1692722/1‑93 (MQ=255)
                                                            atgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtgtg                             >  2:300670/1‑93 (MQ=255)
                                                               aaaCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtag                          <  1:2092662/93‑1 (MQ=255)
                                                               aaaCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtag                          >  1:3700852/1‑93 (MQ=255)
                                                                aaCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtagt                         >  2:5549301/1‑93 (MQ=255)
                                                                aaCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtagt                         >  1:5109770/1‑93 (MQ=255)
                                                                 aCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtagta                        <  1:1029122/93‑1 (MQ=255)
                                                                 aCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtagta                        >  2:2611057/1‑93 (MQ=255)
                                                                 aCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtagta                        <  1:1820289/93‑1 (MQ=255)
                                                                 aCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtagta                        >  1:2007893/1‑93 (MQ=255)
                                                                  ccccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaa                       <  2:2722145/93‑1 (MQ=255)
                                                                  ccccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaa                       >  1:5715075/1‑93 (MQ=255)
                                                                    ccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaata                     >  2:753078/1‑93 (MQ=255)
                                                                    ccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaata                     >  1:3689687/1‑93 (MQ=255)
                                                                    ccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaata                     >  2:1151486/1‑93 (MQ=255)
                                                                     cGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaataa                    >  1:4823807/1‑93 (MQ=255)
                                                                      gCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATaaa                   <  2:2810147/93‑1 (MQ=255)
                                                                       cTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAAt                  >  1:1136105/1‑93 (MQ=255)
                                                                        tttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATc                 <  2:5027910/93‑1 (MQ=255)
                                                                        tttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATc                 >  2:3258177/1‑93 (MQ=255)
                                                                        tttttttATTGATCGCACGCATGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATc                 <  2:3833430/93‑1 (MQ=255)
                                                                           ttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGc              >  1:160623/1‑93 (MQ=255)
                                                                           ttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTATGTGTAGTAAGAAATCAGc              <  1:6556835/93‑1 (MQ=255)
                                                                             ttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGccc            >  1:5018290/1‑93 (MQ=255)
                                                                             ttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGccc            <  2:5259052/93‑1 (MQ=255)
                                                                             ttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGccc            <  1:4297751/93‑1 (MQ=255)
                                                                             ttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGccc            <  2:6335304/93‑1 (MQ=255)
                                                                                 tGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTaaa        <  1:2061127/93‑1 (MQ=255)
                                                                                  gATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGGTGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAAt       <  2:4820112/93‑1 (MQ=255)
                                                                                  gATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAAt       <  2:1527974/93‑1 (MQ=255)
                                                                                   aTCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATg      >  1:5876377/1‑93 (MQ=255)
                                                                                   aTCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATg      <  1:1288972/93‑1 (MQ=255)
                                                                                    tCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATgg     >  2:4309674/1‑93 (MQ=255)
                                                                                    tCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATgg     >  2:3101331/1‑93 (MQ=255)
                                                                                      gCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGt   <  2:2060021/93‑1 (MQ=255)
                                                                                       cacgCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTa  <  1:3314452/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTA  >  NC_000913/1102361‑1102540

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