Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 816,440 G→T 7.0% L70M (CTG→ATG)  ybhK ← putative CofD superfamily transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913816,4400GT 7.0% 180.2 / 2.4 66L70M (CTG→ATG) ybhKputative CofD superfamily transferase
Reads supporting (aligned to +/- strand):  ref base G (23/37);  new base T (2/3);  total (26/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGAATTAAATTGATGGCTTCCAGAGGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTC  >  NC_000913/816292‑816583
                                                                                                                                                    |                                                                                                                                                
cGAATTAAATTGATGGCTTCCAGAGGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAg                                                                                                                                                  <  1:997457/149‑1 (MQ=255)
 gAATTAAATTGATGGCTTCCAGAGGCCGCAGGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTGCGCCATTGCCACGAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGCTTGCGGTagcagc                                                                                                                                                 >  2:892928/1‑149 (MQ=255)
 gAATTAAATTGATGGCTTCCAGAGGCCGCAAGCAGCATCGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGc                                                                                                                                                 >  1:634914/1‑149 (MQ=255)
 gAATTAAATCGCTGGCTTCCAGAGGCAGCACGCTAAGGGGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGc                                                                                                                                                 <  1:552703/149‑1 (MQ=255)
    ttAAATTGATGGCTTCCAGAGGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCTAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTgg                                                                                                                                              >  2:849652/1‑149 (MQ=255)
      aaaTTGATGGCTTCCAGAGGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGtt                                                                                                                                            <  1:1111473/149‑1 (MQ=255)
      aaaTGGATGGCTTCCAGAGGCAGCACGCTAAGGTGATCCAGCGCCGTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCACAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGtt                                                                                                                                            <  2:554319/149‑1 (MQ=255)
        aTTGATGGCTTCCAGAGGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACTGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGa                                                                                                                                          <  1:14072/149‑1 (MQ=255)
        aTTGATGGCTTCCAGAGGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGa                                                                                                                                          <  2:301941/149‑1 (MQ=255)
             tGGCTTCCAGAGGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAg                                                                                                                                     >  2:30213/1‑149 (MQ=255)
                   ccAGAGGGCGCCAGCTAAGTTGCTCCAGCGCCTTTAACATCAAGTTTCCGAGCTTCTGACCGGAAAGTGCGCCATTGCCCCCAAAACGGTATTCACACATCGCGGAGGCGTCGCTCGGTTCCGTTATCAGCTGGTTGCGGCAGTTgcgc                                                                                                                               <  2:1065290/149‑1 (MQ=255)
                     agagGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCat                                                                                                                             >  2:17195/1‑149 (MQ=255)
                      gagGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCCTTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCACCTGGTTGAGGCAGTTGCGCcta                                                                                                                            >  1:1038018/1‑149 (MQ=255)
                      gagGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCata                                                                                                                            >  1:585361/1‑149 (MQ=255)
                      gagGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTAATAAGCTGGTTGAGGCAGTTGCGCata                                                                                                                            >  1:495275/1‑149 (MQ=255)
                       agGCCGCACGCTAAGGGGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCatat                                                                                                                           <  2:213405/149‑1 (MQ=255)
                               cGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAg                                                                                                                  <  1:1232615/149‑1 (MQ=255)
                                    aggggctCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGATCGGAAATTTCGCCATTGCCACCAACACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAAt                                                                                                             <  2:475472/143‑1 (MQ=255)
                                     ggTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGCTTGAGGCAGTTGCGCATATCGCCCCGGGCAATg                                                                                                            >  2:90296/1‑149 (MQ=255)
                                      gTGATACAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGGTCCGTTATCAGCTGGGTGAGGCAGTTGCGCATATCGCCACAGACaagtc                                                                                                           >  2:802562/1‑146 (MQ=255)
                                       tGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCATCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCGATgcc                                                                                                          >  1:238437/1‑149 (MQ=255)
                                       tGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCATCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATgcc                                                                                                          >  1:269634/1‑149 (MQ=255)
                                         aTCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCGGTTGCGCATATCGCCCCAGGCAATgccgc                                                                                                        <  1:791868/149‑1 (MQ=255)
                                            cAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAACCGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTCAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCtt                                                                                                     <  1:839782/149‑1 (MQ=255)
                                              gcgcCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCCGGCAATGCCGCCTTCt                                                                                                   >  2:432625/1‑149 (MQ=255)
                                              gcgcCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCt                                                                                                   >  1:216299/1‑149 (MQ=255)
                                                 catttAACCTCGAGGTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGaa                                                                                                <  1:1133046/147‑1 (MQ=255)
                                                   tttAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAGACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACg                                                                                              >  1:845458/1‑149 (MQ=255)
                                                   tttAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACg                                                                                              >  1:817150/1‑149 (MQ=255)
                                                     tAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGc                                                                                            <  1:735723/149‑1 (MQ=255)
                                                      aaCATCAAGTTGCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCg                                                                                           <  2:906906/149‑1 (MQ=255)
                                                               tttCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGATATCATCTTGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGccc                                                                                  <  1:1214629/149‑1 (MQ=255)
                                                                  ccGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTTAACGGCGAATACGCCCCGt                                                                               >  1:1136357/1‑149 (MQ=255)
                                                                   cGAGATTATGACCGGAGAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGGTGAGGCAGTTGCGGATATCGCCCCAGGCAATGCCGCCTTCGGCACGGCGAATACGCCGCGTc                                                                              >  2:759226/1‑149 (MQ=255)
                                                                   cGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGACGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTc                                                                              >  2:959555/1‑149 (MQ=255)
                                                                    gagaTCATGACCGGAAAGTTCGCCCTCGCCACCAAAACGGTATTCAACCATCGCGGAGGCGCCGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCg                                                                             <  1:70132/149‑1 (MQ=255)
                                                                     agaTTATGACCGGAAAGTTCGCCATGGCCACAAAAACGGTATGCAAACATCGCGGAGGCGACGCTCGGGTCCGTTATCAGCTGTTAGATGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGa                                                                            <  1:249464/149‑1 (MQ=255)
                                                                       aTTAGGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGc                                                                          <  2:551879/149‑1 (MQ=255)
                                                                             aCCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGGCGGGCCACCAt                                                                    >  2:815626/1‑149 (MQ=255)
                                                                                     gtccGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGACAATGCTCGGTTCCGTTATCAGCTGGTATGTGCACTGCACCCATGCGTCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCggtg                                                            <  1:629151/146‑1 (MQ=255)
                                                                                        cGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCggtggtg                                                         <  1:892928/149‑1 (MQ=255)
                                                                                        cGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCggtggtg                                                         <  1:30213/149‑1 (MQ=255)
                                                                                        cGCCATTGCCACCAAAACGGCACTCAAACATCGCGGGGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCggtggtg                                                         <  2:1011131/149‑1 (MQ=255)
                                                                                             ttGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTCCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGAc                                                    <  1:50854/149‑1 (MQ=255)
                                                                                               gCCACCAGAACGGTATTCAAACATCGCGGAGGCGACGATCGGTTCCGTTATCAGCTGGTTGAGGCAGTTACGCATATCGCCCCAGGCAATGCAGCCTTCTCAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGa                                                  <  2:946127/149‑1 (MQ=255)
                                                                                               gCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGCGCCACCATTATCGGTGGTGGTGACGa                                                  <  2:634914/149‑1 (MQ=255)
                                                                                                  accaAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGGCCCGTCGAGCCACTATTATCGGTGGTGGTGACGATAc                                               <  1:849652/149‑1 (MQ=255)
                                                                                                  accaAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCCTTTTCGGCGGCGGTGagggtgc                                               >  2:841456/1‑143 (MQ=255)
                                                                                                     aaaaCGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCg                                            >  2:1119533/1‑149 (MQ=255)
                                                                                                     aaaaCGGTATTCAAACATCGCCGAGGCGCCGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTGGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCg                                            <  2:681809/149‑1 (MQ=255)
                                                                                                          ggTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGACCGGCGAATACGCCCCGTCGAGCCACCATTATCGGTAGTGGTGACGATACACGTTaaa                                       >  1:35226/1‑149 (MQ=255)
                                                                                                          ggTATTCAAACATCGCAGAGGCGACGCTCAGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGAGGTGGTGACGATACCCGTTaaa                                       >  1:821512/1‑149 (MQ=255)
                                                                                                          ggTATTCAAACACCCCGGAGGCGCCGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGGTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTaaa                                       <  1:30055/149‑1 (MQ=255)
                                                                                                                  aaCATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGCTGACGATACCCGTTAAACGAGAAcc                               <  1:652462/149‑1 (MQ=255)
                                                                                                                      gcgcgGAGGCGACGCTCGGTTCCGTTATCAGGTGGTTGAGGAAGTTGCGCTTAACGCACCGTGCAATGCCGAATTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCaaa                           <  2:596498/148‑1 (MQ=255)
                                                                                                                       cgcgGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGGCGATGCCCGTTAAACGGGAACCCAAAg                          >  1:374119/1‑149 (MQ=255)
                                                                                                                              gggACTCTCGGTTCCCTTATGATCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAg                   <  2:238437/147‑1 (MQ=255)
                                                                                                                              gCGACGCTCGGTTCCGTTATCATCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAg                   <  2:269634/149‑1 (MQ=255)
                                                                                                                                 aCGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGtga                <  1:902205/149‑1 (MQ=255)
                                                                                                                                 aCGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGtga                <  1:1041854/149‑1 (MQ=255)
                                                                                                                                  cGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGTATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATAATCCAAGGTATTTACGATACCCGTTAAACGAGAACCCAAAGACGAAAGtgat               <  1:629902/149‑1 (MQ=255)
                                                                                                                                  cGCCCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAAGGCCGCCTTCTGAACGGCGAATACGCCCCGCAATGTCGCCATTTTCGGTCGAGTTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGtgat               <  1:629869/149‑1 (MQ=255)
                                                                                                                                    cTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGtgatga             <  2:718222/149‑1 (MQ=255)
                                                                                                                                     tCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATgag            >  1:580729/1‑149 (MQ=255)
                                                                                                                                     tCGGTCCCGTTATCGGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATgag            <  1:435151/149‑1 (MQ=255)
                                                                                                                                              ttATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATCCGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGt   <  2:1022950/149‑1 (MQ=255)
                                                                                                                                              ttATCAGCTGGTTGAGGCAGTTGCGCATACCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTTTCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGt   <  1:1078542/149‑1 (MQ=255)
                                                                                                                                               tATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATAAGCCCCGTTGAGTCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTc  >  2:328390/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGAATTAAATTGATGGCTTCCAGAGGCCGCACGCTAAGGTGATCCAGCGCCTTTAACATCAAGTTTCCGAGATTATGACCGGAAAGTTCGCCATTGCCACCAAAACGGTATTCAAACATCGCGGAGGCGACGCTCGGTTCCGTTATCAGCTGGTTGAGGCAGTTGCGCATATCGCCCCAGGCAATGCCGCCTTCTGAACGGCGAATACGCCCCGTCGAGCCACCATTATCGGTGGTGGTGACGATACCCGTTAAACGAGAACCCAAAGACGAAAGTGATGAGAGAACGCGTC  >  NC_000913/816292‑816583

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