Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 801,406 A→G E216E (GAA→GAG)  ybhH → putative PrpF family isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913801,4060AG98.4% 238.8 / NA 64E216E (GAA→GAG) ybhHputative PrpF family isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (29/34);  minor base T (0/1);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAACCAAAACCGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTG  >  NC_000913/801258‑801553
                                                                                                                                                    |                                                                                                                                                   
cTGGAACCAAAACCGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCt                                                                                                                                                   >  1:503528/1‑151 (MQ=255)
  ggAACCAAAACCGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTg                                                                                                                                                  <  1:101533/150‑1 (MQ=255)
    aaCCAAAACCGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGAt                                                                                                                                               <  1:487990/151‑1 (MQ=255)
    aaCCAAAACCGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGAt                                                                                                                                               >  1:998568/1‑151 (MQ=255)
     aCCAAAACCGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTgg                                                                                                                                                 >  1:196414/1‑148 (MQ=255)
        aaaaCCGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCg                                                                                                                                           <  1:1024494/151‑1 (MQ=255)
          aaCCGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGAc                                                                                                                                         >  1:528312/1‑151 (MQ=255)
            ccGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACaa                                                                                                                                       <  2:75395/151‑1 (MQ=255)
            ccGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACaa                                                                                                                                       <  1:766585/151‑1 (MQ=255)
               gAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGc                                                                                                                                    <  2:583714/151‑1 (MQ=255)
                  aaGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTCTTCCGGCTGAATATCTGGGAAAAACAGGTTATGACTTACCGGCGGAGCTGGATGCCGACAAAGCatt                                                                                                                                 <  2:785689/151‑1 (MQ=255)
                       ttCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCattatt                                                                                                                              >  1:754375/1‑149 (MQ=255)
                       ttCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCAttatta                                                                                                                             <  2:669859/150‑1 (MQ=255)
                        tCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGATCTGGATGCCGACAAAGCATTATTAg                                                                                                                            <  2:72376/150‑1 (MQ=255)
                             aCTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGc                                                                                                                       >  1:984816/1‑150 (MQ=255)
                                gATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTg                                                                                                                   <  1:801373/151‑1 (MQ=255)
                                     tCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCt                                                                                                              >  2:505839/1‑151 (MQ=255)
                                       aGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTAt                                                                                                            >  1:804057/1‑151 (MQ=255)
                                       aGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCAATATTAGCCCGCATTGAATCTAt                                                                                                            <  1:268871/151‑1 (MQ=255)
                                            gattATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGt                                                                                                        <  1:33753/150‑1 (MQ=255)
                                            gattATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGt                                                                                                        <  2:898295/150‑1 (MQ=255)
                                                  tttGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACa                                                                                                   >  1:719959/1‑149 (MQ=255)
                                                       cGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGTGGGt                                                                                            <  2:147755/151‑1 (MQ=255)
                                                        gATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCggg                                                                                             >  1:711808/1‑149 (MQ=255)
                                                        gATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGt                                                                                            >  1:717948/1‑150 (MQ=255)
                                                          tGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTAACGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTa                                                                                           <  2:849429/149‑1 (MQ=255)
                                                                    aCCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCt                                                                               <  1:210030/151‑1 (MQ=255)
                                                                       tGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCtt                                                                              >  2:751764/1‑149 (MQ=255)
                                                                        gTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTa                                                                             <  1:548060/149‑1 (MQ=255)
                                                                        gTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAg                                                                            <  1:563009/150‑1 (MQ=255)
                                                                          aTCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTg                                                                         >  2:969127/1‑151 (MQ=255)
                                                                              atatGGCGATGCCAGTCGTCATTATTCCGTCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTTGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGAAGt                                                                     >  1:124849/1‑151 (MQ=255)
                                                                                atGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGCGCCGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAATCAATTGGGATAGGGGATg                                                                      >  1:53926/1‑148 (MQ=255)
                                                                                atGGCGATGCCAGTCGTCACTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTc                                                                    >  2:327618/1‑150 (MQ=255)
                                                                                  ggCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGt                                                                 <  2:395706/151‑1 (MQ=255)
                                                                                    ctaTGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGt                                                                 >  2:615506/3‑149 (MQ=255)
                                                                                        gCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAatat                                                            >  1:979158/1‑150 (MQ=255)
                                                                                         ccAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATgg                                                          <  1:579107/151‑1 (MQ=255)
                                                                                           aGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGtt                                                        <  1:1034552/151‑1 (MQ=255)
                                                                                                 cATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCt                                                  <  2:511100/151‑1 (MQ=255)
                                                                                                     attCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTaaa                                               <  2:436758/150‑1 (MQ=255)
                                                                                                       tCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTaaa                                               >  2:17986/1‑148 (MQ=255)
                                                                                                         cGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCtgt                                          >  1:246664/1‑151 (MQ=255)
                                                                                                         cGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCtgt                                          <  2:284003/151‑1 (MQ=255)
                                                                                                                atatCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTAtt                                   <  1:868054/151‑1 (MQ=255)
                                                                                                                atatCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTAtt                                   <  2:156774/151‑1 (MQ=255)
                                                                                                                 tatCTGGGAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTAttt                                  >  1:899738/1‑151 (MQ=255)
                                                                                                                       ggAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAg                            <  2:711771/151‑1 (MQ=255)
                                                                                                                       ggAAAAACAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGGAAAACAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAAACTGTGCTTATTTCTcc                              >  2:134421/1‑149 (MQ=255)
                                                                                                                          aaaaCAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTcc                              >  2:176392/1‑146 (MQ=255)
                                                                                                                          aaaaCAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCa                             <  1:176393/147‑1 (MQ=255)
                                                                                                                            aaCAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAg                       <  2:401755/151‑1 (MQ=255)
                                                                                                                             aCAGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGGTTATTTCTCCAgc                           >  2:649655/1‑146 (MQ=255)
                                                                                                                               aGGTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGGGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGaa                     >  1:3016/1‑150 (MQ=255)
                                                                                                                                ggTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAg                   <  1:168791/151‑1 (MQ=255)
                                                                                                                                ggTTATGAATTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAg                   <  2:611029/151‑1 (MQ=255)
                                                                                                                                       aaTTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCgggg             >  1:20246/1‑150 (MQ=255)
                                                                                                                                        aTTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCa           <  1:292270/151‑1 (MQ=255)
                                                                                                                                        aTTACCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCa           <  1:975525/151‑1 (MQ=255)
                                                                                                                                           aCCGGCGGAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGGAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAAtt        >  1:108754/1‑151 (MQ=255)
                                                                                                                                             cggcggAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAAtt        >  2:497809/1‑149 (MQ=255)
                                                                                                                                             cggcggAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTaa      <  1:640141/151‑1 (MQ=255)
                                                                                                                                              ggcggAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAAt     >  1:339800/1‑151 (MQ=255)
                                                                                                                                                 ggAGCTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAAtgtg  <  2:533755/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGGAACCAAAACCGGAAAAGTTTTCCCGACTGATAATCAGATTGATTATTTTGACGATGTCCCGGTGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTG  >  NC_000913/801258‑801553

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