Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,691,282 G→T 8.3% L194I (CTC→ATC)  PP_3257 ← putative Methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,691,2820GT 8.3% 182.6 / 5.5 60L194I (CTC→ATC) PP_3257putative Methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (25/30);  new base T (3/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

AGCGGTCCAGTACCACCAGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAGGGTGC  >  NC_002947/3691133‑3691424
                                                                                                                                                     |                                                                                                                                              
aGCGGTCCAGTACCACCAGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCgag                                                                                                                                                >  2:452974/1‑150 (MQ=255)
   ggTCCAGTACCACCAGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATc                                                                                                                                             >  2:832307/1‑150 (MQ=255)
           accaccAGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCGGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGATCTcgacgaga                                                                                                                                      >  1:148855/1‑147 (MQ=255)
           accaccAGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGATCTcgacgaga                                                                                                                                      <  2:148855/149‑3 (MQ=255)
             caccaGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCtt                                                                                                                                   <  2:250338/150‑1 (MQ=255)
               ccaGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTtgg                                                                                                                                 <  1:135096/150‑1 (MQ=255)
               ccaGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTtgg                                                                                                                                 >  1:366544/1‑150 (MQ=255)
                caGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGAAGACGCTtggc                                                                                                                                >  1:129198/1‑150 (MQ=255)
                 aGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTtggct                                                                                                                               >  2:1106536/1‑150 (MQ=255)
                     cTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCa                                                                                                                           <  1:452974/150‑1 (MQ=255)
                      ttCCCCGTGCAGACCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCAGCGGGGATGGCCATGTTCGGGGTGCGTCGTTATGTGGATGAATGACAACGGCCATCATGGCCAGGTCTCGTCGAGATCGACGACGCTTGGCTGGCAt                                                                                                                          <  1:1178508/150‑1 (MQ=255)
                      ttCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCAt                                                                                                                          >  2:977346/1‑150 (MQ=255)
                        ccACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGa                                                                                                                        <  2:546924/150‑1 (MQ=255)
                         cACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGAc                                                                                                                       >  2:463963/1‑150 (MQ=255)
                         cACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGATTCGTTGTGTGCCGGGATGACAACGACAGTCATGGTCACGGCTCGTCGAGATCGACGACGCTTGGCTGGCATGAc                                                                                                                       >  2:738358/1‑150 (MQ=255)
                          aCGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTtggct                                                                                                                               <  2:898088/141‑1 (MQ=255)
                          aCGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTtggct                                                                                                                               >  1:898088/1‑141 (MQ=255)
                           cGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCa                                                                                                                     <  2:722969/150‑1 (MQ=255)
                                  tCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTcc                                                                                                              <  2:1059844/150‑1 (MQ=255)
                                   cccGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCa                                                                                                             <  1:302645/150‑1 (MQ=255)
                                    ccGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAg                                                                                                            >  1:1203625/1‑150 (MQ=255)
                                          gctgcGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCaaaaa                                                                                                      >  2:128517/1‑150 (MQ=255)
                                                ccTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCCTATCGCTATTCCAGTGGTGGGTGGCATGACCAGTATTCCAGCAAAAAGTCCGc                                                                                                <  1:1114386/150‑1 (MQ=255)
                                                     cGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGcctcct                                                                                           <  2:810318/150‑1 (MQ=255)
                                                             caccaGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGa                                                                                   <  1:832307/150‑1 (MQ=255)
                                                               ccaGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGaaa                                                                                 <  1:3185/150‑1 (MQ=255)
                                                                  gCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAActg                                                                              >  1:199891/1‑150 (MQ=255)
                                                                   cACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAActgc                                                                             <  2:129198/150‑1 (MQ=255)
                                                                    aCCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAActgct                                                                            <  1:128517/150‑1 (MQ=255)
                                                                    aCCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAActgct                                                                            >  2:1110002/1‑150 (MQ=255)
                                                                    aCCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAActgct                                                                            <  2:366544/150‑1 (MQ=255)
                                                                           gcAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCaa                                                                     <  1:977346/150‑1 (MQ=255)
                                                                               gCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGAt                                                                                                                                              >  2:711231/1‑73 (MQ=255)
                                                                               gCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGAt                                                                                                                                              <  1:711231/73‑1 (MQ=255)
                                                                                 aTCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCa                                                                                                             <  2:559824/104‑1 (MQ=255)
                                                                                 aTCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCa                                                                                                             >  1:559824/1‑104 (MQ=255)
                                                                                       cAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCg                                                         <  1:1068095/150‑1 (MQ=255)
                                                                                          ggcggcGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGt                                                      <  2:199891/150‑1 (MQ=255)
                                                                                              gcggGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAAt                                                  <  1:1084166/150‑1 (MQ=255)
                                                                                                  gCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGc                                              <  1:738358/150‑1 (MQ=255)
                                                                                                   cTTCGTTATGTGCAGGGATGACACCGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCg                                             >  1:1158711/1‑150 (MQ=255)
                                                                                                   cTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCg                                             >  1:280993/1‑150 (MQ=255)
                                                                                                             tgCAGGGATGACAACGGCAATCATGGTCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTcg                                   >  2:488875/1‑150 (MQ=255)
                                                                                                             tgCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTcg                                   >  2:146544/1‑150 (MQ=255)
                                                                                                             tgCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTcg                                   >  2:474493/1‑150 (MQ=255)
                                                                                                                 gggATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGAc                               <  2:1203625/150‑1 (MQ=255)
                                                                                                                 gggATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGAc                               >  2:511565/1‑150 (MQ=255)
                                                                                                                  ggATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGACCGCCTCCTCAt                                                                                        >  1:704586/1‑92 (MQ=255)
                                                                                                                  ggATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGACCGCCTCCTCAt                                                                                        <  2:704586/92‑1 (MQ=255)
                                                                                                                         aaCGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGc                       >  1:747569/1‑150 (MQ=255)
                                                                                                                            ggCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCggg                    <  1:704794/150‑1 (MQ=255)
                                                                                                                            ggCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCggg                    >  2:571820/1‑150 (MQ=255)
                                                                                                                               aaTCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCaa                 <  2:747569/150‑1 (MQ=255)
                                                                                                                                  cATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAAccc              >  2:929593/1‑150 (MQ=255)
                                                                                                                                   aTGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCcg             <  2:1158711/150‑1 (MQ=255)
                                                                                                                                          cGTCTCGTCGAGACCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAgg      <  1:204999/150‑1 (MQ=255)
                                                                                                                                              tcgtcgATCTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAggg     >  1:126448/1‑147 (MQ=255)
                                                                                                                                              tcgtcgATCTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAggg     <  2:126448/147‑1 (MQ=255)
                                                                                                                                              tcgtcgATCTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAGGGTGc  >  1:424785/1‑150 (MQ=255)
                                                                                                                                              tcgtcgAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAAt                                                  <  2:24632/102‑1 (MQ=255)
                                                                                                                                              tcgtcgAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAAt                                                  >  1:24632/1‑102 (MQ=255)
                                                                                                                                                   gagaTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTc                                    >  2:295776/1‑111 (MQ=255)
                                                                                                                                                   gagaTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTc                                    <  1:295776/111‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
AGCGGTCCAGTACCACCAGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAGGGTGC  >  NC_002947/3691133‑3691424

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