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,4510AG100.0% 628.2 / NA 190intergenic (‑32/+435)csgD/b4518CsgD/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (84/106);  total (84/106)

Only 100 of 193 total aligned reads displayed.

GAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTA  >  NC_000913/1102361‑1102540
                                                                                          |                                                                                         
gaagaGTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGcc                                                                                         <  1:3547588/93‑1 (MQ=255)
 aagaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCt                                                                                        <  2:4643751/93‑1 (MQ=255)
 aagaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCt                                                                                        <  2:4396475/93‑1 (MQ=255)
  agaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTg                                                                                       >  1:3960905/1‑93 (MQ=255)
  agaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTg                                                                                       >  2:1678927/1‑93 (MQ=255)
  agaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACGTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTg                                                                                       >  1:3010285/1‑93 (MQ=255)
   gaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGa                                                                                      >  2:2441612/1‑93 (MQ=255)
     tttAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACa                                                                                    >  2:3389126/1‑93 (MQ=255)
       tAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGc                                                                                  <  1:3985517/93‑1 (MQ=255)
        aGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCt                                                                                 <  2:3180380/93‑1 (MQ=255)
         gTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCtc                                                                                >  1:1443440/1‑92 (MQ=255)
         atGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTg                                                                                <  1:1149950/92‑1 (MQ=255)
            aTCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCct                                                                             >  2:4617164/1‑93 (MQ=255)
             tCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCctc                                                                            >  2:3734516/1‑93 (MQ=255)
             tCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCctc                                                                            >  2:4076323/1‑93 (MQ=255)
             tCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCctc                                                                            <  1:656327/93‑1 (MQ=255)
              caacaaTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCctct                                                                           >  2:1991781/1‑93 (MQ=255)
              caacaaTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCctct                                                                           <  1:4282746/93‑1 (MQ=255)
              caacaaTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCctct                                                                           <  1:4368049/93‑1 (MQ=255)
              caacaaTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCctct                                                                           <  2:3232921/93‑1 (MQ=255)
               aacaaTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTa                                                                          >  1:1720472/1‑93 (MQ=255)
                  cataatGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTaaaa                                                                       <  1:3973551/92‑1 (MQ=255)
                  aataatGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTaaaa                                                                       <  2:4575386/93‑1 (MQ=255)
                    taatGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATa                                                                     <  2:4633868/93‑1 (MQ=255)
                    taatGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATa                                                                     <  1:2413082/93‑1 (MQ=255)
                    taatGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATa                                                                     >  2:3426568/1‑93 (MQ=255)
                     aatGTATGACCATGAATACTATAGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAg                                                                    >  2:2991755/1‑93 (MQ=255)
                      atGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGa                                                                   <  1:4476528/93‑1 (MQ=255)
                      atGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGa                                                                   <  1:4475460/93‑1 (MQ=255)
                       tGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGaa                                                                  >  2:3820937/1‑93 (MQ=255)
                         tATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGc                                                                <  1:2806018/93‑1 (MQ=255)
                         tATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGc                                                                >  1:4106614/1‑93 (MQ=255)
                         tATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGc                                                                <  2:3108208/93‑1 (MQ=255)
                             aCCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCa                                                            <  1:1198610/93‑1 (MQ=255)
                              ccATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAg                                                           >  1:2942726/1‑93 (MQ=255)
                              ccATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAg                                                           >  2:3693122/1‑93 (MQ=255)
                               cATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGa                                                          >  2:4733741/1‑93 (MQ=255)
                                aTGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGaa                                                         >  1:2076163/1‑93 (MQ=255)
                                  gAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAga                                                       >  2:3146704/1‑92 (MQ=255)
                                    aTACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTAc                                                     <  1:3573642/93‑1 (MQ=255)
                                    aTACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTAc                                                     >  2:1415134/1‑93 (MQ=255)
                                    aTACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTAc                                                     <  2:544362/93‑1 (MQ=255)
                                     tACTATGGACTTCATTAAACATGATGAACCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACt                                                    <  2:4489265/93‑1 (MQ=255)
                                     tACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACt                                                    <  2:3785296/93‑1 (MQ=255)
                                     aaCTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACt                                                    <  1:3716534/92‑1 (MQ=255)
                                      aCTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTg                                                   <  2:2534461/93‑1 (MQ=255)
                                      aCTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTg                                                   <  1:2226234/93‑1 (MQ=255)
                                        tATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGAc                                                 >  1:592452/1‑93 (MQ=255)
                                        tATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGAc                                                 <  2:2559610/93‑1 (MQ=255)
                                        tATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGAc                                                 >  1:3505874/1‑93 (MQ=255)
                                          tGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAg                                               >  2:2295538/1‑93 (MQ=255)
                                           ggACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGa                                              <  1:1389404/93‑1 (MQ=255)
                                            gACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGAt                                             <  1:1517070/93‑1 (MQ=255)
                                            gACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGAt                                             <  2:2960075/93‑1 (MQ=255)
                                             aCTTCATTAAACATGATGAAGCCCCGCTTTTTTTATTGATCGCACGCCTGCCAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATg                                            <  1:3525012/93‑1 (MQ=255)
                                             aCTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATg                                            >  1:1549216/1‑93 (MQ=255)
                                               ttCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGtt                                          >  2:1572090/1‑93 (MQ=255)
                                                tCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTg                                         >  1:4493789/1‑93 (MQ=255)
                                                tCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTg                                         <  2:4056472/93‑1 (MQ=255)
                                                  aTTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCa                                       <  2:3924791/93‑1 (MQ=255)
                                                      aaCATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActgc                                   <  1:2050061/93‑1 (MQ=255)
                                                      aaCATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActgc                                   >  2:3519539/1‑93 (MQ=255)
                                                      aaCATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActgc                                   <  1:433497/93‑1 (MQ=255)
                                                      aaCATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActgc                                   <  1:1611766/93‑1 (MQ=255)
                                                      aaCATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActgc                                   <  2:2359486/93‑1 (MQ=255)
                                                       aCATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActgct                                  <  1:140795/93‑1 (MQ=255)
                                                         atgatgAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgt                                >  2:1686363/1‑93 (MQ=255)
                                                          tgatgaAACCCCGCTTTTTTTATTGATCGCACGCATGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCGCTGCtgtg                               <  2:4153946/93‑1 (MQ=255)
                                                           gatgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtgt                              <  2:1212006/93‑1 (MQ=255)
                                                           gatgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtgt                              <  1:3341270/93‑1 (MQ=255)
                                                            atgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtgtg                             <  1:2008473/93‑1 (MQ=255)
                                                              gaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgta                           >  2:1285922/1‑93 (MQ=255)
                                                              gaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgta                           >  2:364350/1‑93 (MQ=255)
                                                               aaaCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtag                          >  1:3658951/1‑93 (MQ=255)
                                                                  ccccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaa                       <  2:1019247/93‑1 (MQ=255)
                                                                   cccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaat                      <  2:169740/93‑1 (MQ=255)
                                                                    ccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaata                     >  2:812238/1‑93 (MQ=255)
                                                                    ccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaata                     >  1:1013884/1‑93 (MQ=255)
                                                                     cGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaataa                    >  2:3674550/1‑93 (MQ=255)
                                                                       cTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAAt                  <  1:2599502/93‑1 (MQ=255)
                                                                        tttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATc                 <  1:107029/93‑1 (MQ=255)
                                                                        tttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATc                 >  1:2333248/1‑93 (MQ=255)
                                                                         ttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCa                <  1:1655082/93‑1 (MQ=255)
                                                                         ttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCa                <  1:1923864/93‑1 (MQ=255)
                                                                         ttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCa                <  2:1765957/93‑1 (MQ=255)
                                                                         ttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCa                <  1:778494/93‑1 (MQ=255)
                                                                         ttttttATTGATAGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAAtta                >  2:2293788/1‑91 (MQ=255)
                                                                             ttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGccc            >  2:2126299/1‑93 (MQ=255)
                                                                              tATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCt           <  2:181641/93‑1 (MQ=255)
                                                                               aTTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTa          >  2:4273815/1‑93 (MQ=255)
                                                                                ttGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTaa         >  2:33008/1‑93 (MQ=255)
                                                                                 tGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTaaa        <  1:4345823/93‑1 (MQ=255)
                                                                                 tGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTaaa        <  1:4019418/93‑1 (MQ=255)
                                                                                  gATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAAt       <  1:2565715/93‑1 (MQ=255)
                                                                                  gATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAAt       <  2:689137/93‑1 (MQ=255)
                                                                                    tCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATgg     <  1:2588468/93‑1 (MQ=255)
                                                                                      gCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATggg  >  2:1161746/1‑92 (MQ=255)
                                                                                      gCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGt   <  1:4617164/93‑1 (MQ=255)
                                                                                      gCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGt   <  1:2842954/93‑1 (MQ=255)
                                                                                       cacgCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTa  <  2:4527989/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTA  >  NC_000913/1102361‑1102540

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