Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,473,810 T→C 5.6% N329D (AAC→GAC) ‡ proA ← gamma‑glutamyl phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,473,8100TC 5.6% 212.6 / 2.5 71N329D (AAC→GAC) ‡proAgamma‑glutamyl phosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base C (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

CATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCC  >  NC_002947/5473666‑5473955
                                                                                                                                                |                                                                                                                                                  
cATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAgg                                                                                                                                                >  1:188952/1‑149 (MQ=255)
 aTCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGc                                                                                                                                               >  3:102576/1‑149 (MQ=255)
      cAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCg                                                                                                                                          >  6:305549/1‑149 (MQ=255)
      cAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGACATGAACTGACTGGCTTCAACCTCGTGTTCTATGATTATCTAGTCGGTGTGGTGAGAGCCATATTGGTTGATGTGCTCGATGGCCTGGGTCATGACGTc                                                                                                                                           >  6:349350/1‑148 (MQ=255)
       aaaGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGa                                                                                                                                         <  6:36037/149‑1 (MQ=255)
       aaaGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGACGGCCTGAACCAGGCCATCGa                                                                                                                                         <  8:92774/149‑1 (MQ=255)
        aaGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGGCAGAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGag                                                                                                                                        >  4:264898/1‑148 (MQ=255)
         atcAGGTCGTGGTGTTTTTCATGAACGTCGCCTAATCGATCTCAGCCATTTACTTACGGTCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTTGTGCGAGCCATAGTGGTTGATGTGCTCGACGGCCTGAACCAGGCCATCGagc                                                                                                                                       <  8:231466/147‑3 (MQ=255)
         aGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGAcc                                                                                                                                       <  4:340316/149‑1 (MQ=255)
         aGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGGGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGAcc                                                                                                                                       <  7:331553/149‑1 (MQ=255)
             ggTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCAcgc                                                                                                                                   >  2:128583/1‑149 (MQ=255)
             ggTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCAcgc                                                                                                                                   <  3:356891/149‑1 (MQ=255)
                 ggggTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAAt                                                                                                                               <  5:305549/149‑1 (MQ=255)
                    gtgtTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGa                                                                                                                            <  5:363851/149‑1 (MQ=255)
                    gtgtTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTAGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGa                                                                                                                            >  1:293416/1‑149 (MQ=255)
                       ttGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACa                                                                                                                          <  3:263305/148‑1 (MQ=255)
                       ttGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACa                                                                                                                          <  2:293416/148‑1 (MQ=255)
                          aGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGa                                                                                                                       <  5:175829/148‑1 (MQ=255)
                                 ccGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCg                                                                                                               <  7:282502/149‑1 (MQ=255)
                                  cGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATc                                                                                                                >  3:309056/1‑147 (MQ=255)
                                   gACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAg                                                                                                             <  1:128583/149‑1 (MQ=255)
                                      gCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTa                                                                                                          <  4:102576/149‑1 (MQ=255)
                                      gCCGAATCGACCTCAGCCATGAACTGACGGGCGTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACACGATCGCATCGAGGTa                                                                                                          <  4:44093/149‑1 (MQ=255)
                                        cGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTAtt                                                                                                        >  1:47435/1‑149 (MQ=255)
                                                   cAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAg                                                                                             >  3:10929/1‑149 (MQ=255)
                                                     gCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATTGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTc                                                                                           <  1:243365/149‑1 (MQ=255)
                                                           aaCTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGc                                                                                                <  2:95159/138‑1 (MQ=255)
                                                           aaCTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGc                                                                                                >  1:95159/1‑138 (MQ=255)
                                                             cTGACGGGCTTCACCCTGGTGTTCGCTGGTGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCgg                                                                                   <  4:228369/149‑1 (MQ=255)
                                                                aCGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTg                                                                                 <  2:188952/148‑1 (MQ=255)
                                                                 cGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTTGGTGGc                                                                               <  4:184104/149‑1 (MQ=255)
                                                                           ccTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGgcg                                                                     >  6:234529/1‑149 (MQ=255)
                                                                               gtgtTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtga                                                                 <  1:313440/149‑1 (MQ=255)
                                                                               gtgtTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtga                                                                 <  5:187714/149‑1 (MQ=255)
                                                                                tgtTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTTTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtgat                                                                <  7:56976/149‑1 (MQ=255)
                                                                                  ttCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtgatga                                                              <  7:216944/149‑1 (MQ=255)
                                                                                  ttCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCggatga                                                              >  5:227403/1‑149 (MQ=255)
                                                                                     gCTGTTGATCGTTGCGGGGTGGTGCGAGGCATAGTTGTTGATTTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACATGATCTCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATa                                                            <  3:54140/148‑1 (MQ=255)
                                                                                      cTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGc                                                          >  2:27769/1‑149 (MQ=255)
                                                                                        gatgatCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCt                                                        >  7:315015/1‑149 (MQ=255)
                                                                                        gatgatCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCt                                                        >  2:28572/1‑149 (MQ=255)
                                                                                           gatCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTgg                                                      <  5:234529/148‑1 (MQ=255)
                                                                                             tCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTAGGCGCTGATGATAGCCTGGGTa                                                   >  8:139330/1‑149 (MQ=255)
                                                                                                 gTCGGTGTTGTGCGCGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCTAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCt                                              <  5:349350/149‑1 (MQ=255)
                                                                                                  tCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTc                                             <  2:47435/149‑1 (MQ=255)
                                                                                                    ggTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGAGGATAGCCTGGGTACGCTCg                                            >  4:218134/1‑148 (MQ=255)
                                                                                                      tgtgGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTCGGCGCTGATGATAGCCTGGGTACGCTCGCaa                                         <  7:139330/149‑1 (MQ=255)
                                                                                                          gtgCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACcgc                                     >  2:47403/1‑149 (MQ=255)
                                                                                                              gAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAg                                 >  8:339705/1‑149 (MQ=255)
                                                                                                                   aTAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCa                            >  8:351975/1‑149 (MQ=255)
                                                                                                                   aTAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATAGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCGGGTTTGGCGCTGATGATAGCCTGGTTACGCTCGCAACCGCGCAGCTCCa                            >  8:76990/1‑149 (MQ=255)
                                                                                                                      gTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACg                          <  4:10929/148‑1 (MQ=255)
                                                                                                                           tGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTtc                    >  5:182765/1‑149 (MQ=255)
                                                                                                                             aTGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTtctc                  >  3:335790/1‑149 (MQ=255)
                                                                                                                              tgtgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCCTCGCAACCGCGCAGCTCCACGCCCTtctc                  <  1:28572/149‑1 (MQ=255)
                                                                                                                               gtgGTCGACGGCCTGAACCAGGCCATCGAGCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCa                >  3:147398/1‑149 (MQ=255)
                                                                                                                               gtgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCa                >  5:229853/1‑149 (MQ=255)
                                                                                                                                tgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCa                <  5:126887/148‑1 (MQ=255)
                                                                                                                                tgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCAc               >  3:297002/1‑149 (MQ=255)
                                                                                                                                 gCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTg                                                                                 >  6:76610/1‑83 (MQ=255)
                                                                                                                                 gCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTg                                                                                 <  5:76610/83‑1 (MQ=255)
                                                                                                                                  cTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCa                                  >  1:143550/1‑128 (MQ=255)
                                                                                                                                  cTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCa                                  <  2:143550/128‑1 (MQ=255)
                                                                                                                                   tCGACGGCCTGAACCAGGCCATCGAGCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACg              >  5:154148/1‑147 (MQ=255)
                                                                                                                                     gATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAgc          >  3:262012/1‑149 (MQ=255)
                                                                                                                                     gATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTAGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAgc          >  7:1789/1‑149 (MQ=255)
                                                                                                                                           cTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCggg    >  1:328811/1‑149 (MQ=255)
                                                                                                                                            tGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTgatgat                                                             >  7:271782/1‑92 (MQ=255)
                                                                                                                                            tGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTgatgat                                                             <  8:271782/92‑1 (MQ=255)
                                                                                                                                            tGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGc   >  6:231933/1‑149 (MQ=255)
                                                                                                                                              gTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACTCCCTTCTCCACGAAGCGGCGGGcc  <  4:297002/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCC  >  NC_002947/5473666‑5473955

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