Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,953,408 T→C 8.0% T45A (ACC→GCC) ‡ mnmC ← methyltransferase/FAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,953,4080TC 8.0% 208.3 / 7.3 75T45A (ACC→GCC) ‡mnmCmethyltransferase/FAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base C (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 = 5.12e-01

GCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATGGTGGGCTCGGTT  >  NC_002947/1953259‑1953553
                                                                                                                                                     |                                                                                                                                                 
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                                                                                                                                 >  1:45962/1‑149 (MQ=255)
                   gcgTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTg                                                                                                                                 >  7:119675/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:268223/1‑149 (MQ=255)
                                                            gAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGAcgt                                                                                        >  8:282748/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                                                           <  2:42056/143‑1 (MQ=255)
                                                                                               cGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGccc                                                           >  1:42056/1‑143 (MQ=255)
                                                                                                    aCCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATggggg                                                                 <  2:17161/132‑1 (MQ=255)
                                                                                                    aCCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATggggg                                                                 >  1:17161/1‑132 (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                                                                                                            <  1:271178/68‑1 (MQ=255)
                                                                                                                         aGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTc                                                                                                            >  2:271178/1‑68 (MQ=255)
                                                                                                                           gTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTggagga                         <  7:268223/149‑1 (MQ=255)
                                                                                                                                cAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAg                                               >  2:95582/1‑122 (MQ=255)
                                                                                                                                 aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTgg                                                        <  1:67124/112‑1 (MQ=255)
                                                                                                                                 aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTgg                                                        >  2:67124/1‑112 (MQ=255)
                                                                                                                                 aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTgg                   >  5:321128/1‑149 (MQ=255)
                                                                                                                                 aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTgg                   >  1:304170/1‑149 (MQ=255)
                                                                                                                                      ccTCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTg                                                         >  6:247586/2‑106 (MQ=255)
                                                                                                                                      ccTCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTg                                                         <  5:247586/105‑1 (MQ=255)
                                                                                                                                       cTCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGggc                                                                >  5:29626/1‑98 (MQ=255)
                                                                                                                                       cTCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGggc                                                                <  6:29626/98‑1 (MQ=255)
                                                                                                                                          ggCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGc                                <  5:292893/127‑1 (MQ=255)
                                                                                                                                          ggCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGc                                >  6:292893/1‑127 (MQ=255)
                                                                                                                                                  gggTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATGGTGGGCTCGGtt  <  6:147511/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATGGTGGGCTCGGTT  >  NC_002947/1953259‑1953553

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: