Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,051,604 T→C 6.6% K106E (AAA→GAA)  PP_0909 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,051,6040TC 6.6% 173.7 / 2.2 60K106E (AAA→GAA) PP_0909hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base C (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GATGTCCGGCTTGCTGTCGAACAGGATGTACTCATAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTCCTTGCCTTCGCAA  >  NC_002947/1051464‑1051746
                                                                                                                                             |                                                                                                                                              
gATGTCCGGCTTGCTGTCGAACAGGATGTACTCATAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAg                                                                                                                                        >  4:254522/1‑149 (MQ=255)
      cGGCTTGCTGTCGAACAGGATGTACTCATAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAg                                                                                                                                   >  4:223135/1‑148 (MQ=255)
      cGGCTTGCTGTCGAACAGGATGTACTCATAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAg                                                                                                                                   >  2:270292/1‑148 (MQ=255)
             cTGTCGAACAGGATGTACTCATAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTAc                                                                                                                           >  6:181647/1‑149 (MQ=255)
              tGTCGAACAGGATGTACTCATAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTAcc                                                                                                                          >  2:178939/1‑149 (MQ=255)
              tGTCGAACAGGATGTACTCATAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTATCAGGCCAGAATTGTAcc                                                                                                                          >  4:345847/1‑149 (MQ=255)
                       ggATGTACTCATAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGGCAGGCCAGAATTGTACCTGCCAGgcg                                                                                                                 <  8:27488/149‑1 (MQ=255)
                        gATGTACTCATAGGCCCACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGgcgc                                                                                                                >  4:172327/1‑149 (MQ=255)
                            tACTCATAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGt                                                                                                             <  5:155754/148‑1 (MQ=255)
                                 aTAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCtt                                                                                                       >  3:118536/1‑149 (MQ=255)
                                  tAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTc                                                                                                      <  4:118536/149‑1 (MQ=255)
                                   aGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCtt                                                                                                       >  6:66168/1‑147 (MQ=255)
                                    ggCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCa                                                                                                     <  3:107080/148‑1 (MQ=255)
                                     gCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTGCCACTACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAg                                                                                                    <  1:111879/149‑1 (MQ=255)
                                        gACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCaa                                                                                                >  2:217063/1‑149 (MQ=255)
                                            aGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGCCCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAAGTGTACCTGCCAGGCGCGGTTGCCTTCAGGCCACggg                                                                                            >  4:28874/1‑149 (MQ=255)
                                            aGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACggg                                                                                            <  1:270292/149‑1 (MQ=255)
                                               cACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAAt                                                                                                                                <  3:80071/110‑1 (MQ=255)
                                               cACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAAt                                                                                                                                >  4:80071/1‑110 (MQ=255)
                                                  gTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGACGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTATTCTTGTAAGGCCAGAATTGTAACTGCCAGAAGAGGTTTCCTTCAGCAAAAGGTCAGAc                                                                                       >  5:94506/1‑148 (MQ=255)
                                                      cACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCTACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTa                                                                                  <  2:189286/149‑1 (MQ=255)
                                                        attcGCTGGACTTGCCCAACGACTCAGCTTGAACCGGCTTGTCGCAGTTGATGAGTTG‑CTCGACCTGACGCCCGACCATGTTTTTATTGTCTGTACAGAATTGTACCTGAAAGGCGCGGTTGCATTAAGCCAACGGTACGACCATTatc                                                                                <  6:94506/146‑3 (MQ=255)
                                                             cTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTCCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTg                                                                           >  2:81792/1‑149 (MQ=255)
                                                             cTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTCCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTg                                                                           >  2:81794/1‑149 (MQ=255)
                                                              tGGACTTGCCCAACGCCTCAGCTTCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGc                                                                          >  7:54952/1‑149 (MQ=255)
                                                                gACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTCCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGc                                                                          <  1:81794/147‑1 (MQ=255)
                                                                gACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTCCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGc                                                                          <  1:81792/147‑1 (MQ=255)
                                                                   ttGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTc                                                                     >  2:264831/1‑149 (MQ=255)
                                                                        cAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGccc                                                                <  5:181647/149‑1 (MQ=255)
                                                                          aCGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAg                                                              <  3:254522/149‑1 (MQ=255)
                                                                              cTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGtt                                                                                                            <  5:245459/99‑1 (MQ=255)
                                                                              cTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGtt                                                                                                            >  6:245459/1‑99 (MQ=255)
                                                                               tCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCttt                                                         >  2:44610/1‑149 (MQ=255)
                                                                                 aGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTATCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCttttt                                                       <  3:345847/149‑1 (MQ=255)
                                                                                       aaCCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCAt                                                 <  3:223135/149‑1 (MQ=255)
                                                                                        aCCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCAtt                                                >  1:84116/1‑149 (MQ=255)
                                                                                        aCCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCAtt                                                >  1:88947/1‑149 (MQ=255)
                                                                                                          aCCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGc                              >  2:102742/1‑149 (MQ=255)
                                                                                                           ccAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCt                             <  1:44610/149‑1 (MQ=255)
                                                                                                           ccAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCt                             <  2:88947/149‑1 (MQ=255)
                                                                                                           ccAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCt                             <  2:84116/149‑1 (MQ=255)
                                                                                                             aGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCttt                           >  3:105038/1‑149 (MQ=255)
                                                                                                              gTTG‑CTTGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTg                          >  7:116123/1‑149 (MQ=255)
                                                                                                               ttG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCAGTGAGGAAGTCGGCGCGGGCTTTGt                         >  6:253296/1‑149 (MQ=255)
                                                                                                                    tCGACCTGACGCCCGACCAGGTTTTTCATGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAg                     >  1:196335/1‑149 (MQ=255)
                                                                                                                     cGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAgg                    <  8:54952/149‑1 (MQ=255)
                                                                                                                      gACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAgg                    <  1:178939/148‑1 (MQ=255)
                                                                                                                            aCGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTcc              <  8:116123/148‑1 (MQ=255)
                                                                                                                             cGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGGTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTCCtt            <  2:263673/149‑1 (MQ=255)
                                                                                                                              gCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTCCTTg           >  4:33719/1‑149 (MQ=255)
                                                                                                                              gCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTCCTTg           <  1:264831/149‑1 (MQ=255)
                                                                                                                               cccGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCgg                                >  5:365241/1‑127 (MQ=255)
                                                                                                                               cccGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCgg                                <  6:365241/127‑1 (MQ=255)
                                                                                                                               cccGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTCCTTGc          <  3:172327/149‑1 (MQ=255)
                                                                                                                                 cGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTCCTTGCCt        <  4:105038/149‑1 (MQ=255)
                                                                                                                                   aCCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGGTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTCCTTGCCTTc      <  3:28874/149‑1 (MQ=255)
                                                                                                                                    ccAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTCCTTGCCTTCg     >  8:115072/1‑149 (MQ=255)
                                                                                                                                       ggTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTCCTTGCCTTCGCaa  <  2:196335/149‑1 (MQ=255)
                                                                                                                                        gTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTTGGCGCGGGCt                             >  7:7228/1‑121 (MQ=255)
                                                                                                                                        gTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTTGGCGCGGGCt                             <  8:7228/121‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GATGTCCGGCTTGCTGTCGAACAGGATGTACTCATAGGCCGACAAGCCACGTACCACTACGCTGGACTTGCCCAACGCCTCAGCATCAACCGGCTTGTCGCCGTTGACCAGTTG‑CTCGACCTGACGCCCGACCAGGTTTTTCTTGTCAGGCCAGAATTGTACCTGCCAGGCGCGGTTGCCTTCAGCCAACGGGCCGACCAGTAGCGGTTGCAGCTCGGCCCAGGCTTTTTGTGCATTGAGGAAGTCGGCGCGGGCTTTGTCCAGGCTTTCCTTGCCTTCGCAA  >  NC_002947/1051464‑1051746

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