Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,785,965 C→A 6.6% F323L (TTC→TTA) ‡ PP_5068 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,785,9650CA 6.6% 190.9 / 3.7 61F323L (TTC→TTA) ‡PP_5068hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base A (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

CCAACTTCATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAAC  >  NC_002947/5785817‑5786113
                                                                                                                                                    |                                                                                                                                                       
ccAACTTCATCTCAAAGCACTCGGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCTTCGCCCACTGTAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTc                                                                                                                                                         <  8:241342/149‑1 (MQ=255)
    cTTCATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGc                                                                                                                                                     >  7:269541/1‑149 (MQ=255)
        aTCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGc                                                                                                                                                  <  6:34762/148‑1 (MQ=255)
            gCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTAc                                                                                                                                             >  2:4891/1‑149 (MQ=255)
             cAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTg                                                                                                                                                <  8:269541/145‑1 (MQ=255)
             cAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTa                                                                                                                                              <  2:30499/147‑1 (MQ=255)
                 cACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGc                                                                                                                                        <  6:120726/149‑1 (MQ=255)
                      cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGCCCaac                                                                                                                                   <  4:195591/149‑3 (MQ=255)
                                        cAGCTACGCCAACCCTGTGCCCATCGCCCACTGTCACCTCAGCGCGCTGGCCCAGGAACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATgcgc                                                                                                                  <  5:293472/148‑1 (MQ=255)
                                         aGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCg                                                                                                                >  3:134368/1‑149 (MQ=255)
                                            tACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcg                                                                                                             >  8:108099/1‑149 (MQ=255)
                                            tACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcg                                                                                                             >  1:159101/1‑149 (MQ=255)
                                              cGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCg                                                                                                                                                   <  4:156693/109‑1 (MQ=255)
                                              cGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCg                                                                                                                                                   >  3:156693/1‑109 (MQ=255)
                                              cGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcgc                                                                                                            >  4:9720/1‑148 (MQ=255)
                                               gCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTTGGCCTGTTCCTGCCGCTGTAACAGGCCCATTACCTGGACCTGATGCGACGgcgca                                                                                                           >  6:244348/1‑147 (MQ=255)
                                               gCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGa                                                                                                                            <  1:126027/131‑1 (MQ=255)
                                               gCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGa                                                                                                                            >  2:126027/1‑131 (MQ=255)
                                                 cAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGCCCaa                                                                                                                                    >  5:263299/1‑120 (MQ=38)
                                                 cAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGCCCaa                                                                                                                                    <  6:263299/121‑2 (MQ=38)
                                                     cAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCt                                                                                                    <  4:134368/149‑1 (MQ=255)
                                                      aGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCt                                                                                                                               >  2:165251/1‑121 (MQ=255)
                                                      aGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCt                                                                                                                               <  1:165251/121‑1 (MQ=255)
                                                       ggTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGa                                                                                                  <  1:4891/149‑1 (MQ=255)
                                                             cATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                   <  1:155112/142‑1 (MQ=255)
                                                             cATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                   >  2:155112/1‑142 (MQ=255)
                                                               tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTa                                                                                                                                  <  2:148964/109‑1 (MQ=255)
                                                               tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTa                                                                                                                                  >  1:148964/1‑109 (MQ=255)
                                                               tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccgc                                                                                           >  7:115976/1‑148 (MQ=255)
                                                               tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccgc                                                                                           <  8:115976/148‑1 (MQ=255)
                                                                cGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCAGCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAggg                                                                                                      <  4:67698/136‑1 (MQ=255)
                                                                cGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAggg                                                                                                      >  3:67698/1‑136 (MQ=255)
                                                                 gCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGc                                                                                                                                                     >  3:94231/1‑88 (MQ=255)
                                                                 gCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGc                                                                                                                                                     <  4:94231/88‑1 (MQ=255)
                                                                  cccACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGaa                                                                                       >  1:344833/1‑149 (MQ=255)
                                                                       tGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATgc                                                                                                                    >  4:98220/1‑115 (MQ=255)
                                                                       tGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATgc                                                                                                                    <  3:98220/115‑1 (MQ=255)
                                                                            ccTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTAcc                                                                                                                                >  1:139433/1‑98 (MQ=255)
                                                                            ccTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTAcc                                                                                                                                <  2:139433/98‑1 (MQ=255)
                                                                            ccTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACa                                                                          >  3:344653/1‑149 (MQ=255)
                                                                             cTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑tgatga                                                                            >  2:217648/1‑146 (MQ=255)
                                                                             cTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑tgatga                                                                            <  1:217648/146‑1 (MQ=255)
                                                                               cAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATgg                                                                       >  5:280011/1‑149 (MQ=255)
                                                                                  cgcgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGc                                                                    <  7:108099/149‑1 (MQ=255)
                                                                                    cgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTg                                                                  <  6:280011/149‑1 (MQ=255)
                                                                                         gCCCCGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGa                                                             <  4:346855/149‑1 (MQ=255)
                                                                                           ccAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAgc                                                           <  3:9720/149‑1 (MQ=255)
                                                                                             aGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGAc                                                                                                                           >  1:222101/1‑86 (MQ=255)
                                                                                             aGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGAc                                                                                                                           <  2:222101/86‑1 (MQ=255)
                                                                                                cACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCctg                                                      <  4:344653/149‑1 (MQ=255)
                                                                                                       caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                   <  8:201756/100‑1 (MQ=255)
                                                                                                       caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                   >  7:201756/1‑100 (MQ=255)
                                                                                                       caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGt                                               <  2:159101/149‑1 (MQ=255)
                                                                                                       caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGCGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCAGACAACCGCCGAA‑‑GA‑TGATGACATGGTGCTTGTCGTGCGGCTGCTGCAGt                                               >  6:203901/1‑149 (MQ=255)
                                                                                                         ccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCTAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGtgc                                             <  5:244348/149‑1 (MQ=255)
                                                                                                            tGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGtgcatg                                          >  3:92324/1‑149 (MQ=255)
                                                                                                                 gAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAAcg                                     >  6:89036/1‑149 (MQ=255)
                                                                                                                       gAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAATAAGTGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGAgg                                  >  1:294156/1‑149 (MQ=255)
                                                                                                                             cccatcgaagCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTg                           <  1:305358/141‑1 (MQ=255)
                                                                                                                             cTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGa                         >  7:219081/1‑149 (MQ=255)
                                                                                                                                    aggtaatGGGCCTGGTACAGCGGCAGGAACAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGc                  >  7:145827/7‑149 (MQ=38)
                                                                                                                                       cGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTg               <  5:89036/149‑1 (MQ=255)
                                                                                                                                                    ccTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAAc  >  5:107847/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
CCAACTTCATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAAC  >  NC_002947/5785817‑5786113

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