Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 464,358 T→G 15.1% intergenic (+50/‑44) ybaV → / → fadM putative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913464,3580TG15.1% 436.2 / 69.8 167intergenic (+50/‑44)ybaV/fadMputative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III
Reads supporting (aligned to +/- strand):  ref base T (64/76);  new base G (10/15);  total (75/92)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-01

Only 100 of 167 total aligned reads displayed.

ATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTC  >  NC_000913/464261‑464455
                                                                                                 |                                                                                                 
aTGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                 >  2:3488437/1‑100 (MQ=255)
   ggCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGa                                                                                              >  2:193676/1‑100 (MQ=255)
   ggCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGa                                                                                              >  1:1310817/1‑100 (MQ=255)
     cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGAcc                                                                                            <  2:796077/100‑1 (MQ=255)
     cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGAcc                                                                                            <  2:2408484/100‑1 (MQ=255)
       aTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctc                                                                                          <  1:3600064/100‑1 (MQ=255)
        ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                         <  1:3102725/100‑1 (MQ=255)
        ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                         <  1:1211387/100‑1 (MQ=255)
        ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                         <  2:3647614/100‑1 (MQ=255)
         tCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctctt                                                                                        <  1:1234577/100‑2 (MQ=39)
          cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGCCTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                       >  1:3527824/1‑100 (MQ=255)
          cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                       >  1:2067266/1‑100 (MQ=255)
          cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                       >  1:2522906/1‑100 (MQ=255)
          cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttc                                                                                       <  2:2964103/100‑3 (MQ=38)
          cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttc                                                                                       <  1:3193550/100‑3 (MQ=38)
           gCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTa                                                                                      <  1:1133156/100‑1 (MQ=255)
           gCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttca                                                                                      <  2:3070652/100‑4 (MQ=38)
           gCTAGTTGAACATAATCTGGCCATATGAACCCATGAATTAATTTCCATACTGGCAATTTTTTGCCAGACTGAAGAAGTCATAACTGCTATGACctcttca                                                                                      <  2:1322071/100‑4 (MQ=255)
            cTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTTCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTAc                                                                                     >  2:3463789/1‑100 (MQ=255)
             tggtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACt                                                                                    <  2:1001244/100‑1 (MQ=255)
              ggtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACtt                                                                                   <  1:1994033/100‑1 (MQ=255)
              ggtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACtt                                                                                   >  1:649878/1‑100 (MQ=255)
                 ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                                <  1:194478/100‑1 (MQ=255)
                 ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                                <  1:88366/100‑1 (MQ=255)
                 ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                                <  1:2686723/100‑1 (MQ=255)
                   aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaac                                                                              <  1:6197/100‑1 (MQ=255)
                   aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaa                                                                               <  1:909861/99‑1 (MQ=255)
                    aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                             >  1:3721636/1‑100 (MQ=255)
                    aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                             <  1:1084223/100‑1 (MQ=255)
                    aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                             <  1:1402059/100‑1 (MQ=255)
                     cGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaagat                                                                            >  1:605978/1‑97 (MQ=255)
                     cGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaa                                                                            <  2:3079735/100‑1 (MQ=255)
                      gTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaac                                                                           <  2:2744732/100‑1 (MQ=255)
                      gTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaac                                                                           <  2:1426099/100‑1 (MQ=255)
                       tAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                          <  1:3411111/100‑1 (MQ=255)
                       tAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                          >  1:3448416/1‑100 (MQ=255)
                       tAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                          >  1:70518/1‑100 (MQ=255)
                        aaTCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaacaa                                                                         <  2:2259249/100‑1 (MQ=255)
                           cTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGt                                                                      >  1:530321/1‑100 (MQ=255)
                             ggCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                    >  2:3538718/1‑100 (MQ=255)
                             ggCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                    >  2:1848857/1‑100 (MQ=255)
                              gCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                   <  1:1848867/100‑1 (MQ=255)
                                  tATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATCCCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTa                                                               <  2:1867324/100‑1 (MQ=255)
                                  tATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTAGCCAGAATGAAGAGGTAATACCAGATATGACCTCTGTACTTATAACAACAACGTacgctt                                                                >  2:1764771/1‑94 (MQ=255)
                                     tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGATCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTg                                                            <  2:2538044/100‑1 (MQ=255)
                                       aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcg                                                          >  2:460353/1‑100 (MQ=255)
                                       aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcg                                                          <  2:2053694/100‑1 (MQ=255)
                                           tGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                       <  1:80005/99‑1 (MQ=255)
                                           tGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                      <  1:3679383/100‑1 (MQ=255)
                                           tGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                      <  2:3136560/100‑1 (MQ=255)
                                             taattaatTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGc                                                    >  1:2238940/1‑100 (MQ=255)
                                              aattaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                   <  1:395033/100‑1 (MQ=255)
                                              aattaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                   <  1:1863975/100‑1 (MQ=255)
                                               attaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                   >  1:3186893/1‑99 (MQ=255)
                                                ttaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                 >  1:3797978/1‑100 (MQ=255)
                                                 taattTGCATAGTGGCAATTTTTTGCCAGACTGAACAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                <  2:784249/100‑1 (MQ=255)
                                                    tttGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                             >  2:2662168/1‑100 (MQ=255)
                                                    tttGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                             >  2:675077/1‑100 (MQ=255)
                                                      tGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaaa                                           >  2:1114287/1‑100 (MQ=255)
                                                       tccgatctgCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAAt                                          <  2:1766476/92‑1 (MQ=255)
                                                        cATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATc                                         >  1:2804607/1‑100 (MQ=255)
                                                           aGTGGCAATGTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaaa                                      <  2:592218/100‑1 (MQ=255)
                                                            gTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAg                                     <  1:2206308/100‑1 (MQ=255)
                                                             tggcaaaaaattGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGt                                    >  2:2201346/11‑100 (MQ=255)
                                                              ggCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGtt                                   >  2:1263847/1‑100 (MQ=255)
                                                               gCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAAtagatcgga                                  >  1:1766486/1‑92 (MQ=255)
                                                               gCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTc                                  >  2:2624402/1‑100 (MQ=255)
                                                                   ttttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTgg                              <  2:3507251/100‑1 (MQ=255)
                                                                     ttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGat                            <  2:307400/100‑1 (MQ=255)
                                                                         gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTATTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatca                        <  2:3205698/100‑1 (MQ=255)
                                                                         gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCCAAGTTCGTGGATatca                        >  2:3353245/1‑100 (MQ=255)
                                                                         gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatca                        <  2:3425003/100‑1 (MQ=255)
                                                                          ccAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatcat                       >  1:2462925/1‑100 (MQ=255)
                                                                             gACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtctc                    >  1:1196547/1‑100 (MQ=255)
                                                                              aCTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCg                   >  2:1622196/1‑100 (MQ=255)
                                                                               gtacagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGa                  >  1:2901076/5‑100 (MQ=38)
                                                                               gtacagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGa                  >  2:3641837/5‑100 (MQ=38)
                                                                                 gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg                <  1:1959207/100‑1 (MQ=255)
                                                                                 gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg                >  2:146358/1‑100 (MQ=255)
                                                                                 gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg                >  2:774391/1‑100 (MQ=255)
                                                                                 acagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGGTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg                >  1:1281328/3‑100 (MQ=38)
                                                                                  cagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGt               >  1:3485520/2‑100 (MQ=39)
                                                                                   agagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGtt              >  2:3560450/1‑100 (MQ=255)
                                                                                     agGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTa            <  1:2530529/100‑1 (MQ=255)
                                                                                      ggTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAc           >  1:1077953/1‑100 (MQ=255)
                                                                                      ggTCATACCAGTTATGACCTCTGTACGTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAc           >  2:754129/1‑100 (MQ=255)
                                                                                       gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAcc          <  1:2893614/100‑1 (MQ=255)
                                                                                       gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAcc          <  1:2757234/100‑1 (MQ=255)
                                                                                        tCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa         <  2:2716050/100‑1 (MQ=255)
                                                                                        tCATACCAGTTATGACCTCTGTACTTAGAACAACAACGTAAGGTTATTGCTCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa         <  2:2824457/100‑1 (MQ=255)
                                                                                         cATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAg        >  1:3324218/1‑100 (MQ=255)
                                                                                         cATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAg        >  2:3605996/1‑100 (MQ=255)
                                                                                           tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa      <  1:3423335/100‑1 (MQ=255)
                                                                                           tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa      <  2:2956608/100‑1 (MQ=255)
                                                                                           tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa      >  2:3739118/1‑100 (MQ=255)
                                                                                           tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa      <  1:2624421/100‑1 (MQ=255)
                                                                                             ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACg    <  1:1482720/100‑1 (MQ=255)
                                                                                             ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACg    >  1:3820075/1‑100 (MQ=255)
                                                                                              cAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGt   >  1:749585/1‑100 (MQ=255)
                                                                                               aGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTc  <  1:155229/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
ATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTC  >  NC_000913/464261‑464455

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