Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,200,184 (A)7→5 intergenic (+90/‑96) yjaG → / → hupA DUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,200,1790A.98.6% 316.8 / ‑2.2 73intergenic (+85/‑102)yjaG/hupADUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (28/44);  total (28/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
*NC_0009134,200,1800A.98.6% 316.9 / ‑2.4 73intergenic (+86/‑101)yjaG/hupADUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (28/44);  total (28/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAG  >  NC_000913/4200081‑4200272
                                                                                                  ||                                                                                            
tttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaa                                                                                              <  1:780937/100‑1 (MQ=255)
  tttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaa                                                                                            >  2:828452/1‑100 (MQ=255)
    tCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaag                                                                                          <  1:2288892/100‑2 (MQ=255)
     cAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaag                                                                                          >  2:1101606/1‑98 (MQ=255)
      aGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaagtg                                                                                        <  1:773311/100‑4 (MQ=255)
      aGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaagtg                                                                                        >  2:2019844/1‑97 (MQ=255)
      aGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaagtg                                                                                        <  1:1590138/100‑4 (MQ=255)
      aGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaa                                                                                            >  2:2345244/1‑96 (MQ=255)
       gCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaagtgg                                                                                       <  2:2373494/100‑5 (MQ=255)
       gCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaagtgg                                                                                       >  1:976219/1‑96 (MQ=255)
        cAATAAGACCAGAAAACGTTATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaagtggc                                                                                      <  1:2020632/100‑6 (MQ=255)
        cAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaagtggc                                                                                      >  1:1898822/1‑95 (MQ=255)
        cAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaagtggc                                                                                      <  2:1464879/100‑6 (MQ=255)
         aaTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaagtggct                                                                                     >  2:234565/1‑94 (MQ=255)
          aTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaagtggcta                                                                                    >  1:2856950/1‑93 (MQ=255)
            aaGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATc                                                                                <  1:1773595/100‑1 (MQ=255)
            aaGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATc                                                                                >  2:825983/1‑100 (MQ=255)
              gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCgg                                                                              >  2:207182/1‑100 (MQ=255)
                  ngAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCg                                                                          >  1:291515/2‑100 (MQ=255)
                  aGAAAACGTGCTTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCg                                                                          <  2:1898808/100‑1 (MQ=255)
                  aGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCg                                                                          >  1:1376393/1‑100 (MQ=255)
                    aaaaCGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCgtg                                                                        <  1:940683/100‑1 (MQ=255)
                    aaaaCGTGATTTAACGCCTGATGTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCgtg                                                                        <  2:670483/100‑1 (MQ=255)
                       aCGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTAt                                                                     <  1:1605771/100‑1 (MQ=255)
                        cGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATg                                                                    <  2:1049718/100‑1 (MQ=255)
                         gTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGc                                                                   <  2:126523/100‑1 (MQ=255)
                         gTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGc                                                                   <  2:580540/100‑1 (MQ=255)
                         gTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGc                                                                   <  2:264146/100‑1 (MQ=255)
                          tGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCa                                                                  <  1:250876/100‑1 (MQ=255)
                          tGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTGCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑ACAAAGTGGCTATCGGTGCGTGTATGCa                                                                  <  1:1758975/100‑1 (MQ=255)
                           gATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAg                                                                 >  1:56535/1‑100 (MQ=255)
                           gATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAg                                                                 <  1:710300/100‑1 (MQ=255)
                            aTTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAgg                                                                <  1:571826/100‑1 (MQ=255)
                            aTTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAgg                                                                <  2:508498/100‑1 (MQ=255)
                             tttAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGga                                                               >  2:786014/1‑100 (MQ=255)
                              ttAACGCCTGAGTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGgag                                                              <  2:698188/100‑1 (MQ=255)
                                  cGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTg                                                          >  1:211592/1‑100 (MQ=255)
                                  cGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTg                                                          >  2:516564/1‑100 (MQ=255)
                                   gCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGc                                                         <  2:491367/100‑1 (MQ=255)
                                   gCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGc                                                         <  1:1299169/100‑1 (MQ=255)
                                      tGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTAt                                                      <  2:976217/100‑1 (MQ=255)
                                      tGATTTGTCGTACCTGGAGTCTTCCCTTTCACCGCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGAAGAGTGCTAt                                                      >  2:1896664/1‑100 (MQ=255)
                                        aTTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTc                                                    >  1:284469/1‑100 (MQ=255)
                                         tttGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGTGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTCTTc                                                    >  2:2377482/1‑99 (MQ=255)
                                           tGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTgg                                                 >  1:2389191/1‑100 (MQ=255)
                                           tGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTgg                                                 >  2:1367260/1‑100 (MQ=255)
                                              cGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAt                                              <  2:2335321/100‑1 (MQ=255)
                                              cGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAt                                              <  1:2531697/100‑1 (MQ=255)
                                               gTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAtt                                             <  1:234566/100‑1 (MQ=255)
                                                tACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAt                                              <  2:1934063/98‑1 (MQ=255)
                                                         tCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCAc                                   <  1:1081276/100‑1 (MQ=255)
                                                          cTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACt                                  >  1:2193734/1‑100 (MQ=255)
                                                            tCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCg                                <  1:2097385/100‑1 (MQ=255)
                                                             cccTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTcaa                               >  2:1759428/1‑98 (MQ=255)
                                                             cccTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGa                               >  2:2711193/1‑100 (MQ=255)
                                                              ccTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAt                              >  2:1732809/1‑100 (MQ=255)
                                                              ccTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAt                              <  1:721025/100‑1 (MQ=255)
                                                               cTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATg                             >  2:1507444/1‑100 (MQ=255)
                                                                 ttCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCt                           >  2:221044/1‑100 (MQ=255)
                                                                  tCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCtt                          >  1:1375788/1‑100 (MQ=255)
                                                                   ccccccccGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGGGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTa                         <  1:2546492/98‑1 (MQ=255)
                                                                   cGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTa                         <  1:1507452/100‑1 (MQ=255)
                                                                    gCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAg                        >  1:147421/1‑100 (MQ=255)
                                                                    gCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAg                        <  1:863573/100‑1 (MQ=255)
                                                                    gCCCCCCGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAg                        <  1:12684/100‑1 (MQ=255)
                                                                     ccccccGTCTGGTCTACACTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGc                       <  2:831924/100‑1 (MQ=255)
                                                                      cccccGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCa                      <  2:2594014/100‑1 (MQ=255)
                                                                      cccccGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCa                      <  1:2884342/100‑1 (MQ=255)
                                                                          cGTCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCg                  <  1:185058/100‑1 (MQ=255)
                                                                           gTCTGGTCTACATTTGGGGGGCTAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGc                   <  2:696642/100‑1 (MQ=255)
                                                                            tCTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGAt                <  1:825984/100‑1 (MQ=255)
                                                                             cTGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATa               >  1:1130059/1‑100 (MQ=255)
                                                                              tGGTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATaa              >  2:2566849/1‑100 (MQ=255)
                                                                               tgTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATaaa             <  2:2350943/99‑1 (MQ=255)
                                                                               ggTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATaaa             >  2:1031640/1‑100 (MQ=255)
                                                                                tagaccagacggggggCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAac            <  2:2364759/90‑1 (MQ=255)
                                                                                gTCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATaaaa            >  1:2350958/1‑99 (MQ=255)
                                                                                 tCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAaca           <  2:284468/100‑1 (MQ=255)
                                                                                 tCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAaca           >  2:2144027/1‑100 (MQ=255)
                                                                                 tCTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAaca           >  1:2283338/1‑100 (MQ=255)
                                                                                  cTACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAacac          <  1:2103012/100‑1 (MQ=255)
                                                                                   tACATTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAAcaca         <  2:2594128/100‑1 (MQ=255)
                                                                                      aTTTGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTg      <  1:1504200/100‑1 (MQ=255)
                                                                                        ttGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAAGACATTg      <  2:2792354/98‑1 (MQ=255)
                                                                                         tGGGGGGCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTaa   >  2:2310940/1‑100 (MQ=255)
                                                                                          ggggggCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAg  <  2:1130055/100‑1 (MQ=255)
                                                                                          ggggggCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAg  <  1:2254201/100‑1 (MQ=255)
                                                                                          ggggggCG‑‑AAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAg  <  1:1305096/100‑1 (MQ=255)
                                                                                                  ||                                                                                            
TTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAG  >  NC_000913/4200081‑4200272

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