Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 165,038 T→G 10.8% P103P (CCT→CCG mrcB → fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913165,0380TG10.8% 127.9 / 5.4 61P103P (CCT→CCGmrcBfused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (20/34);  new base G (5/2);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

GCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATT  >  NC_000913/164890‑165185
                                                                                                                                                    |                                                                                                                                                   
gcaaaggGCGTACGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAActgcct                                                                                                                                                     <  2:640632/149‑1 (MQ=255)
gcaaaggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAActgcct                                                                                                                                                     >  1:372156/1‑149 (MQ=255)
    aggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCGATCGTTTTTGCCGTGCTGATCGCCATTTCCGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCgg                                                                                                                                                 <  1:779080/149‑1 (MQ=255)
    aggGCGTAAGCCTCGGGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACGGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGAGGGCAAGGTCTGGCAACTGCCTGCgg                                                                                                                                                 >  2:538119/1‑149 (MQ=255)
    aggGCGTAAGCCGCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTTTCGCTTTTGCCGTGCTGTTCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCgg                                                                                                                                                 <  1:422150/149‑1 (MQ=255)
        cGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGt                                                                                                                                             <  2:372156/149‑1 (MQ=255)
         gTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGGGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGCAGGCAGtt                                                                                                                                            <  2:572599/149‑1 (MQ=255)
             gCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATg                                                                                                                                        <  1:221885/149‑1 (MQ=255)
              ccTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATgg                                                                                                                                       <  1:970607/149‑1 (MQ=255)
              ccTCGTGGCAAACGCGGCTGGCTATGGCTACAGCTAAAACTGGCAATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGAATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATgg                                                                                                                                       <  2:571258/149‑1 (MQ=255)
               cTCGTGGCAAACGCGGCTGGCTATGGCAACTGCTAAAACTGGATATCGTTTTTGCCGTGCTGATCGACATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGc                                                                                                                                      <  2:938284/149‑1 (MQ=255)
                      cAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGGTTTTGCCGGGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATgg                                                                                                                               <  2:557375/149‑1 (MQ=255)
                       aaaCGCGGCTGGCTATGGCTACTGCTAAAACTGGCTCTCGTTTTTGCCGTGCTGAACGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGt                                                                                                                              <  2:567421/149‑1 (MQ=255)
                       aaaCGCGGCGGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGt                                                                                                                              <  2:958798/149‑1 (MQ=255)
                        aaCGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTc                                                                                                                             >  1:547012/1‑149 (MQ=255)
                        aaCGCCGCTGGCCACGCCTACTGCCAACACTCGCAAGCGCTTTTCGCGGGCCGATCGCCAATTACGGCGTGTAACTCGATCAAAAACTTCGTAGCCGTATTGCTGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTc                                                                                                                             <  2:943953/149‑1 (MQ=255)
                                tggctttggctaCTGCTAAAACTGGCTATCGTTTTTGCCGGGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGTCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGa                                                                                                                     <  1:138483/149‑1 (MQ=255)
                                    tatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCa                                                                                                                 >  1:176/1‑149 (MQ=255)
                                    tatggctaCTGCTAAAACCGGCTATCGTTTTTGCCGTGCTCATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCa                                                                                                                 <  2:458144/149‑1 (MQ=255)
                                           aCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGa                                                                                                          <  1:553506/149‑1 (MQ=255)
                                               cTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCCTAGCCGTAATGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCAt                                                                                                      <  1:801880/149‑1 (MQ=255)
                                               cTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTAACTCGATCAAAAAATACGTAGCCGTATTGATCGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTAAATCTTGAGCCAGACATGACCAt                                                                                                      >  1:479135/1‑149 (MQ=255)
                                                       ggCTATCGTTTTTGCCGTTCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTAAGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGa                                                                                              <  1:52231/149‑1 (MQ=255)
                                                       ggCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGa                                                                                              <  1:909171/149‑1 (MQ=255)
                                                            tCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACgag                                                                                         >  1:101347/1‑149 (MQ=255)
                                                              gTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGAt                                                                                       <  1:462942/149‑1 (MQ=255)
                                                               tttttGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATg                                                                                      <  2:240748/149‑1 (MQ=255)
                                                                   tGCCGTGCTGATCGCCCTTTACCGCGTTTATGTCGATCAAAAAATTCGTAGCCGTATTTATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAAGGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGa                                                                                  <  2:607775/149‑1 (MQ=255)
                                                                      cGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAAGCTTGAGCCAGCCATGACCATCAGCAAGAACGAGATGGTGAAgc                                                                               >  2:618247/1‑149 (MQ=255)
                                                                       gTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGTAAGTTTATGGCCGAATGGTCAATCTTGGGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgct                                                                              >  1:623793/1‑149 (MQ=255)
                                                                                 ccATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCg                                                                    >  1:488267/1‑149 (MQ=255)
                                                                                 ccATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCg                                                                    >  2:893477/1‑149 (MQ=255)
                                                                                   aTTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGAc                                                                  >  1:139527/1‑149 (MQ=255)
                                                                                       aCGGCGTTTATCTCGATCAAAACATTCGAAGCCGTATTGATGGCACGTCCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGGCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAg                                                              <  1:459515/149‑1 (MQ=255)
                                                                                       aCGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGATCCAGACATGACAATCAGCAAGAACGAGATGATGAAGCTGCTGGAGGCGACCCAg                                                              >  2:684117/1‑149 (MQ=255)
                                                                                           cGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATc                                                          <  1:286095/149‑1 (MQ=255)
                                                                                           cGTTTATCTCGAACAAAAAATTCGTAGCCGTATTGATGGCAAGTTCTGGCAACTTCCTGAGTCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATc                                                          <  1:652302/149‑1 (MQ=255)
                                                                                                aTCTCTATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATTGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAg                                                     <  1:511107/149‑1 (MQ=255)
                                                                                                  ctcGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGgt                                                   <  1:627103/149‑1 (MQ=255)
                                                                                                    cGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGgtgt                                                 >  1:641237/1‑149 (MQ=255)
                                                                                                     gATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTc                                                <  1:893477/149‑1 (MQ=255)
                                                                                                     gATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTc                                                >  2:56565/1‑149 (MQ=255)
                                                                                                      aTCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCg                                               >  1:201984/1‑149 (MQ=255)
                                                                                                         aaaaaaTGATTAGACTTATTGATGGCATGGACTGGCAACTGCATGCGGCAGGTTATGGCCGAAGGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGAGTGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGaaa                                            <  2:88190/149‑1 (MQ=255)
                                                                                                              atacGTAGCCGTATTGATGTCAAGGGCTGGCAACTGCCTGCGGCAGTTTTTGGCCGAATGGTCAATCTTGAGCCAGACATGAGCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGAc                                       <  1:273480/146‑1 (MQ=255)
                                                                                                                   tAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTAAATCTTGAGCCAGCCATTACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTc                                  >  2:383316/1‑149 (MQ=255)
                                                                                                                    aGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAATTGACCCGTcc                                 >  1:981044/1‑149 (MQ=255)
                                                                                                                    aGCCGTATTGATGCCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGGGAAGCTGCTGGATGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTcc                                 <  2:192566/149‑1 (MQ=255)
                                                                                                                      ccGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTg                               >  1:72431/1‑149 (MQ=255)
                                                                                                                               tGGCTAGGTCTGGCAACTGCCTGCGGCAGATTATGGCCGAATGGTGAATGTTGAGCCAGACATGACCAGCAGCAAGAACGCGATGGTGAAGCAGCTTGAGGCGACCCACTATGGTCAGGTGTCGAAAATGACCCGGCCTGGCGAATTTa                      <  2:1008210/149‑1 (MQ=255)
                                                                                                                                ggCAAGGTCTGGCAACTGCCTGCGGCGGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATAAGCAAGAACGAGATGGTGCAGCAGCTGGAGGCGACCCAGAATCGTCAGGTGTCGAAAATGCCCCGTCCAGGCGAATTTAc                     >  2:467850/1‑149 (MQ=255)
                                                                                                                                 gCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTAcc                    <  1:841400/149‑1 (MQ=255)
                                                                                                                                        cTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAgg             <  1:923777/149‑1 (MQ=255)
                                                                                                                                         tGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGc            >  2:696754/1‑149 (MQ=255)
                                                                                                                                         tGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGc            <  2:547012/149‑1 (MQ=255)
                                                                                                                                         cataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGc            >  1:882301/5‑149 (MQ=255)
                                                                                                                                         cataaaCTGCCGCAGGCAGGTTATGGCCGAATGGTCAATCTTGATCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCAAGTATCTTCAGGTGTCGAAAATGCCCCGTACTGGCGAATTAACCGTGCAGGc            >  2:789335/5‑149 (MQ=255)
                                                                                                                                          ataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGCTGGTGAAGCTGCTGGAGGCGACACAGGATCGTCAGGTGTAGAAAATGACCAGACCTAGCGAAATTACCCTGCAGGcc           >  2:1005369/4‑149 (MQ=255)
                                                                                                                                          ataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAAt                                          >  2:65191/4‑118 (MQ=255)
                                                                                                                                          ataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAAt                                          <  1:65191/115‑1 (MQ=255)
                                                                                                                                           gCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGATAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCa          >  2:358068/1‑149 (MQ=255)
                                                                                                                                            cAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCATACATGACCATCAGCAAGAACGAGATGGTGAAGCTCCTGGAGGCGACCCAGTATCGGCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCaa         <  2:602311/149‑1 (MQ=255)
                                                                                                                                              aCTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGGGAAGCTGCTGGAGGCGACCCAGTTTCGTCAGGTGTCGAAAATGCCCCGTCCTGGCGAAGTTACCGTGCAGGccacca       >  1:73984/1‑149 (MQ=255)
                                                                                                                                               ctgcctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAg      >  2:822912/1‑149 (MQ=255)
                                                                                                                                               ctgcctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACCAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAg      >  2:526134/1‑149 (MQ=255)
                                                                                                                                                tgcctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGc     >  2:860214/1‑149 (MQ=255)
                                                                                                                                                tgcctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTAGAAAATGACCCGTCCTGGCGTATTTACCGTGCAGGCCAACAGc     >  1:55679/1‑149 (MQ=255)
                                                                                                                                                 gcctgcGGCTGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCCTCAGCAAGTACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCa    >  1:390343/1‑149 (MQ=255)
                                                                                                                                                   ctgcGGCTGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCGGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCAtt  <  1:240272/149‑1 (MQ=255)
                                                                                                                                                   ctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCAtt  <  2:2458/149‑1 (MQ=255)
                                                                                                                                                   ctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCAtt  >  1:291104/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAAAGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATT  >  NC_000913/164890‑165185

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