Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 165,039 G→C 9.5% A104P (GCG→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,0390GC 9.5% 156.6 / 8.1 61A104P (GCG→CCG) mrcBfused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (21/34);  new base C (4/2);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATT  >  NC_000913/164894‑165185
                                                                                                                                                 |                                                                                                                                                  
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:547012/149‑1 (MQ=255)
                                                                                                                                     tGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGc            >  2:696754/1‑149 (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  >  1:291104/1‑149 (MQ=255)
                                                                                                                                               ctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCAtt  <  2:2458/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATT  >  NC_000913/164894‑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: ‑