Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,498,172 A→C E439D (GAA→GAC)  opgD → OPG biosynthetic periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,498,1720AC98.4% 227.8 / ‑2.9 64E439D (GAA→GAC) opgDOPG biosynthetic periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (29/34);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

CTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGAAGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGGAAGCGAAGCAAATCGAAATT  >  NC_000913/1498028‑1498319
                                                                                                                                                |                                                                                                                                                   
cTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTgg                                                                                                                                               >  2:339848/1‑151 (MQ=255)
   cTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGgcg                                                                                                                                            >  1:482671/1‑151 (MQ=255)
   cTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGgcg                                                                                                                                            >  2:736828/1‑151 (MQ=255)
   cTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGgc                                                                                                                                             <  1:724892/150‑1 (MQ=255)
   cTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGAAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTggg                                                                                                                                              <  1:294/149‑1 (MQ=255)
    ttCTGGCAGCCAGAAAAAGTTGTAAAAGCAGGTGATGATTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATCGCCCATTAGCGCGCGTTATGGCGACTCGTACCGGCATGGGCGGTTTCTCGGACGGt                                                                                                                                                  <  2:471346/144‑1 (MQ=255)
    ttCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGGGGGCGAAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGgc                                                                                                                                             <  2:325585/149‑1 (MQ=255)
    ttCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGgc                                                                                                                                             <  1:293522/149‑1 (MQ=255)
     tCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGcc                                                                                                                                          >  2:555135/1‑151 (MQ=255)
        ggCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAg                                                                                                                                        <  2:212289/150‑1 (MQ=255)
        ggCAGCAAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAg                                                                                                                                        <  1:262732/150‑1 (MQ=255)
           aGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTg                                                                                                                                     >  2:282964/1‑150 (MQ=255)
           aGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTg                                                                                                                                     <  2:147108/150‑1 (MQ=255)
            gCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGa                                                                                                                                    >  1:200379/1‑150 (MQ=255)
             ccAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAac                                                                                                                                  >  2:548299/1‑151 (MQ=255)
             ccAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAac                                                                                                                                  >  1:183218/1‑151 (MQ=255)
               aGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAaca                                                                                                                                 <  1:41255/150‑1 (MQ=255)
                 aaaaaGCTGAAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTa                                                                                                                              <  1:150134/151‑1 (MQ=255)
                  aaaaGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTa                                                                                                                              >  2:328969/1‑150 (MQ=255)
                    aaGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATcc                                                                                                                           >  1:833530/1‑151 (MQ=255)
                         gTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGGGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGaa                                                                                                                       <  2:602726/150‑1 (MQ=255)
                           aaaaGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGaaaa                                                                                                                     <  1:558979/150‑1 (MQ=255)
                                  ggTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGgc                                                                                                               <  1:116488/149‑1 (MQ=255)
                                    tgatgaGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGcgcg                                                                                                            <  1:382290/150‑1 (MQ=255)
                                    tgatgaGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGcgcg                                                                                                            >  1:704955/1‑150 (MQ=255)
                                        gaGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGcgtcgt                                                                                                        <  2:150518/150‑1 (MQ=255)
                                         aGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGttt                                                                                                      >  1:314755/1‑151 (MQ=255)
                                           tttGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGcgtcgt                                                                                                        <  2:703008/147‑1 (MQ=255)
                                           tttGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGtttt                                                                                                     <  2:901737/150‑1 (MQ=255)
                                            ttGCATTCCAGTATCGTCTGTACGGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTg                                                                                                    <  2:661721/150‑1 (MQ=255)
                                              gCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCg                                                                                                 >  1:202141/1‑151 (MQ=255)
                                                aTTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTc                                                                                               >  2:459998/1‑151 (MQ=255)
                                                 ttCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCg                                                                                              >  1:383322/1‑151 (MQ=255)
                                                 ttCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCg                                                                                              >  2:819744/1‑151 (MQ=255)
                                                  tCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCg                                                                                              >  1:535750/1‑150 (MQ=255)
                                                   ccAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGAt                                                                                            >  2:281953/1‑151 (MQ=255)
                                                      gTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTc                                                                                         <  1:806539/151‑1 (MQ=255)
                                                       tATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTc                                                                                         <  1:432932/150‑1 (MQ=255)
                                                                       gTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTg                                                                          <  1:787608/149‑1 (MQ=255)
                                                                          cgcAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGAAGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGc                                                                     >  1:389190/1‑151 (MQ=255)
                                                                          cgcAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCCGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTg                                                                          <  2:262800/146‑1 (MQ=255)
                                                                             aaCCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGc                                                                  >  1:403604/1‑151 (MQ=255)
                                                                                 gCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCcgcg                                                              >  2:288573/1‑151 (MQ=255)
                                                                                            gCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCATCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCAt                                                   <  2:223995/151‑1 (MQ=255)
                                                                                             cccATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCAtt                                                  <  1:141499/151‑1 (MQ=255)
                                                                                              ccATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTg                                                 >  1:807515/1‑151 (MQ=255)
                                                                                                   aGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCg                                            <  1:856707/151‑1 (MQ=255)
                                                                                                      gcgcgTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTg                                         <  1:563645/151‑1 (MQ=255)
                                                                                                       cgcgTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTg                                         >  2:652029/1‑150 (MQ=255)
                                                                                                               ggCGACGCGTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCt                                 <  2:383314/150‑1 (MQ=255)
                                                                                                                      cgTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTg                         <  1:212075/151‑1 (MQ=255)
                                                                                                                       gTACCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGGGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTgg                        >  2:309605/1‑151 (MQ=255)
                                                                                                                         aCCGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTgggg                      >  2:715481/1‑151 (MQ=255)
                                                                                                                           cGGCATGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGga                     >  2:580553/1‑150 (MQ=255)
                                                                                                                               aTGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGgaagcga                <  2:313054/151‑1 (MQ=255)
                                                                                                                               aTGGGCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGgaagcga                <  1:42361/151‑1 (MQ=255)
                                                                                                                               aTGGGCGGTTTCTCGGACGGTTGGGCGCCAGGGGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGCGATCTGAAAGCTGCCGCGCCCAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGgaa                    >  1:611431/1‑147 (MQ=255)
                                                                                                                                 gggCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGgaagcgaag              <  2:319011/151‑1 (MQ=255)
                                                                                                                                 gggCGGTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGgaagcgaa               <  1:130928/150‑1 (MQ=255)
                                                                                                                                     ggTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGGAAGCGAAGCaaa          <  2:833512/151‑1 (MQ=255)
                                                                                                                                      gTTTCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGGAAGCGAAGCAAAt         <  1:48040/151‑1 (MQ=255)
                                                                                                                                       tttCTCGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGGAAGCGAAGCAAATc        >  1:410770/1‑151 (MQ=255)
                                                                                                                                         tctcGGACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGGAAGCGAAGCAAATCg       >  2:833144/1‑150 (MQ=255)
                                                                                                                                             ggACGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGGAAGCGAAGCAAATCGAAAtt  >  2:231415/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCATTAGCGCGCGTTATGGCGACGCGTACCGGCATGGGCGGTTTCTCGGAAGGTTGGGCGCCAGGTGAACACTATCCCGAAAAATGGGCGCGTCGTTTTGCCGTCGATTTCGTTGGTGGTGATCTGAAAGCTGCCGCGCCAAAAGGCATTGAGCCGGTGATTACGCTTTCCAGTGGGGAAGCGAAGCAAATCGAAATT  >  NC_000913/1498028‑1498319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑