Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,477,547 G→A 6.4% A140V (GCC→GTC)  PP_3088 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,477,5470GA 6.4% 184.8 / 2.3 62A140V (GCC→GTC) PP_3088hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/25);  new base A (2/2);  total (35/27)
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.01e-01

TCTCGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGAAGTGTTC  >  NC_002947/3477410‑3477695
                                                                                                                                         |                                                                                                                                                    
tctcGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGt                                                                                                                                            <  2:21309/148‑1 (MQ=255)
tctcGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTg                                                                                                                                           >  3:210716/1‑149 (MQ=255)
     gCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGt                                                                                                                                      >  8:27696/1‑149 (MQ=255)
      cGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTg                                                                                                                                     <  5:87439/149‑1 (MQ=255)
      cGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTg                                                                                                                                     <  5:249380/149‑1 (MQ=255)
           aTCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCt                                                                                                                                 >  5:117756/1‑148 (MQ=255)
                gCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCtt                                                                                                                                <  6:117756/144‑1 (MQ=255)
                   ttCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCACTGGATAATGGGCAAGCCACCGGCAATTTAGGTGAGCTTGGTATTCAGCCAg                                                                                                                         >  7:112334/1‑148 (MQ=255)
                               tGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCa                                                                                                                              >  1:258076/1‑131 (MQ=255)
                               tGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCa                                                                                                                              <  2:258076/131‑1 (MQ=255)
                                gTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCACCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGGCTTCAGCCAGGCGAGGcgcgccgc                                                                                                           >  2:285001/1‑149 (MQ=255)
                                     cATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCt                                                                                                      >  5:306603/1‑149 (MQ=255)
                                       tCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTAATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCg                                                                                                    >  1:234196/1‑149 (MQ=255)
                                           cTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCg                                                                                                                                              >  6:262237/1‑103 (MQ=255)
                                           cTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCg                                                                                                                                              <  5:262237/103‑1 (MQ=255)
                                              ggggTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCatg                                                                                              <  7:27696/148‑3 (MQ=255)
                                                 gTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGt                                                                                          >  4:283754/1‑149 (MQ=255)
                                                      gCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTc                                                                                      >  4:175999/1‑148 (MQ=255)
                                                       cAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTc                                                                                      >  2:274376/1‑147 (MQ=255)
                                                        aaGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAg                                                                                   >  2:154778/1‑149 (MQ=255)
                                                          ggCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCa                                                                                                                              <  6:263092/104‑1 (MQ=255)
                                                          ggCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCa                                                                                                                              >  5:263092/1‑104 (MQ=255)
                                                             gAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTgg                                                                              >  7:119929/1‑149 (MQ=255)
                                                              aGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTggg                                                                             <  3:283754/149‑1 (MQ=255)
                                                               ggTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTgggg                                                                            <  6:289516/149‑1 (MQ=255)
                                                                 ttCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTa                                                                          <  6:306603/149‑1 (MQ=255)
                                                                      tACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTc                                                                                                                                  <  8:131495/88‑1 (MQ=255)
                                                                      tACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTc                                                                                                                                  >  7:131495/1‑88 (MQ=255)
                                                                      tACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACa                                                                        <  7:230628/146‑1 (MQ=255)
                                                                          tGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTc                                                                 <  7:29253/149‑1 (MQ=255)
                                                                            gCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTTGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTcc                                                                >  7:275187/1‑148 (MQ=255)
                                                                            gCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGt                                                                           <  8:242587/137‑1 (MQ=255)
                                                                            gCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGt                                                                           >  7:242587/1‑137 (MQ=255)
                                                                            gCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTccc                                                               <  1:274376/149‑1 (MQ=255)
                                                                            gCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTcc                                                                >  1:157005/1‑148 (MQ=255)
                                                                                 aTTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTTGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGaa                                                           <  8:275187/148‑1 (MQ=255)
                                                                                 aTTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCg                                                                                     <  6:58372/122‑1 (MQ=255)
                                                                                 aTTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCg                                                                                     >  5:58372/1‑122 (MQ=255)
                                                                                  ttCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTc                                                                                                     <  3:219197/105‑1 (MQ=255)
                                                                                  ttCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTc                                                                                                     >  4:219197/1‑105 (MQ=255)
                                                                                     gTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCaa                                                      >  8:13712/1‑149 (MQ=255)
                                                                                           gcgcAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAg                                                  >  8:221338/1‑147 (MQ=255)
                                                                                                    gCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTccc                                                               <  6:275316/125‑1 (MQ=255)
                                                                                                    gCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTccc                                                               >  5:275316/1‑125 (MQ=255)
                                                                                                    gCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAAc                                                     >  4:97532/1‑135 (MQ=255)
                                                                                                    gCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAAc                                                     <  3:97532/135‑1 (MQ=255)
                                                                                                                 cccTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGc                           >  8:160221/1‑148 (MQ=255)
                                                                                                                 cccTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGc                           >  1:72547/1‑148 (MQ=255)
                                                                                                                  ccTGGATAATGCGCAAGCCACCGGCAATGTCGGTCAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCACGACCAGGCCGTCATCGAGGGTGGGGTACAGTGTTTCCCAGAACGGAACCACGAGCTGTTCGGTCAGCTCGACGcc                          >  4:281845/1‑148 (MQ=255)
                                                                                                                        tAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAg                                               <  3:116736/121‑1 (MQ=255)
                                                                                                                        tAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAg                                               >  4:116736/1‑121 (MQ=255)
                                                                                                                            gggCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGc               >  2:224397/1‑149 (MQ=255)
                                                                                                                             ggCAAGCCACCGGCAATGTCGGGGAGCGTGGTCTTCCGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGGACAGTGTTTCCCAGAACGCAACCAGGAGCGGTTCGGGCAGCTCGAGGCCAAAGGTGAGGcc              >  8:69338/1‑149 (MQ=255)
                                                                                                                                cgctcACCGACATTGCCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTc           >  2:35557/5‑149 (MQ=255)
                                                                                                                                 gctcACCGACATTGCCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTc           >  5:173913/4‑148 (MQ=255)
                                                                                                                                  ctcACCGACATTGCCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTGTGTCAGCTCGAGGCCAAAGGTGATGCCTTCGa         <  1:35557/147‑1 (MQ=255)
                                                                                                                                   ccACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGaa        <  8:112334/149‑1 (MQ=255)
                                                                                                                                   ccACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGaa        <  2:234196/149‑1 (MQ=255)
                                                                                                                                     aCCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGAAgt      >  1:233830/1‑149 (MQ=255)
                                                                                                                                     aCCGACATTGCCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGaa        <  6:173913/147‑1 (MQ=255)
                                                                                                                                      ccGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGaa        <  6:82038/146‑1 (MQ=255)
                                                                                                                                      ccGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGaa        >  5:82038/1‑146 (MQ=255)
                                                                                                                                         gCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGAAGTGTTc  <  7:13712/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
TCTCGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGAAGTGTTC  >  NC_002947/3477410‑3477695

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