Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,393 A→G 5.1% F230L (TTC→CTC)  yggE ← oxidative stress defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,067,3930AG 5.1% 245.1 / 3.2 78F230L (TTC→CTC) yggEoxidative stress defense protein
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base G (2/2);  total (38/40)
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.94e-01

GGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCC  >  NC_000913/3067247‑3067540
                                                                                                                                                  |                                                                                                                                                   
ggTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGATGCAGGTGTTTTAGCGGGTTGTTGATCAAAAGGTTCTAACTGGaa                                                                                                                                                     >  8:309369/1‑147 (MQ=255)
      tACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcaca                                                                                                                                              >  1:375994/1‑148 (MQ=255)
            aTGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAcc                                                                                                                                        <  7:298902/148‑1 (MQ=255)
               gTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTg                                                                                                                                      <  2:328089/147‑1 (MQ=255)
               gTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGa                                                                                                                                     >  7:198162/1‑148 (MQ=255)
               gTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGa                                                                                                                                     >  7:226180/1‑148 (MQ=255)
                       gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAGGACCACATCGACCTGATCATCaaa                                                                                                                             >  4:275505/1‑148 (MQ=255)
                       gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaaa                                                                                                                             <  5:93330/148‑1 (MQ=255)
                           attacagaaaaAGGAGGCAATATCGGGTAAAGGCAGTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGa                                                                                                                         <  3:273729/148‑1 (MQ=255)
                                agaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGc                                                                                                                    <  8:198162/148‑1 (MQ=255)
                                     aaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGa                                                                                                                                                   >  4:111039/1‑112 (MQ=255)
                                     aaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGa                                                                                                                                                   <  3:111039/112‑1 (MQ=255)
                                        gaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAGGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCt                                                                                                            >  3:330165/1‑148 (MQ=255)
                                             aaTATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaa                                                                                                                              >  7:237388/1‑125 (MQ=255)
                                             aaTATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaa                                                                                                                              <  8:237388/125‑1 (MQ=255)
                                              atatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAg                                                                                                      >  5:136692/1‑148 (MQ=255)
                                               tatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGc                                                                                                                    >  7:97908/1‑133 (MQ=255)
                                               tatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGc                                                                                                                    <  8:97908/133‑1 (MQ=255)
                                                 tCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGtt                                                                                                    <  2:170983/147‑1 (MQ=255)
                                                     gTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAGGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTg                                                                                               <  3:275505/148‑1 (MQ=255)
                                                     gTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTg                                                                                               >  3:368114/1‑148 (MQ=255)
                                                         agtcATTAGCCCGACGAATACGTCGGGCTAAAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACATGATCATCAAACTGAATAGCGGCCTGCTCGTCAGTTTCCTCGGCg                                                                                           <  6:341758/145‑1 (MQ=255)
                                                          ggCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCgg                                                                                          <  4:368114/148‑1 (MQ=255)
                                                            cATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGa                                                                                                                                          <  2:328796/98‑1 (MQ=255)
                                                            cATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGa                                                                                                                                          >  1:328796/1‑98 (MQ=255)
                                                             cttatgccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGaca                                                                                       <  2:216650/142‑1 (MQ=255)
                                                             aTTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGaca                                                                                       <  3:335185/148‑1 (MQ=255)
                                                               tAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACAcc                                                                                     >  6:359222/1‑148 (MQ=255)
                                                               tAGCCCGACGAATACGTCGGGCTACAAATATTAGTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACAcc                                                                                     >  3:220035/1‑148 (MQ=255)
                                                                  cccGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgc                                                                                  >  8:31021/1‑148 (MQ=255)
                                                                     gacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTTGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGacac                                                                                      <  5:312027/141‑1 (MQ=255)
                                                                     gacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGcgcg                                                                               >  1:366002/1‑148 (MQ=255)
                                                                      acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGa                                                                                                                         <  6:299794/105‑1 (MQ=255)
                                                                      acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGa                                                                                                                         >  5:299794/1‑105 (MQ=255)
                                                                      acgaAAACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGCGCTGAGGATCAGACTGACTAGCGGCCGGCTCGTAAGTTGCCTGGGGGGACACGgcgctgg                                                                             >  6:312027/1‑145 (MQ=255)
                                                                         aaTACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCAt                                                                           >  3:353197/1‑148 (MQ=255)
                                                                               tCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCtt                                                                     <  2:375994/148‑1 (MQ=255)
                                                                               tCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCtt                                                                     <  3:290189/148‑1 (MQ=255)
                                                                               tCGGGCTACAAATATTATTGTGCGTCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACTCCGGCGCGGCATCGGCtt                                                                     <  3:104890/148‑1 (MQ=255)
                                                                                 gggCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCttt                                                                    <  2:400905/147‑1 (MQ=255)
                                                                                 gggCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCttt                                                                    <  7:31021/147‑1 (MQ=255)
                                                                                  ggCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTtca                                                                  <  6:89514/148‑1 (MQ=255)
                                                                                       cAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTcatcat                                                              <  4:220035/147‑1 (MQ=255)
                                                                                        aaaTATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATcc                                                            >  7:195443/1‑148 (MQ=255)
                                                                                         aaTATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCg                                                           >  1:176435/1‑148 (MQ=255)
                                                                                           tattatTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAAc                                                                                                                            >  7:35229/1‑81 (MQ=255)
                                                                                           tattatTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAAc                                                                                                                            <  8:35229/81‑1 (MQ=255)
                                                                                                   ttctgcAGGTGTTTTAGCGGGTTGTTGATCCACAGGTGCTAACTGGAAGACCACATCGACCTTATCATCAAACTGAATAGCGGCCTGCTCTTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtggg                                                 <  4:54871/146‑1 (MQ=255)
                                                                                                     ctgcAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATtgggct                                               >  6:195562/1‑148 (MQ=255)
                                                                                                           gtgtTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTgggctgggct                                          >  8:85217/1‑147 (MQ=255)
                                                                                                             gtTTTAGCGGGTTGTTTATACACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCTGCCTTCTAGTAAGTTTCCTGTGCGGACACCGGCGCGACATCGGCATTCATCATCCGCACCATAGGGCTGGGCTGa                                        >  4:56933/1‑147 (MQ=255)
                                                                                                              ttttAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGAt                                       <  7:309369/147‑1 (MQ=255)
                                                                                                               tttAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAg                                     >  1:56777/1‑148 (MQ=255)
                                                                                                               tatAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAg                                     >  2:9108/3‑148 (MQ=255)
                                                                                                                ttAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAg                                     <  8:226180/147‑1 (MQ=255)
                                                                                                                  aGCGGGTTGTTGATCCACAGGTTCTAACTGGAGGACCACATCGACCTGATCATCAAACTGAATAGCTGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCGCCACCATTGGGCTGGGCTGATAGtt                                   <  4:330165/147‑1 (MQ=255)
                                                                                                                   gCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTg                                  <  1:9108/147‑1 (MQ=255)
                                                                                                                     gggTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGc                                                                             >  4:184985/1‑102 (MQ=255)
                                                                                                                     gggTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGc                                                                             <  3:184985/102‑1 (MQ=255)
                                                                                                                           ttgATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAAc                             >  4:340868/1‑144 (MQ=255)
                                                                                                                           ttgATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAAc                             <  3:340868/144‑1 (MQ=255)
                                                                                                                              aTCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAg                      >  8:61014/1‑148 (MQ=255)
                                                                                                                               tCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAgc                     >  8:110149/1‑148 (MQ=255)
                                                                                                                                 cacaGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAgcgc                   <  6:136692/148‑1 (MQ=255)
                                                                                                                                 cacaGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAgcgc                   <  1:309476/148‑1 (MQ=255)
                                                                                                                                 cacaGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAgcgc                   <  2:366002/148‑1 (MQ=255)
                                                                                                                                 cacaGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAgcg                    >  2:309476/1‑147 (MQ=255)
                                                                                                                                  acaGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTTGGCTGGGCTGATAGTTGGAAACATGGTAGCGCa                  >  6:57980/1‑148 (MQ=255)
                                                                                                                                  acaGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGAt                                       <  8:398685/127‑1 (MQ=255)
                                                                                                                                  acaGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGAt                                       >  7:398685/1‑127 (MQ=255)
                                                                                                                                       ttCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCta             >  6:322495/1‑148 (MQ=255)
                                                                                                                                         cTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCtat            <  2:56777/147‑1 (MQ=255)
                                                                                                                                           aaCTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTa                                                                                                        <  3:109184/53‑1 (MQ=255)
                                                                                                                                           aaCTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTa                                                                                                        >  4:109184/1‑53 (MQ=255)
                                                                                                                                            aCTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATAc        <  4:353197/148‑1 (MQ=255)
                                                                                                                                             cTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTcc                                                                                                 <  3:157102/58‑1 (MQ=255)
                                                                                                                                             cTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTcc                                                                                                 >  4:157102/1‑58 (MQ=255)
                                                                                                                                             cTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCtatata         >  3:68926/1‑146 (MQ=255)
                                                                                                                                               ggAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCgg     >  4:214192/1‑148 (MQ=255)
                                                                                                                                                  aGACCACATCGACCTAATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGccc  >  7:373502/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCC  >  NC_000913/3067247‑3067540

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