Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 779,919 A→C 6.2% G184G (GGT→GGG)  PP_0670 ← bile acid/Na+ symporter family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947779,9190AC 6.2% 201.2 / 3.6 64G184G (GGT→GGG) PP_0670bile acid/Na+ symporter family transporter
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base C (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CACCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAG  >  NC_002947/779778‑780065
                                                                                                                                             |                                                                                                                                                  
cACCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGgcgc                                                                                                                                             <  2:112092/149‑1 (MQ=255)
 aCCGGCAGCAAACGGTCGAAATCACGCAGTACTACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCt                                                                                                                                            >  8:79178/1‑149 (MQ=255)
 aCCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGTGCCCGgcgc                                                                                                                                             <  7:180695/148‑1 (MQ=255)
  ccGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTg                                                                                                                                           <  5:107045/149‑1 (MQ=255)
   cGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTgg                                                                                                                                          <  2:202823/149‑1 (MQ=255)
     gcagcaAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGt                                                                                                                                         >  3:33350/1‑148 (MQ=255)
                tCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGt                                                                                                                                         <  3:51811/137‑1 (MQ=255)
                tCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGt                                                                                                                                         >  4:51811/1‑137 (MQ=255)
                tCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAg                                                                                                                             >  6:154144/1‑149 (MQ=255)
                  gcaaTCACGCAGTACCACGTCAATCTGCGAAGCGCTCACCACCATCAGCACCAAGGCCCTGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAg                                                                                                                             <  6:22534/145‑1 (MQ=255)
                  gAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAggg                                                                                                                           >  6:70507/1‑149 (MQ=255)
                  gAAATCACGCACTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAgg                                                                                                                            <  4:208742/148‑1 (MQ=255)
                   aaaTCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGc                                                                                                                          >  5:329251/1‑149 (MQ=255)
                     aTCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCa                                                                                                                        >  4:183152/1‑149 (MQ=255)
                      tCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGACGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTCGTGAGCACCGCGAGGGCCAt                                                                                                                       >  4:180091/1‑149 (MQ=255)
                          gCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGc                                                                                                                   <  1:326092/149‑1 (MQ=255)
                             gTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCATGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTTAGCACCGCGAGGGCCATCGGCAac                                                                                                                <  2:181551/149‑1 (MQ=255)
                              tACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAaca                                                                                                               <  7:191123/149‑1 (MQ=255)
                              tACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAaca                                                                                                               >  2:201461/1‑149 (MQ=255)
                                    gTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATc                                                                                                          >  3:202698/1‑148 (MQ=255)
                                         tCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGa                                                                                                                              <  5:267755/123‑1 (MQ=255)
                                         tCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGa                                                                                                                              >  6:267755/1‑123 (MQ=255)
                                            gCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCAc                                                                                                  >  3:77665/1‑148 (MQ=255)
                                             cGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTGCCAGGGCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCa                                                                                                   <  6:293047/146‑1 (MQ=255)
                                                    aTCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTc                                                                                         <  4:33350/149‑1 (MQ=255)
                                                    aTCACCACCATCAGCACC‑AGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTc                                                                                         >  4:267151/1‑148 (MQ=255)
                                                     tCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTcc                                                                                        <  6:329251/149‑1 (MQ=255)
                                                       accaccATCAGCACC‑AGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAg                                                                                      <  3:267151/148‑1 (MQ=255)
                                                         caccaTCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGa                                                                                    <  5:70507/149‑1 (MQ=255)
                                                             aTCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATc                                                                                                          >  6:290425/1‑123 (MQ=255)
                                                             aTCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATc                                                                                                          <  5:290425/123‑1 (MQ=255)
                                                              tCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACggg                                                                               >  6:118815/1‑149 (MQ=255)
                                                                          ccATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGACGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGAtt                                                                   >  4:135150/1‑149 (MQ=255)
                                                                           cATGGCCGGCTCCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTg                                                                  >  1:140858/1‑149 (MQ=255)
                                                                            aTGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGc                                                                 >  8:212306/1‑149 (MQ=255)
                                                                            aTGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGc                                                                 >  1:156097/1‑149 (MQ=255)
                                                                                       cggcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGAtt                                                                   <  8:8704/136‑1 (MQ=255)
                                                                                       cggcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGAtt                                                                   >  7:8704/1‑136 (MQ=255)
                                                                                              cACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCAcc                                               >  7:32437/1‑149 (MQ=255)
                                                                                              cACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCAcc                                               >  3:199690/1‑149 (MQ=255)
                                                                                                cGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCa                                             <  5:118815/149‑1 (MQ=255)
                                                                                                 gCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATgg                                                                        <  8:165545/121‑1 (MQ=255)
                                                                                                 gCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATgg                                                                        >  7:165545/1‑121 (MQ=255)
                                                                                                 gCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAg                                            >  3:230666/1‑149 (MQ=255)
                                                                                                   ccAGGCGTCGTTCCTGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCa                                          >  5:231329/1‑149 (MQ=255)
                                                                                                     aGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGATGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATc                                        <  3:135150/149‑1 (MQ=255)
                                                                                                      ggcggcGGTCCAGCTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGACCGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCa                                       <  8:124969/149‑1 (MQ=255)
                                                                                                      ggCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCa                                       <  4:77665/149‑1 (MQ=255)
                                                                                                        cgtcgtTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGTCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTTATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAgg                                     <  1:201461/149‑1 (MQ=255)
                                                                                                        cgtcgtTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCATCAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAgg                                     <  6:10818/149‑1 (MQ=255)
                                                                                                                cAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGt                              <  7:212306/148‑1 (MQ=255)
                                                                                                                   gTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTacac                          <  3:180091/149‑1 (MQ=255)
                                                                                                                    tCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACAcc                         <  4:230666/149‑1 (MQ=255)
                                                                                                                    tCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACAcc                         <  4:202698/149‑1 (MQ=255)
                                                                                                                      gcgACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCgg                       <  3:183152/149‑1 (MQ=255)
                                                                                                                        gACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCgg                       >  2:68650/1‑147 (MQ=255)
                                                                                                                                accaGCGCCGGTGCCCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCAc              >  3:75611/4‑148 (MQ=255)
                                                                                                                                  caGCGCCGGTGCCCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAg           <  4:75611/148‑1 (MQ=255)
                                                                                                                                  caGCGCCGGTGCCCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAg           >  4:71835/2‑149 (MQ=255)
                                                                                                                                      gcCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTgc       <  1:278391/149‑1 (MQ=255)
                                                                                                                                        cGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAgg                                     <  6:287899/117‑1 (MQ=255)
                                                                                                                                        cGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAgg                                     >  5:287899/1‑117 (MQ=255)
                                                                                                                                        cGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTgcgc     <  4:199690/149‑1 (MQ=255)
                                                                                                                                            cACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCTGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAg  >  7:3886/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CACCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTGCGCCAG  >  NC_002947/779778‑780065

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