Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,434,402 T→C 58.8% intergenic (‑178/+556) ynaE ← / ← ttcC cold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,434,4020TC58.8% 29.7 / 96.3 90intergenic (‑178/+556)ynaE/ttcCcold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (20/17);  new base C (33/20);  total (53/37)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

GTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGG  >  NC_000913/1434303‑1434501
                                                                                                   |                                                                                                   
gTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGt                                                                                                     <  1:2499792/100‑1 (MQ=255)
gTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGt                                                                                                     >  2:2040829/1‑100 (MQ=40)
  tttGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTc                                                                                                   <  1:344484/100‑1 (MQ=40)
  tttGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTc                                                                                                   <  1:3261411/100‑1 (MQ=255)
  tttGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTc                                                                                                   >  2:1883901/1‑100 (MQ=255)
   ttGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCa                                                                                                  <  1:1138725/100‑1 (MQ=40)
    tGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCat                                                                                                 <  2:304827/100‑1 (MQ=255)
     gTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCata                                                                                                <  2:2578897/100‑1 (MQ=255)
       cTGGGTTATTGGATATGCCCATGCGTAATATTATTAGAGCGTGAACAACAATGCTGTATATGATAAGGAGCGAAATAATTATAAGTATGCTGCTCataca                                                                                              >  2:1555781/1‑98 (MQ=9)
        tGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCAtata                                                                                              <  2:270357/99‑1 (MQ=255)
         gggTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGc                                                                                            >  2:3496304/1‑100 (MQ=255)
         gggTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGc                                                                                            <  2:3947179/100‑1 (MQ=255)
         gggTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGc                                                                                            >  1:711040/1‑100 (MQ=255)
           gTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCt                                                                                          >  2:1293292/1‑100 (MQ=255)
             tATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCg                                                                                        <  1:1210817/100‑1 (MQ=255)
             tATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCg                                                                                        >  1:1208766/1‑100 (MQ=255)
             tATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCg                                                                                        <  1:2328507/100‑1 (MQ=255)
                 ggATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAACCTCAAAtt                                                                                    >  2:453057/1‑100 (MQ=255)
                 ggATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAAtt                                                                                    <  1:3452137/100‑1 (MQ=255)
                   atatGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAg                                                                                  <  1:3496326/100‑1 (MQ=255)
                    tatGCCAATCGCTAATTTTATTAGAGCATGACGAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGt                                                                                 <  2:3434244/100‑1 (MQ=255)
                    tatGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTTAATATTATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGt                                                                                 <  2:3460423/100‑1 (MQ=255)
                    tatGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGt                                                                                 <  2:3460432/100‑1 (MQ=40)
                      tGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTaa                                                                               >  1:1345701/1‑100 (MQ=255)
                      tGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTa                                                                                <  1:2233416/99‑1 (MQ=255)
                       gCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAAt                                                                              <  2:1153132/100‑1 (MQ=40)
                       gCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAAt                                                                              >  1:1095495/1‑100 (MQ=255)
                       gCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAAt                                                                              >  2:3316188/1‑100 (MQ=255)
                       gCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAAt                                                                              >  2:3813223/1‑100 (MQ=255)
                       gCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAAt                                                                              <  2:745418/100‑1 (MQ=255)
                        ccAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAAtg                                                                             >  2:2960888/1‑100 (MQ=255)
                        ccAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAAtg                                                                             >  2:996452/1‑100 (MQ=255)
                           aTCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATgtgt                                                                          <  2:269278/100‑1 (MQ=255)
                              ccTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTtat                                                                       >  2:2151440/1‑100 (MQ=255)
                              ccTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTtat                                                                       >  2:1593621/1‑100 (MQ=255)
                               cTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCGGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTtata                                                                      >  1:3843387/1‑100 (MQ=255)
                                tAATTTTATTAGAGCATGACTAAATATGCTGAATATGATAAGGAGCGAAGTGATTATCTGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTtat                                                                       <  2:2809819/98‑1 (MQ=255)
                                tAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTtatat                                                                     >  1:326585/1‑100 (MQ=255)
                                tAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTtata                                                                      <  2:1756334/99‑1 (MQ=40)
                                   ttttATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATg                                                                  <  2:1124373/100‑1 (MQ=255)
                                     ttattaGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGat                                                                <  2:1023655/100‑1 (MQ=255)
                                      tattaGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGata                                                               <  2:1229646/100‑1 (MQ=255)
                                       attaGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGatat                                                              >  1:57810/1‑100 (MQ=255)
                                       attaGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGatat                                                              >  1:1612053/1‑100 (MQ=255)
                                         taGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAg                                                            >  1:1340629/1‑100 (MQ=255)
                                          agagCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTAACAGTCTGCTGCGCATATAGCCTCGACTTAGTAATGTGTTATATATGATATAGt                                                           >  2:2981335/1‑100 (MQ=14)
                                            agCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTg                                                         <  1:2221974/100‑1 (MQ=40)
                                            aaaataaatAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGtt                                                          <  2:2435135/91‑1 (MQ=255)
                                             gCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGa                                                        >  2:3021595/1‑100 (MQ=38)
                                             gCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTGATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGa                                                        >  2:2642953/1‑100 (MQ=255)
                                             gCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGa                                                        >  1:3041472/1‑100 (MQ=255)
                                             gCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTAATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGa                                                        >  2:2643445/1‑100 (MQ=255)
                                              cATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGAc                                                       >  2:3507541/1‑100 (MQ=25)
                                                 gACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAAt                                                    <  1:3109556/100‑1 (MQ=38)
                                                   cTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAAttt                                                  >  2:261768/1‑100 (MQ=38)
                                                   cTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAAttt                                                  >  1:3997310/1‑100 (MQ=255)
                                                   cTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAAttt                                                  >  1:2760136/1‑100 (MQ=255)
                                                    tAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAAtttt                                                 >  2:2050110/1‑100 (MQ=38)
                                                    tAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAAtttt                                                 >  2:1100709/1‑100 (MQ=38)
                                                    tAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAAtttt                                                 >  1:2373199/1‑100 (MQ=255)
                                                     aaaaaTGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAAttttt                                                >  1:1058190/1‑100 (MQ=38)
                                                      aaaaTGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTa                                               <  1:3609747/100‑1 (MQ=255)
                                                           gCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATcttg                                          >  1:2217643/1‑97 (MQ=255)
                                                            cTGAATATGATAAGGCGCGAAGTGATTATCAGTATGCTGCTCATATTGCCTCGAATTAGTAATGTGTGATATATGATATAGTTGACAATTTTTATCTTgg                                         <  2:3861801/100‑1 (MQ=11)
                                                            cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                         >  2:2683331/1‑100 (MQ=18)
                                                            cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATcttag                                         >  2:2496754/1‑96 (MQ=255)
                                                             tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                        <  2:3286316/100‑1 (MQ=12)
                                                             tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                        >  1:1676227/1‑100 (MQ=255)
                                                             tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                        >  1:2873583/1‑100 (MQ=255)
                                                              gAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgt                                       >  1:1332736/1‑100 (MQ=18)
                                                              gAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgt                                       >  1:3205690/1‑100 (MQ=25)
                                                              gAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgt                                       >  1:641831/1‑100 (MQ=255)
                                                               aaTATGATAAGGAGCGAAGTGATTATCAGTATTCTGTTTATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTCGACAATTTTTATCTTGGgtg                                      <  2:2525473/100‑1 (MQ=11)
                                                                atatGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAAGTAGAAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgtgt                                     >  2:1999815/1‑100 (MQ=255)
                                                                   tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                  >  1:3607025/1‑100 (MQ=18)
                                                                   tGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                  <  1:3525301/100‑1 (MQ=255)
                                                                   tGATAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                  <  2:2974295/100‑1 (MQ=9)
                                                                      tAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTaa                               >  2:1137768/1‑100 (MQ=18)
                                                                      tAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTaa                               >  2:3839308/1‑100 (MQ=17)
                                                                      tAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTaa                               <  1:2180037/100‑1 (MQ=255)
                                                                      tAAGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTaa                               <  1:1416453/100‑1 (MQ=9)
                                                                       aaGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTATTaaa                              >  2:2755112/1‑100 (MQ=12)
                                                                        aGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAg                             >  2:1387635/1‑100 (MQ=255)
                                                                        aGGAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAg                             <  2:1535810/100‑1 (MQ=12)
                                                                         ggAGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGCATTCTATATGATATAATTGACAATTTTTATCTTGAGTGTTCTTaaaat                            >  2:429972/1‑98 (MQ=11)
                                                                          naGCGAAGTGATTATCAGTATTCTTTTCATATAGCCTCGAATTAGTAATGTGTTATTTATGAATTAGTTGACAATTTTTATCTTGGGTGTTTTTATAGtt                           >  2:3980886/2‑100 (MQ=11)
                                                                           aGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGtt                           <  1:1703362/99‑1 (MQ=9)
                                                                           aGCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTc                          >  1:2073268/1‑100 (MQ=9)
                                                                            gCGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                         >  1:2826419/1‑100 (MQ=2)
                                                                             cGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCGt                        >  2:2988026/1‑100 (MQ=0)
                                                                             cGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCGt                        >  1:1220436/1‑100 (MQ=2)
                                                                             cGAAGTGATTATCAGTATGCTGCTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCGt                        <  1:1117275/100‑1 (MQ=0)
                                                                                                   ttCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCTgg  >  2:1745174/1‑100 (MQ=17)
                                                                                                   |                                                                                                   
GTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGG  >  NC_000913/1434303‑1434501

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