Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,248,904 C→G 6.6% G375R (GGC→CGC) ‡ PP_3724 ← acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,248,9040CG 6.6% 192.0 / 4.0 61G375R (GGC→CGC) ‡PP_3724acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base G (2/2);  total (29/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 = 3.63e-01

GCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTGCTCGA  >  NC_002947/4248760‑4249047
                                                                                                                                                |                                                                                                                                               
gcaACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCcgcg                                                                                                                                             <  6:185003/149‑1 (MQ=255)
 caACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCgcgg                                                                                                                                            <  4:222874/149‑1 (MQ=255)
     aGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGGGCgcggcgat                                                                                                                                        <  5:363534/149‑3 (MQ=255)
           gggAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAg                                                                                                                                  >  6:196386/1‑149 (MQ=255)
            ggAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAgg                                                                                                                                 <  7:240610/149‑1 (MQ=255)
            ggAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAgg                                                                                                                                 <  1:76213/149‑1 (MQ=255)
              aaaGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGgt                                                                                                                               <  6:105182/149‑1 (MQ=255)
                 gcggatgTGACGAAGCCGATCTTGCGGTTGCGTTCAACGGCATCGCAGTACCTCCGCAGGAGCAGGTCGCATTCCTTGCGTTCGGCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCTCCCGCGGCGCCGCTCAGGGtgtt                                                                                                                            <  6:5832/143‑1 (MQ=255)
                 gCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGtgtt                                                                                                                            >  4:183800/1‑149 (MQ=255)
                    gAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTc                                                                                                                         <  5:32581/149‑1 (MQ=255)
                          aCGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCAt                                                                                                                   >  4:204899/1‑149 (MQ=255)
                            gAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGGGCgcggcga                                                                                                                                         >  3:52001/1‑124 (MQ=255)
                            gAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGGGCgcggcga                                                                                                                                         <  4:52001/125‑2 (MQ=255)
                            gAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGc                                                                                                                 >  5:3715/1‑149 (MQ=255)
                               gCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGtt                                                                                                              >  4:325843/1‑149 (MQ=255)
                                 cGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTg                                                                                                             >  3:223789/1‑148 (MQ=255)
                                        gCCGTGGCGTTCAACGGCATCGCAGTCCAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGggcg                                                                                                      <  1:263412/148‑1 (MQ=255)
                                        gCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGggcg                                                                                                      <  8:124279/148‑1 (MQ=255)
                                          cGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCt                                                                                                   >  2:149826/1‑149 (MQ=255)
                                            tGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCt                                                                                                   >  2:319978/1‑147 (MQ=255)
                                            tGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCt                                                                                                   <  1:319978/147‑1 (MQ=255)
                                              gCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCAt                                                                                               <  5:196386/149‑1 (MQ=255)
                                                 ttCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGc                                                                                            >  8:182048/1‑149 (MQ=255)
                                                  tCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGc                                                                                            <  2:309670/148‑1 (MQ=255)
                                                  tCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGc                                                                                            >  1:309670/1‑148 (MQ=255)
                                                      cGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCCGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCGTCATCGccc                                                                                          >  8:159958/1‑146 (MQ=255)
                                                        gCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGAt                                                                                     <  3:325843/149‑1 (MQ=255)
                                                          aTCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGtgtt                                                                                                                            <  6:12057/108‑1 (MQ=255)
                                                          aTCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGtgtt                                                                                                                            >  5:12057/1‑108 (MQ=255)
                                                          aTCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTAGCAGGGCGGCTTCATCGCCCCAGATgg                                                                                   >  7:14159/1‑149 (MQ=255)
                                                                gTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAgtgg                                                                             <  4:223789/149‑1 (MQ=255)
                                                                  aCAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAgg                                                                                                                                 <  2:113136/95‑1 (MQ=255)
                                                                  aCAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAgg                                                                                                                                 >  1:113136/1‑95 (MQ=255)
                                                                   cAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAgtggtg                                                                           >  1:133602/1‑148 (MQ=255)
                                                                    aGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGaa                                                                         <  1:149826/149‑1 (MQ=255)
                                                                           gcagcagGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGt                                                                  <  7:182048/149‑1 (MQ=255)
                                                                             agcagGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTgg                                                                <  3:204899/149‑1 (MQ=255)
                                                                             agcagGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTgg                                                                >  7:68685/1‑149 (MQ=255)
                                                                                     gCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGAcc                                                         <  8:14159/148‑1 (MQ=255)
                                                                                      cATTCCTTGCGTTCGCCGGTAAACGCCACCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAg                                                       >  1:11296/1‑149 (MQ=255)
                                                                                           cTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTc                                                                      >  4:255195/1‑129 (MQ=255)
                                                                                            ttGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTgg                                                  >  6:315051/1‑148 (MQ=255)
                                                                                                gTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTc                                                                      <  6:306190/124‑1 (MQ=255)
                                                                                                   cGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGccc                                          >  5:154969/1‑149 (MQ=255)
                                                                                                    gCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCg                                         <  5:315051/149‑1 (MQ=255)
                                                                                                      cGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGc                                       <  8:47763/149‑1 (MQ=255)
                                                                                                       ggTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAg                                             <  6:192015/142‑1 (MQ=255)
                                                                                                       ggTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTGGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAg                                             >  5:192015/1‑142 (MQ=255)
                                                                                                           aaCGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATa                                              >  2:350647/1‑137 (MQ=255)
                                                                                                           aaCGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATa                                              <  1:350647/137‑1 (MQ=255)
                                                                                                              gccgccGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTc                               >  8:294348/1‑149 (MQ=255)
                                                                                                                      gTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTg                       <  7:294348/149‑1 (MQ=255)
                                                                                                                         gTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTgccg                    <  2:133602/149‑1 (MQ=255)
                                                                                                                                accagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCAc                <  7:133112/146‑1 (MQ=255)
                                                                                                                                accagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCAc                >  8:133112/1‑146 (MQ=255)
                                                                                                                                  cagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGcc           >  6:53211/1‑149 (MQ=255)
                                                                                                                                  cagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGcc           <  6:154969/149‑1 (MQ=255)
                                                                                                                                       gcgccgcGGGCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTGc      >  3:23802/2‑149 (MQ=255)
                                                                                                                                          ccgcGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTGCTCg   >  4:197170/1‑149 (MQ=255)
                                                                                                                                          ccgcGCCCGCGGCGCCACTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTGCTCg   >  8:203441/1‑149 (MQ=255)
                                                                                                                                           cgcgCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTGCTCGa  <  8:68685/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTGCTCGA  >  NC_002947/4248760‑4249047

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: ‑