Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,953,403 G→A 8.0% R46R (CGC→CGT) ‡ mnmC ← methyltransferase/FAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,953,4030GA 8.0% 217.8 / 6.0 75R46R (CGC→CGT) ‡mnmCmethyltransferase/FAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base A (2/4);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

TCGACGCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGG  >  NC_002947/1953254‑1953536
                                                                                                                                                     |                                                                                                                                     
tcGACGCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAgg                                                                                                                                       >  3:11503/1‑149 (MQ=255)
     gCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGt                                                                                                                                  <  4:11503/149‑1 (MQ=255)
        gACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCt                                                                                                                                >  8:37231/1‑149 (MQ=255)
              gTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGACCCGCCTGGCCGa                                                                                                                          <  8:295819/149‑1 (MQ=255)
              gTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGACCCGCCTGGCCGa                                                                                                                          <  5:105782/149‑1 (MQ=255)
               tGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAg                                                                                                                         <  7:37231/149‑1 (MQ=255)
                 cAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAgg                                                                                                                        <  6:237917/148‑1 (MQ=255)
                  aGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGaa                                                                                                                      >  1:313583/1‑149 (MQ=255)
                     cgcgcgTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAacac                                                                                                                   <  2:313583/149‑1 (MQ=255)
                       cgcgTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGt                                                                                                                 >  4:272763/1‑149 (MQ=255)
                       cgcgTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGt                                                                                                                 >  5:104367/1‑149 (MQ=255)
                        gcgTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTg                                                                                                                >  7:119675/1‑149 (MQ=255)
                        gcgTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTg                                                                                                                >  1:45962/1‑149 (MQ=255)
                          gTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCt                                                                                                              <  8:253653/149‑1 (MQ=255)
                             gCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTg                                                                                                            <  8:119675/148‑1 (MQ=255)
                              cTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGt                                                                                                          >  5:118402/1‑149 (MQ=255)
                              cTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGt                                                                                                          >  4:328852/1‑149 (MQ=255)
                               tgtgtgCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTg                                                                                                            >  2:246122/1‑146 (MQ=255)
                               tgtgtgCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTCGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTg                                                                                                            <  1:246122/146‑1 (MQ=255)
                                   tgcgtgcTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTctct                                                                                                      <  2:45962/148‑1 (MQ=255)
                                     cgtgcTGGTCTAACAGCTGCCAGTCTCAGAAAAAATTCATGCCGGTGCCAAAAGCTGTTTCGCCTATCACCGTGCAGGCGTGCTGCGCCAGGTTGGCAAAACGCTCGGCCATGCGGGTCTTTCGGAGGAACTCGTGGTTGGTTTCTTCg                                                                                                   <  6:26714/149‑1 (MQ=255)
                                         cTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACg                                                                                                                  >  2:196197/1‑130 (MQ=255)
                                         cTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACg                                                                                                                  <  1:196197/130‑1 (MQ=255)
                                          tGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGcc                                                                                              >  5:169467/1‑149 (MQ=255)
                                            gTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCtt                                                                                            >  4:286006/1‑149 (MQ=255)
                                             tCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTc                                                                                           >  6:252964/1‑149 (MQ=255)
                                             tCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTc                                                                                           <  3:328852/149‑1 (MQ=255)
                                               gAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCg                                                                                          <  3:272763/148‑1 (MQ=255)
                                                   aGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGtt                                                                                        <  6:92139/146‑1 (MQ=255)
                                                   aGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGtt                                                                                        >  5:92139/1‑146 (MQ=255)
                                                                 gAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGAcgt                                                                       >  8:282748/1‑149 (MQ=255)
                                                                 gAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGAcgt                                                                       >  8:268223/1‑149 (MQ=255)
                                                                        ttCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATa                                                                <  6:169467/149‑1 (MQ=255)
                                                                           aTGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGTGGGTGTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTg                                                             >  7:269667/1‑149 (MQ=255)
                                                                           aTGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGt                                                              <  7:282748/148‑1 (MQ=255)
                                                                               cGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCcgc                                                         <  6:104367/149‑1 (MQ=255)
                                                                               cGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCcgc                                                         >  8:116953/1‑149 (MQ=255)
                                                                                 gTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACt                                                                                    <  6:61017/120‑1 (MQ=255)
                                                                                 gTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACt                                                                                    >  5:61017/1‑120 (MQ=255)
                                                                                     cAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATgg                                                   <  5:252964/149‑1 (MQ=255)
                                                                                      aaaaCCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATggg                                                  >  4:315832/1‑149 (MQ=255)
                                                                                      aaaaCCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAAAcgt                                                                       >  8:219533/1‑128 (MQ=255)
                                                                                                  gCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGTGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTgg                                       <  8:269667/148‑1 (MQ=255)
                                                                                                  gCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgt                                      >  5:147511/1‑149 (MQ=255)
                                                                                                    cGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGccc                                          >  1:42056/1‑143 (MQ=255)
                                                                                                    cGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGccc                                          <  2:42056/143‑1 (MQ=255)
                                                                                                         aCCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATggggg                                                >  1:17161/1‑132 (MQ=255)
                                                                                                         aCCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATggggg                                                <  2:17161/132‑1 (MQ=255)
                                                                                                         aCCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTcc                                 <  3:315832/147‑1 (MQ=255)
                                                                                                          ccGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCTGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTcc                                 >  7:50811/1‑146 (MQ=255)
                                                                                                          ccGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCTGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCa                               <  8:50811/148‑1 (MQ=255)
                                                                                                          ccGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATa                                                                            >  5:96986/1‑103 (MQ=255)
                                                                                                          ccGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATa                                                                            <  6:96986/103‑1 (MQ=255)
                                                                                                          ccGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATgggg                                                 <  3:256243/130‑1 (MQ=255)
                                                                                                          ccGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATgggg                                                 >  4:256243/1‑130 (MQ=255)
                                                                                                           cGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGt                             >  1:135342/1‑149 (MQ=255)
                                                                                                            gTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTc                            <  3:286006/149‑1 (MQ=255)
                                                                                                                  gCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgtc                                     <  1:312711/134‑1 (MQ=255)
                                                                                                                  gCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgtc                                     >  2:312711/1‑134 (MQ=255)
                                                                                                                   cGTGCGGCGCCAGGTGGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTgg                     <  4:251862/149‑1 (MQ=255)
                                                                                                                         gcgcCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGc               <  6:118402/149‑1 (MQ=255)
                                                                                                                             cAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTgga           <  7:116953/149‑1 (MQ=255)
                                                                                                                              aGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTc                                                                                           >  2:271178/1‑68 (MQ=255)
                                                                                                                              aGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTc                                                                                           <  1:271178/68‑1 (MQ=255)
                                                                                                                                gTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTggagga        <  7:268223/149‑1 (MQ=255)
                                                                                                                                     cAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAg                              >  2:95582/1‑122 (MQ=255)
                                                                                                                                      aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTgg                                       >  2:67124/1‑112 (MQ=255)
                                                                                                                                      aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTgg                                       <  1:67124/112‑1 (MQ=255)
                                                                                                                                      aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTgg  >  1:304170/1‑149 (MQ=255)
                                                                                                                                      aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTgg  >  5:321128/1‑149 (MQ=255)
                                                                                                                                           ccTCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTg                                        <  5:247586/105‑1 (MQ=255)
                                                                                                                                           ccTCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTg                                        >  6:247586/2‑106 (MQ=255)
                                                                                                                                            cTCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGggc                                               <  6:29626/98‑1 (MQ=255)
                                                                                                                                            cTCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGggc                                               >  5:29626/1‑98 (MQ=255)
                                                                                                                                               ggCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGc               >  6:292893/1‑127 (MQ=255)
                                                                                                                                               ggCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGc               <  5:292893/127‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                     
TCGACGCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGG  >  NC_002947/1953254‑1953536

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