Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 790,960 T→A 9.8% Q46L (CAG→CTG)  PP_0677 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947790,9600TA 9.8% 238.5 / 13.4 81Q46L (CAG→CTG) PP_0677lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base A (3/5);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

AGCGCTCGTGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACAAAGCCC  >  NC_002947/790814‑791109
                                                                                                                                                  |                                                                                                                                                     
aGCGCTCGTGGTCTTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTggc                                                                                                                                                    <  1:293703/150‑1 (MQ=255)
      cGTGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTTACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccgg                                                                                                                                              <  1:1347844/150‑3 (MQ=255)
      cGTGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccgg                                                                                                                                              <  1:388844/150‑3 (MQ=255)
        tGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccgggc                                                                                                                                            <  2:289865/150‑5 (MQ=255)
             ttGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCa                                                                                                                                       <  2:1082028/150‑1 (MQ=255)
                gcgcTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGcc                                                                                                                                    >  2:454560/1‑150 (MQ=255)
                    tGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTcc                                                                                                                                <  1:800234/150‑1 (MQ=255)
                         ttcGCTGTTCTGCATGGCGTGAAACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGaa                                                                                                                           >  1:1343780/1‑150 (MQ=255)
                               gTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGTAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCTCCCCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATc                                                                                                                     <  2:507346/150‑1 (MQ=255)
                                ttCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATcc                                                                                                                    >  2:1116800/1‑150 (MQ=255)
                                   tGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccggg                                                                                                                                             >  2:1184743/1‑119 (MQ=255)
                                   tGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCAGGCGGccggg                                                                                                                                             <  1:1184743/122‑4 (MQ=255)
                                     cATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATa                                                                                                               <  2:227161/150‑1 (MQ=255)
                                         gCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCa                                                                                                           <  1:42163/150‑1 (MQ=255)
                                         gCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCa                                                                                                           <  1:250723/150‑1 (MQ=255)
                                            ggACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCagga                                                                                                        <  2:756189/149‑1 (MQ=255)
                                              acacTAGCGGAACCCATCGTCCCAGCCGCTTTCGTACCGTCGCTCGTGCAGGTAGCGCCGTACGTCCTTGCGGAAGTTGCCCACCATCGGCACGGCCGCATGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCaggagg                                                                                                      <  2:541599/150‑1 (MQ=255)
                                                acTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTTCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGt                                                                                                    >  2:100500/1‑150 (MQ=255)
                                                            cATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGa                                                                                        >  1:847136/1‑150 (MQ=255)
                                                             aTCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGAc                                                                                       >  1:1099559/1‑150 (MQ=255)
                                                             aTCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGAc                                                                                       >  1:1164341/1‑150 (MQ=255)
                                                              tCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTggcggc                                                                                                                                                 <  1:1060691/91‑1 (MQ=255)
                                                              tCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTggcggc                                                                                                                                                 >  2:1060691/1‑91 (MQ=255)
                                                                 tCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGaa                                                                                                                           <  2:172037/110‑1 (MQ=255)
                                                                 tCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGaa                                                                                                                           >  1:172037/1‑110 (MQ=255)
                                                                    cAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCGTGGAAGCCATCCGCATAGGCagga                                                                                                        >  2:956078/1‑126 (MQ=255)
                                                                    cAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCGTGGAAGCCATCCGCATAGGCagga                                                                                                        >  2:955991/1‑126 (MQ=255)
                                                                    cAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCagga                                                                                                        <  1:956078/126‑1 (MQ=255)
                                                                    cAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCagga                                                                                                        <  1:955991/126‑1 (MQ=255)
                                                                    cAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTg                                                                                             >  2:1137885/1‑137 (MQ=255)
                                                                    cAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTg                                                                                             <  1:1137885/137‑1 (MQ=255)
                                                                     aGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCt                                                                               >  1:289843/1‑150 (MQ=255)
                                                                      gCCGGTTGCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTc                                                                              >  1:1128277/1‑150 (MQ=255)
                                                                       ccGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACGGTGGTTGACCATCTGCTCg                                                                             >  2:934751/1‑150 (MQ=255)
                                                                               gTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCgtgggttgc                                                                     >  2:180202/1‑147 (MQ=255)
                                                                                 aCTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCAt                                                                                    <  2:29670/133‑1 (MQ=255)
                                                                                 aCTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCAt                                                                                    >  1:29670/1‑133 (MQ=255)
                                                                                  cTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTc                                                                  <  2:951771/150‑1 (MQ=255)
                                                                                   tGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCaggagg                                                                                                      >  1:1361636/1‑113 (MQ=255)
                                                                                   tGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCaggagg                                                                                                      <  2:1361636/113‑1 (MQ=255)
                                                                                      cGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCaa                                                              <  2:851321/150‑1 (MQ=255)
                                                                                       gCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAAc                                                             <  1:707099/150‑1 (MQ=255)
                                                                                            tGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGAc                                                        <  2:1128277/150‑1 (MQ=255)
                                                                                             gCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACa                                                       <  2:847136/150‑1 (MQ=255)
                                                                                              cAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCt                                                                               <  1:211941/125‑1 (MQ=255)
                                                                                              cAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCt                                                                               >  2:211941/1‑125 (MQ=255)
                                                                                                   aGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGc                                                     <  2:735855/146‑1 (MQ=255)
                                                                                                   aGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGc                                                     >  1:735855/1‑146 (MQ=255)
                                                                                                          tACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAAt                                          <  1:379503/150‑1 (MQ=255)
                                                                                                             gTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGca                                       >  2:1162984/1‑150 (MQ=255)
                                                                                                               tCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCag                                                                                                          >  1:963589/1‑81 (MQ=255)
                                                                                                               tCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCag                                                                                                          <  2:963589/81‑1 (MQ=255)
                                                                                                                 ttGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAg                                                                                                              <  1:415726/75‑1 (MQ=255)
                                                                                                                 ttGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAg                                                                                                              >  2:415726/1‑75 (MQ=255)
                                                                                                                   gCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATg                                         <  2:543550/142‑1 (MQ=255)
                                                                                                                   gCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATg                                         >  1:543550/1‑142 (MQ=255)
                                                                                                                    cGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCg                                >  2:607794/1‑150 (MQ=255)
                                                                                                                       aaGTTGCCCCCCATCAGCCCCTCCGCCGGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATGCGCATAGGCAGGTGGGTAACCGTTGTTGGCCATCTGCTCGTCGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAAc                             <  2:391000/150‑1 (MQ=255)
                                                                                                                       aaGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAAc                             <  2:679941/150‑1 (MQ=255)
                                                                                                                       aaGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAAc                             <  2:919031/150‑1 (MQ=255)
                                                                                                                            gCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTa                        <  2:1343780/150‑1 (MQ=255)
                                                                                                                              ccaccaTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACa                                                       >  1:48602/1‑117 (MQ=255)
                                                                                                                              ccaccaTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACa                                                       <  2:48602/117‑1 (MQ=255)
                                                                                                                               caccaTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACcg                     >  2:751791/1‑150 (MQ=255)
                                                                                                                                accaTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACcgc                    >  2:1219365/1‑150 (MQ=255)
                                                                                                                                  caTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGTAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGa                  <  1:1162984/150‑1 (MQ=255)
                                                                                                                                  caTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGa                  >  2:713211/1‑150 (MQ=255)
                                                                                                                                   aTCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGAc                 >  1:837241/1‑150 (MQ=255)
                                                                                                                                    tCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACa                <  2:289843/150‑1 (MQ=255)
                                                                                                                                      gcagtggCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATa              >  2:1190806/6‑150 (MQ=255)
                                                                                                                                       gCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGAc                 <  2:98827/146‑1 (MQ=255)
                                                                                                                                       gCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGAc                 >  1:98827/1‑146 (MQ=255)
                                                                                                                                       gCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAg             >  2:856408/1‑150 (MQ=255)
                                                                                                                                        agtggCCGCCAGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCacac                                    <  2:822589/123‑1 (MQ=255)
                                                                                                                                        agtggCCGCCAGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCacac                                    >  1:822589/4‑126 (MQ=255)
                                                                                                                                         gtggCCGCCAGGCGGCCACTGCTGCAGCCGTCATGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGcc           >  1:1341676/3‑150 (MQ=255)
                                                                                                                                            gccgccTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGGAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCcaca        >  1:1100446/1‑150 (MQ=255)
                                                                                                                                             ccgccTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACaa       >  1:843024/1‑150 (MQ=255)
                                                                                                                                              cgccTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACaaa      <  1:1116800/150‑1 (MQ=255)
                                                                                                                                                 cTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACAAAGcc   <  1:1190806/150‑1 (MQ=255)
                                                                                                                                                 cTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACAAAGcc   <  2:169528/150‑1 (MQ=255)
                                                                                                                                                  tGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACAAAGccc  >  1:23238/1‑150 (MQ=255)
                                                                                                                                                  tGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACAAAGccc  <  2:259630/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCGCTCGTGGTATTGGCGCTGTTCTTCGCTGTTCTGCATGGCGTGACACTGGCGGAAGCCATCGTCCCAGCCGGTTTCGTACTGACGCTCGTGCAGGTAGCGCGGTACGTTCTTGCGGAAGTTGCCCACCATCAGCCCGGCCGCCTGGCGGCCACTGCTGCAGCCGTCCTGGAAGCCATCCGCATAGGCAGGAGGGTAACCGTGGTTGACCATCTGCTCGTGGGTGGTTTCGCAACCCGACAGCAACAACGCAATGCACACGCCGAACCAGTACCGCGACATAGCCACAAAGCCC  >  NC_002947/790814‑791109

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