Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 204,870 T→G D126E (GAT→GAG rnhB → ribonuclease HII, degrades RNA of DNA‑RNA hybrids

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913204,8700TG98.4% 225.3 / NA 64D126E (GAT→GAGrnhBribonuclease HII, degrades RNA of DNA‑RNA hybrids
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (39/24);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGACGAGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGATAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATCTGGAAAAACTG  >  NC_000913/204724‑205014
                                                                                                                                                  |                                                                                                                                                
aTCGACGAGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGcc                                                                                                                                              >  2:781952/1‑151 (MQ=255)
  cgacgaGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCcgc                                                                                                                                            >  2:527584/1‑151 (MQ=255)
     cgaGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGt                                                                                                                                          >  1:379350/1‑150 (MQ=255)
     cgaGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTa                                                                                                                                         <  1:609251/151‑1 (MQ=255)
      gaGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTAc                                                                                                                                        >  1:553368/1‑151 (MQ=255)
       aGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTAcc                                                                                                                                       >  2:758656/1‑151 (MQ=255)
       aGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTAcc                                                                                                                                       >  2:147895/1‑151 (MQ=255)
        gCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTAcc                                                                                                                                       <  2:233616/150‑1 (MQ=255)
          tGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCgg                                                                                                                                     >  2:636249/1‑150 (MQ=255)
             aCATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAAt                                                                                                                                 <  1:913665/151‑1 (MQ=255)
               aTTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCa                                                                                                                               <  1:11795/151‑1 (MQ=255)
                    tcATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGcc                                                                                                                          >  2:462850/1‑151 (MQ=255)
                     cATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCcg                                                                                                                         <  1:132007/151‑1 (MQ=255)
                        gCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCcgcg                                                                                                                       <  2:39703/150‑1 (MQ=255)
                                gCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTg                                                                                                              >  1:404790/1‑151 (MQ=255)
                                  tGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGc                                                                                                            >  1:319465/1‑151 (MQ=255)
                                   ggCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCg                                                                                                           <  2:248767/151‑1 (MQ=255)
                                      gATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGaa                                                                                                         >  1:40876/1‑150 (MQ=255)
                                                  cGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAgtgacgcgtg                                                                                              >  1:589254/1‑149 (MQ=255)
                                                     cGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGcc                                                                                         <  1:869731/151‑1 (MQ=255)
                                                      gCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCg                                                                                        >  1:358183/1‑151 (MQ=255)
                                                              gCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATggcg                                                                                <  2:319762/151‑1 (MQ=255)
                                                              gCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATgg                                                                                  <  1:101303/149‑1 (MQ=255)
                                                                atatTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATggcggc                                                                              <  1:855689/151‑1 (MQ=255)
                                                                     gcgcCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTg                                                                          <  1:742766/150‑1 (MQ=255)
                                                                      cgcCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTgg                                                                         >  1:304855/1‑150 (MQ=255)
                                                                        ccGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGata                                                                      <  1:383891/151‑1 (MQ=255)
                                                                          ggAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATAtt                                                                    >  1:791304/1‑151 (MQ=255)
                                                                          ggAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATAtt                                                                    >  2:62183/1‑151 (MQ=255)
                                                                           gAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTg                                                                   >  1:801196/1‑151 (MQ=255)
                                                                             atatGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGtt                                                                 >  2:835876/1‑151 (MQ=255)
                                                                             atatGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGtt                                                                 <  2:219467/151‑1 (MQ=255)
                                                                             atatGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGtt                                                                 >  1:242590/1‑151 (MQ=255)
                                                                                tgtgTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCCGGATATTGTTTTc                                                              >  1:2283/1‑151 (MQ=255)
                                                                                 gtgtTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTcc                                                             >  2:860491/1‑151 (MQ=255)
                                                                                  tgtTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTccc                                                            >  2:92193/1‑151 (MQ=255)
                                                                                   gtTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTccc                                                            <  1:304566/150‑1 (MQ=255)
                                                                                     tgattgatGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGc                                                          >  1:752012/1‑150 (MQ=255)
                                                                                     tgattgatGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGc                                                          <  1:579277/150‑1 (MQ=255)
                                                                                            tGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATgg                                                   >  2:35163/1‑150 (MQ=255)
                                                                                            tGGAAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAatat                                                     >  1:715829/1‑148 (MQ=255)
                                                                                             ggTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGtt                                                 >  2:404654/1‑151 (MQ=255)
                                                                                               tAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGtttt                                               >  1:137921/1‑151 (MQ=255)
                                                                                                aaCCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGtttt                                               <  1:859415/150‑1 (MQ=255)
                                                                                                 aCCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGc                                             <  1:669094/151‑1 (MQ=255)
                                                                                                     cTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCa                                          >  2:83277/1‑150 (MQ=255)
                                                                                                       gCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAac                                        <  2:759814/150‑1 (MQ=255)
                                                                                                       gCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAac                                        >  1:60123/1‑150 (MQ=255)
                                                                                                         ccGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAAcaca                                     <  2:819279/151‑1 (MQ=255)
                                                                                                            aaaTTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAg                                  <  2:80296/151‑1 (MQ=255)
                                                                                                                 aCCGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTTGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTa                              <  1:607592/150‑1 (MQ=255)
                                                                                                                   cGATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTAccc                           >  2:191666/1‑151 (MQ=255)
                                                                                                                    gATGCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTAccc                           >  2:451647/1‑150 (MQ=255)
                                                                                                                       gCCTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAAcc                       <  1:198289/151‑1 (MQ=255)
                                                                                                                        ccTGCGATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCg                      >  1:396631/1‑151 (MQ=255)
                                                                                                                             gATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCtttt                 >  1:638086/1‑151 (MQ=255)
                                                                                                                             gATGGCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCtttt                 >  2:252673/1‑151 (MQ=255)
                                                                                                                                ggCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCAt              >  2:276476/1‑151 (MQ=255)
                                                                                                                                 gCTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATc             >  2:183949/1‑151 (MQ=255)
                                                                                                                                  cTGTGGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATCt            >  1:352251/1‑151 (MQ=255)
                                                                                                                                   tgtgGTGAAAGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATCTg           >  2:162189/1‑151 (MQ=255)
                                                                                                                                          aaaGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATCTGGAAAAAc    >  1:701593/1‑151 (MQ=255)
                                                                                                                                            aGGCGAGAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATCTGGAAAAACTg  <  2:490129/151‑1 (MQ=255)
                                                                                                                                             ggCGATAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATCTGGAAAAACTg  <  1:562884/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATCGACGAGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGATAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATCTGGAAAAACTG  >  NC_000913/204724‑205014

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: