Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,421,944 G→C 5.5% G78R (GGC→CGC) ‡ moeA → molybdenum::molybdopterin ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,421,9440GC 5.5% 236.8 / 3.2 72G78R (GGC→CGC) ‡moeAmolybdenum::molybdopterin ligase
Reads supporting (aligned to +/- strand):  ref base G (38/30);  new base C (2/2);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCG  >  NC_002947/2421800‑2422091
                                                                                                                                                |                                                                                                                                                     
gCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGc                                                                                                                                                     >  7:167764/1‑147 (MQ=255)
gCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGc                                                                                                                                                     <  8:167764/147‑1 (MQ=255)
 cACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGcc                                                                                                                                                    <  8:54669/147‑1 (MQ=255)
 cACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGcc                                                                                                                                                    <  7:103385/147‑1 (MQ=255)
  aCCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAggg                                                                                                                                                      >  5:34155/1‑144 (MQ=255)
  aCCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAggg                                                                                                                                                      <  6:34155/144‑1 (MQ=255)
   ccGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGcc                                                                                                                                                <  5:87608/149‑1 (MQ=255)
      aTCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCt                                                                                                                                             >  3:74786/1‑149 (MQ=255)
         cgcgAGACCGAAACCCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                           <  5:301630/148‑1 (MQ=255)
          gcgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCg                                                                                                                                         >  7:16666/1‑149 (MQ=255)
              gaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTca                                                                                                                                    <  4:22168/149‑1 (MQ=255)
                 cGAAAACCTGCCGCTGGCCGAAGCCGATGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTcagc                                                                                                                                  >  7:62276/1‑148 (MQ=255)
                            cGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTAt                                                                                                                          >  3:175095/1‑145 (MQ=255)
                            cGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTAt                                                                                                                          <  4:175095/145‑1 (MQ=255)
                                  ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTc                                                                                                                      >  5:302478/1‑143 (MQ=255)
                                  ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTc                                                                                                                      <  6:302478/143‑1 (MQ=255)
                                  ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGc                                                                                                                >  6:134050/1‑149 (MQ=255)
                                  ccGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGc                                                                                                                >  8:239602/1‑149 (MQ=255)
                                   cGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGc                                                                                                                >  2:206928/1‑148 (MQ=255)
                                   cGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGc                                                                                                                >  4:100675/1‑148 (MQ=255)
                                    gAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTca                                                                                                                                    >  2:250570/1‑127 (MQ=38)
                                    gAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTca                                                                                                                                    <  1:250570/127‑1 (MQ=38)
                                          gAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGcc                                                                                                                                                >  3:15214/1‑110 (MQ=255)
                                          gAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGcc                                                                                                                                                <  4:15214/110‑1 (MQ=255)
                                            gggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                           <  7:106845/113‑1 (MQ=255)
                                            gggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                           >  8:106845/1‑113 (MQ=255)
                                                 gcgTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCg                                                                                                 >  1:205246/1‑149 (MQ=255)
                                                   gTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                               <  4:68539/149‑1 (MQ=255)
                                                   gTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                               >  7:252231/1‑149 (MQ=255)
                                                      cTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGc                                                                                             >  6:305732/1‑148 (MQ=255)
                                                         gCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTc                                                                                                                                     >  6:201883/1‑105 (MQ=255)
                                                         gCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTc                                                                                                                                     <  5:201883/105‑1 (MQ=255)
                                                           cACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCgg                                                                                       <  8:62276/149‑1 (MQ=255)
                                                           cACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCgg                                                                                       <  8:16666/149‑1 (MQ=255)
                                                           cACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCgg                                                                                       <  3:216273/149‑1 (MQ=255)
                                                              cGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                               <  5:30505/138‑1 (MQ=255)
                                                              cGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                               >  6:30505/1‑138 (MQ=255)
                                                              cGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTAc                                                                                   >  5:46513/1‑149 (MQ=255)
                                                                     gtggCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTgcgc                                                                             >  3:79691/1‑149 (MQ=255)
                                                                        gCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTgcgcgc                                                                           <  5:134050/148‑1 (MQ=255)
                                                                        gCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTgcgcgc                                                                           <  7:239602/148‑1 (MQ=255)
                                                                        gCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGcgcgcg                                                                          <  4:74786/149‑1 (MQ=255)
                                                                         ccTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCTCAGCCGCTGCAACCGGGTACCTGCGCGCgg                                                                         <  1:206928/149‑1 (MQ=255)
                                                                          ctcgctcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGa                                                                        >  4:240417/1‑149 (MQ=255)
                                                                          ctcgctcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGa                                                                        >  3:254797/1‑149 (MQ=255)
                                                                             gctcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCt                                                                    >  6:50916/1‑149 (MQ=255)
                                                                                 gACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTc                                                                  <  5:302907/147‑1 (MQ=255)
                                                                                         gccCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGgcgcgc                                                        >  8:60004/1‑149 (MQ=255)
                                                                                            cTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCtg                                                     >  3:67993/1‑149 (MQ=255)
                                                                                                 cGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCg                                                <  8:252231/149‑1 (MQ=255)
                                                                                                    aCAGCGCAATGGACGGTTATGCCTTGCGCCAGGCTGACTGGCAGGGCCAGCCGCTGACGGTCAGCCAGCGTATTTTCGCCGGCCATGAGCCCCAGCCGCTGAAACCGGGTACCTGCGCGCGCATCTTCACCGGCGCGCCTGTGCCGGAg                                             >  2:77615/1‑149 (MQ=255)
                                                                                                      aGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAggg                                           >  2:16100/1‑149 (MQ=255)
                                                                                                       gcgcCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTACAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAggg                                           >  1:94524/1‑149 (MQ=255)
                                                                                                           cATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCACCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGcc                                       >  2:90024/1‑148 (MQ=255)
                                                                                                                aCGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGc                                 <  4:67993/149‑1 (MQ=255)
                                                                                                                   gTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTc                              <  4:254797/149‑1 (MQ=255)
                                                                                                                     tATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGAc                                    <  2:145565/141‑1 (MQ=255)
                                                                                                                     tATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGAc                                    >  1:145565/1‑141 (MQ=255)
                                                                                                                         ccTTGCGCCTGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                        <  6:46513/148‑1 (MQ=255)
                                                                                                                          cTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                        >  4:181644/1‑148 (MQ=255)
                                                                                                                          cTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                        >  8:129331/1‑148 (MQ=255)
                                                                                                                             gcgcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAg                     <  7:60004/148‑1 (MQ=255)
                                                                                                                               gcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGa                   >  3:247860/1‑148 (MQ=255)
                                                                                                                                 cTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGaaaa                >  4:24595/1‑149 (MQ=255)
                                                                                                                                 cTGGCGGACTGGCAGGGCCAGCCGCTGCCGGGCATCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTTCCTGCGCGCGGATTATTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGaaa                 <  1:90024/149‑1 (MQ=255)
                                                                                                                                      tGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCGGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACAcc            >  7:217600/1‑148 (MQ=255)
                                                                                                                                      tGACCGGCAGCGGCTGGCCCTGCCAGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACAcc            >  2:50868/1‑148 (MQ=255)
                                                                                                                                        aCTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGaa         <  3:240417/149‑1 (MQ=255)
                                                                                                                                        aCTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGa          <  5:305732/148‑1 (MQ=255)
                                                                                                                                        aCTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGa          <  5:50916/148‑1 (MQ=255)
                                                                                                                                          tGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAg        >  1:142746/1‑148 (MQ=255)
                                                                                                                                            gCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTAc     >  3:117776/1‑149 (MQ=255)
                                                                                                                                            gCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTAc     >  6:283912/1‑149 (MQ=255)
                                                                                                                                              aGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTc   >  7:139969/1‑149 (MQ=255)
                                                                                                                                               gggCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCg  <  8:139969/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTCG  >  NC_002947/2421800‑2422091

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