Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,978,570 (G)6→7 coding (71/1386 nt) yifK → putative amino acid/amine transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,978,5641.G98.3% 555.7 / 1.8 176coding (65/1386 nt)yifKputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (86/87);  total (88/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 213 total aligned reads displayed.

AAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCAT  >  NC_000913/3978472‑3978651
                                                                                              |                                                                                        
aaaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                                                           >  1:3828142/1‑93 (MQ=255)
aaaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                                                           <  2:5577821/93‑1 (MQ=255)
aaaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                                                           <  1:5625347/93‑1 (MQ=255)
aaaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                                                           >  2:5530954/1‑93 (MQ=255)
aaaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                                                           >  1:1865377/1‑93 (MQ=255)
 aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                                                                         <  1:105423/93‑1 (MQ=255)
 aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                                                                         >  1:2465846/1‑93 (MQ=255)
 aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                                                                         <  2:548406/93‑1 (MQ=255)
 aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                                                                         >  1:4457268/1‑93 (MQ=255)
 aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                                                                         <  2:3845078/93‑1 (MQ=255)
  aGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgg                                                                                        <  1:3298574/93‑1 (MQ=255)
    cATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                                                                     <  2:6046870/93‑1 (MQ=255)
     aTGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggg                                                                                    <  2:5907025/93‑1 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                                                                   >  1:1723187/1‑93 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                                                                   >  2:3882213/1‑93 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                                                                   >  1:3082197/1‑93 (MQ=255)
        cAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggc                                                                                 <  1:221584/93‑3 (MQ=255)
        cAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggc                                                                                 >  1:4601725/1‑91 (MQ=255)
            cacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                                                              <  1:4857324/93‑1 (MQ=255)
            cacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                                                              >  2:2982818/1‑93 (MQ=255)
            cacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                                                              <  2:6007456/93‑1 (MQ=255)
            cacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                                                              >  1:2295037/1‑93 (MQ=255)
             acaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                                                             <  1:2707448/93‑1 (MQ=255)
             acaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                                                             <  1:1903208/93‑1 (MQ=255)
               aacacGAGGATTTCTGGCCGATAACAAACCAGGCACAGAGCGTGGGCTGGAAGCGCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                                                           <  2:3835814/93‑1 (MQ=255)
               aacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                                                           >  2:3761670/1‑93 (MQ=255)
               aacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                                                           <  2:1970249/93‑1 (MQ=255)
                acacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                                                                          <  1:3352847/93‑1 (MQ=255)
                acacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                                                                          >  2:4646019/1‑93 (MQ=255)
                 cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                                                                         <  2:2645241/93‑1 (MQ=255)
                   cGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                                                                       >  2:120142/1‑93 (MQ=255)
                      ggATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                                                                    >  2:4549592/1‑93 (MQ=255)
                         tttATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                                                                 >  2:3827294/1‑93 (MQ=255)
                           tATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                                                               <  1:3191983/93‑1 (MQ=255)
                             tGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                                                             >  1:5730641/1‑93 (MQ=255)
                              ggCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            >  1:5079684/1‑93 (MQ=255)
                               gCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                                                           >  1:4561300/1‑93 (MQ=255)
                                cAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg                                                          <  2:186756/93‑1 (MQ=255)
                                 aGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg                                                         <  1:1660422/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGGCTGTTTATggg                                                        >  2:4036561/1‑93 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                        <  1:1503311/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                        <  2:5691859/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                        >  2:6225572/1‑93 (MQ=255)
                                   aTAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                       >  1:5903132/1‑93 (MQ=255)
                                   aTAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                       <  1:2503193/93‑1 (MQ=255)
                                    tAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg                                                      >  2:2172821/1‑93 (MQ=255)
                                     aaCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                     <  1:38796/93‑1 (MQ=255)
                                     aaCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                     <  1:3720439/93‑1 (MQ=255)
                                      aCAAACCAGAGCTACGGCGTGGGCTGGAAGCTCGACATATCGAACTCACCGCCCGGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc                                                    <  2:1160063/93‑1 (MQ=255)
                                      aCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc                                                    >  1:1695123/1‑93 (MQ=255)
                                        aacccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCCCCATTGGCGTCGGCCTGTTTATGGGGgccgc                                                  <  1:3596372/90‑1 (MQ=255)
                                         aaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc                                                 >  2:4417988/1‑93 (MQ=255)
                                          aCCAGAGCTACAGCGTGGTCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                <  1:813798/93‑1 (MQ=255)
                                          aCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                <  2:4335899/93‑1 (MQ=255)
                                           ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                               <  2:4823677/93‑1 (MQ=255)
                                           ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                               <  1:5929192/93‑1 (MQ=255)
                                           ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                               >  1:2772439/1‑93 (MQ=255)
                                            cAGAGCTAAAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt                                              <  1:1180968/93‑1 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc                                          <  1:6022694/93‑1 (MQ=255)
                                                 cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGTGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCt                                         >  1:2404874/1‑93 (MQ=255)
                                                   aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGa                                       <  1:1261706/93‑1 (MQ=255)
                                                     aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGaaa                                     <  1:6154446/93‑1 (MQ=255)
                                                     aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGaaa                                     >  2:2095949/1‑93 (MQ=255)
                                                     aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGaaa                                     >  1:4409141/1‑93 (MQ=255)
                                                      gCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                    >  2:1031397/1‑93 (MQ=255)
                                                       cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATg                                   <  2:2355832/93‑1 (MQ=255)
                                                       cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATg                                   <  1:2677017/93‑1 (MQ=255)
                                                        gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgg                                  <  1:4533203/93‑1 (MQ=255)
                                                         tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATggg                                 >  1:693319/1‑93 (MQ=255)
                                                          gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggc                                >  2:504979/1‑93 (MQ=255)
                                                          gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggc                                >  2:4490806/1‑93 (MQ=255)
                                                           ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggcc                               <  2:1822458/93‑1 (MQ=255)
                                                                gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggcc                          <  2:4867284/93‑1 (MQ=255)
                                                                 aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCa                         >  2:672718/1‑93 (MQ=255)
                                                                   gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATc                       <  1:801869/93‑1 (MQ=255)
                                                                    cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATcc                      <  1:5523/93‑1 (MQ=255)
                                                                    cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATcc                      >  2:1521501/1‑93 (MQ=255)
                                                                       gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGtt                   >  2:5248892/1‑92 (MQ=255)
                                                                       gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTa                   >  1:3990283/1‑93 (MQ=255)
                                                                        aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAt                  >  1:5738198/1‑93 (MQ=255)
                                                                        aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAt                  >  1:4797052/1‑93 (MQ=255)
                                                                        aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAt                  <  1:3969989/93‑1 (MQ=255)
                                                                         cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAtt                 <  1:53593/93‑1 (MQ=255)
                                                                          atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttg                <  2:4973564/93‑1 (MQ=255)
                                                                            atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt              >  1:3975394/1‑93 (MQ=255)
                                                                             tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgttg             <  1:2575869/93‑1 (MQ=255)
                                                                              cGAACTCATCGCCCTGGGGGGGCGCCTTTGGCGTCGCCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTgg            <  1:2648421/93‑1 (MQ=255)
                                                                              cGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTgg            <  1:5588440/93‑1 (MQ=255)
                                                                              cGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTgg            <  2:2914431/93‑1 (MQ=255)
                                                                                aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc          >  2:2736345/1‑93 (MQ=255)
                                                                                aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc          >  2:5298928/1‑93 (MQ=255)
                                                                                   tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTCTTGTTGGCCtat       >  1:4714973/1‑93 (MQ=255)
                                                                                    cATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtata      <  2:3606301/93‑1 (MQ=255)
                                                                                    cATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtata      <  2:888456/93‑1 (MQ=255)
                                                                                     aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTCTTGTTGGCCTatat     >  2:1256916/1‑93 (MQ=255)
                                                                                      tCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatc    >  1:2306028/1‑93 (MQ=255)
                                                                                       cGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatca   <  2:5424821/93‑1 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCCTCCGTAGTGTTGGCCTATatcat  >  2:2079624/1‑93 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat  <  1:5812757/93‑1 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat  <  2:859243/93‑1 (MQ=255)
                                                                                              |                                                                                        
AAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCAT  >  NC_000913/3978472‑3978651

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