Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,317 A→C 15.2% N21H (AAT→CAT)  yffL → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3170AC15.2% 493.2 / 75.2 189N21H (AAT→CAT) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (72/88);  new base C (16/13);  total (88/101)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 189 total aligned reads displayed.

AGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTT  >  NC_000913/2560220‑2560415
                                                                                                 |                                                                                                  
     aaTAAATTATTCAATATCAAGAATTTAAACCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACTCAAATAAAAGCCCCCGGAAtgctg                                                                                              <  2:805020/99‑1 (MQ=255)
     aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAGAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                             >  1:2079165/1‑100 (MQ=255)
     aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                             >  1:565743/1‑100 (MQ=255)
     aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                             <  2:1607962/100‑1 (MQ=255)
      aTAAATTATTCAATATGAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTAGAAAATAGGATCAACAAAAATAAAAGCACCCGGAAtgttgta                                                                                            <  2:1861651/100‑3 (MQ=255)
      aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                            <  2:3704315/100‑1 (MQ=255)
      aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACCAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctg                                                                                              <  2:2827023/98‑1 (MQ=255)
      aTAAATTATTAAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                            <  2:2720759/100‑1 (MQ=255)
       tAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                                                                                           >  1:3254847/1‑100 (MQ=255)
         aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGTAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                                                                                           >  2:2500773/1‑98 (MQ=255)
            tattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGc                                                                                      <  1:1148066/100‑1 (MQ=255)
              ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                    <  1:30767/100‑1 (MQ=255)
               tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGTATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgg                                                                                   <  2:15904/100‑1 (MQ=255)
                cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                  <  1:3137197/100‑1 (MQ=38)
                cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                  >  2:2075335/1‑100 (MQ=255)
                cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                  <  1:2438152/100‑1 (MQ=255)
                 aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                 <  2:1751561/100‑1 (MQ=38)
                  atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                <  2:1927339/100‑1 (MQ=38)
                    atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGct                                                                              <  2:2582786/100‑2 (MQ=38)
                    atCAAAAACTTAAGTCCGCGTTTGCAACGAAAGATCCAGGGTTGGAAAATAGGATCAACACAAATAAAAGCCCCCTGCATGTTGCAGCATTCCggggg                                                                                <  2:2393602/98‑1 (MQ=25)
                     tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                             <  2:302222/100‑1 (MQ=255)
                     tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                             >  1:2161839/1‑100 (MQ=255)
                      caacaaCTTCAATCCCTCTTTGCATCTTATGATCCTACGTCCACAAACAGGATCACCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                             <  2:4080503/99‑1 (MQ=255)
                      cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGACATGTTGCAGCATTCCGGGGGcttt                                                                            <  1:2107368/100‑4 (MQ=37)
                          aaCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                        <  1:1516337/100‑1 (MQ=255)
                          aaCTTAAACCCCTATTTGCAACGAAAGATCACGAGTCCGAAAATAGGATCACCCCAACTAACAGCCCCCGGCATGCTGCAACATGCCGGGGGCGGTATga                                                                        <  1:1220780/100‑1 (MQ=255)
                            tttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                      <  1:3839942/99‑1 (MQ=255)
                            ctcagaTCCATGTTTGCAACGAAACATCCATAGCTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTCTga                                                                        <  2:3078636/93‑1 (MQ=255)
                            cTTAAATCCATGTTTGTAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                       <  1:1853100/99‑1 (MQ=255)
                            cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                      <  1:1614272/100‑1 (MQ=255)
                            cTTAAATCCATGTTTGCAACGAAAGATCAAGAGTTAGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAAATGCCGGGGGCGGTATgaga                                                                      >  2:553095/1‑100 (MQ=255)
                                  tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                >  2:2293038/1‑100 (MQ=255)
                                  tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                >  1:3484669/1‑100 (MQ=255)
                                   ccATGTTTGCAACGTAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                               <  1:2500787/100‑1 (MQ=255)
                                   ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                               <  2:3396432/100‑1 (MQ=255)
                                    cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACGAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                               >  1:2380008/1‑99 (MQ=255)
                                      tGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaa                                                            <  1:3964255/100‑1 (MQ=255)
                                       gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaca                                                           >  2:2514571/1‑98 (MQ=255)
                                       gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                           >  1:15971/1‑100 (MQ=255)
                                       gTTTGCAACGAAAGATCCAGAGCTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCAGGGGCGGTATGAGAAACAAGCTaa                                                            >  2:849528/1‑99 (MQ=255)
                                         ttGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                          <  2:90818/99‑1 (MQ=255)
                                            cAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                      >  2:3378621/1‑100 (MQ=255)
                                             aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                     >  1:3650712/1‑100 (MQ=255)
                                             aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                     <  1:2419027/100‑1 (MQ=255)
                                              aCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGCATGAGAAACAAGCTAAAGGAGAAt                                                    >  2:1435061/1‑100 (MQ=255)
                                                gAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAttt                                                  >  1:2457343/1‑100 (MQ=255)
                                                  aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                >  2:499818/1‑100 (MQ=255)
                                                   cttccgatctttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                               <  1:1129350/90‑1 (MQ=255)
                                                   aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                               <  1:3031671/100‑1 (MQ=255)
                                                   aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                               <  1:1927834/100‑1 (MQ=255)
                                                    tccgatctgTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                              <  1:2565936/92‑1 (MQ=255)
                                                        cAGAGTTCGAAAATAGGATCAACACAGATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                          >  1:2512211/1‑100 (MQ=255)
                                                        cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                          >  1:1380201/1‑100 (MQ=255)
                                                        cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                          <  1:561455/100‑1 (MQ=255)
                                                         agagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCt                                         >  1:2416671/1‑100 (MQ=255)
                                                            gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTTCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                      <  2:722613/100‑1 (MQ=255)
                                                            gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                      >  2:3271259/1‑100 (MQ=255)
                                                            gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                      <  2:1577294/100‑1 (MQ=255)
                                                            gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                      <  2:2151890/100‑1 (MQ=255)
                                                            gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGan                                      <  2:115950/100‑2 (MQ=255)
                                                             ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCagatcggaag                                     >  2:1129343/1‑91 (MQ=255)
                                                              tCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGatat                                    <  2:769476/100‑1 (MQ=255)
                                                                gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                  >  1:2683777/1‑100 (MQ=255)
                                                                  aaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTa                                >  1:3048785/1‑100 (MQ=255)
                                                                     tAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCa                             <  1:836880/100‑1 (MQ=255)
                                                                        gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                          >  2:3508798/1‑100 (MQ=255)
                                                                        gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                          <  2:1543125/100‑1 (MQ=255)
                                                                          tCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGa                        >  2:3195418/1‑100 (MQ=255)
                                                                            aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                      >  2:313615/1‑100 (MQ=255)
                                                                            aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                      >  2:2932807/1‑100 (MQ=255)
                                                                             acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                     <  2:2416657/100‑1 (MQ=255)
                                                                             acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                     >  1:2883643/1‑100 (MQ=255)
                                                                              cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                    >  2:329518/1‑100 (MQ=255)
                                                                              cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                    <  2:2369845/100‑1 (MQ=255)
                                                                               acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                   >  1:2117365/1‑100 (MQ=255)
                                                                                caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                  >  1:4028872/1‑100 (MQ=255)
                                                                                caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                  <  1:175026/100‑1 (MQ=255)
                                                                                 aaaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaa                 <  2:3015600/100‑1 (MQ=255)
                                                                                  tcataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                >  1:249072/8‑100 (MQ=25)
                                                                                  aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCGAAAGGAAAATTTCTCATGTCTGATATGCTAGCAGTCGAACAAGaaa                <  2:562238/100‑1 (MQ=255)
                                                                                   cataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac               >  1:1587780/7‑100 (MQ=25)
                                                                                   aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac               <  2:3769489/100‑1 (MQ=255)
                                                                                    ataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaaca              >  2:3829602/6‑100 (MQ=25)
                                                                                     taccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTAGTATGCTAGCAGTAGGACAAGaaaaaa             >  2:2733983/5‑100 (MQ=25)
                                                                                       ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac           >  1:1493733/3‑100 (MQ=38)
                                                                                       aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac           <  2:3619615/100‑1 (MQ=255)
                                                                                       aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac           <  1:1393148/100‑1 (MQ=255)
                                                                                        aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca          >  1:2323754/1‑100 (MQ=255)
                                                                                        aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca          >  2:623562/1‑100 (MQ=255)
                                                                                        aGCCCCCGGAATGCTGCAAAATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGGATTTCTCATATCTGATATGCTAGCAGTAGAACAATAaacaaaca          >  2:1370733/1‑100 (MQ=255)
                                                                                         gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa         >  1:458121/1‑100 (MQ=38)
                                                                                         gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTACCAGTAGAACAAGAaataaaca          >  2:3982801/1‑99 (MQ=37)
                                                                                         gCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa         >  1:3883326/1‑100 (MQ=255)
                                                                                          cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAaaaaacaa         >  1:1464475/1‑99 (MQ=37)
                                                                                          cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat        >  2:1703203/1‑100 (MQ=38)
                                                                                          cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat        >  1:742203/1‑100 (MQ=38)
                                                                                          cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat        >  1:1048252/1‑100 (MQ=38)
                                                                                            cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga      <  1:3198629/100‑1 (MQ=255)
                                                                                              cGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgatg    <  2:749749/100‑1 (MQ=255)
                                                                                                gAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGtt  >  1:1739783/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
AGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTT  >  NC_000913/2560220‑2560415

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: