Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,263,103 T→C 6.2% H437R (CAC→CGC) ‡ PP_1105 ← ATP‑dependent DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,263,1030TC 6.2% 241.1 / 4.2 81H437R (CAC→CGC) ‡PP_1105ATP‑dependent DNA ligase
Reads supporting (aligned to +/- strand):  ref base T (40/36);  new base C (2/3);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

TGCGGATGATGGCGTCGACCTTGCGCATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGT  >  NC_002947/1262954‑1263249
                                                                                                                                                     |                                                                                                                                                  
tGCGGATGATGGCGTCGACCTTGCGCATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCCGCCGGTGCATCCCATCCGGCGAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGt                                                                                                                                                    <  2:1642978/150‑1 (MQ=255)
        atgGCGTCGACCTTGCGCATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCACGGCcgcag                                                                                                                                            <  2:1639557/150‑3 (MQ=255)
         tgGCGTCGACCTTGCGCATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGtt                                                                                                                                           <  2:1136145/150‑1 (MQ=255)
           gCGTCGACCTTGCGCATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTgg                                                                                                                                         <  2:1083233/150‑1 (MQ=255)
                     tGCGCATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCa                                                                                                                               <  2:1903137/150‑1 (MQ=255)
                      gcgcATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAg                                                                                                                              <  2:1690020/150‑1 (MQ=255)
                        gcATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCa                                                                                                                            >  1:1448855/1‑150 (MQ=255)
                         cATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCAc                                                                                                                           >  2:1566738/1‑150 (MQ=255)
                         cATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCAc                                                                                                                           >  1:1457481/1‑150 (MQ=255)
                          aTTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCAcc                                                                                                                          >  1:1776455/1‑150 (MQ=255)
                             tcttcGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCa                                                                                                                       <  2:838327/150‑1 (MQ=255)
                               ttcGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATc                                                                                                                     >  1:309423/1‑150 (MQ=255)
                               ttcGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCACGGCcgcagg                                                                                                                                           >  2:153354/1‑125 (MQ=255)
                               ttcGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGACCTGCGGCCACGGCcgcagg                                                                                                                                           <  1:153354/128‑4 (MQ=255)
                                   tcgCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACg                                                                                                                 >  1:1580397/1‑150 (MQ=255)
                                         aGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGaaa                                                                                                           <  2:1580397/150‑1 (MQ=255)
                                           cccTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCCGCACCGCATCGACGCTGAAAgg                                                                                                         >  1:97776/1‑150 (MQ=255)
                                            ccTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAggg                                                                                                        <  1:1335799/150‑1 (MQ=255)
                                            ccTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAggg                                                                                                        >  1:129164/1‑150 (MQ=255)
                                             cTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGt                                                                                                       <  1:1566738/150‑1 (MQ=255)
                                             cTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGt                                                                                                       >  1:1603571/1‑150 (MQ=255)
                                             cTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGt                                                                                                       <  1:791836/150‑1 (MQ=255)
                                              tGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGTTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGAGTGGCCCTGCGGCCGTgg                                                                                                                                                  <  2:1606561/106‑1 (MQ=255)
                                              tGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGAGTGGCCCTGCGGCCGTgg                                                                                                                                                  >  1:1606561/1‑106 (MQ=255)
                                              tGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAGAGGGTc                                                                                                      >  1:169324/1‑150 (MQ=255)
                                                   aGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATct                                                                                                 <  2:1459844/150‑1 (MQ=255)
                                                     gCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGGATCGACGCTGAAAGGGTCGATcttc                                                                                               <  1:1485794/150‑1 (MQ=255)
                                                     gCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATcttc                                                                                               <  2:129164/150‑1 (MQ=255)
                                                     gCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATcttc                                                                                               <  2:97776/150‑1 (MQ=255)
                                                     gCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATcttc                                                                                               <  1:1485793/150‑1 (MQ=255)
                                                       cTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGAtgcgg                                                                                               >  2:774689/1‑144 (MQ=255)
                                                       cTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGAtgcgg                                                                                               <  1:774689/148‑5 (MQ=255)
                                                             aaaGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATcttccacttcca                                                                                       >  2:1727772/1‑150 (MQ=255)
                                                                    aCCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATa                                                                                                                                    <  1:663568/98‑1 (MQ=255)
                                                                    aCCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATa                                                                                                                                    >  2:663568/1‑98 (MQ=255)
                                                                      cAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGcc                                                                              >  1:661938/1‑150 (MQ=255)
                                                                      cAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGcc                                                                              >  1:944627/1‑150 (MQ=255)
                                                                       aaCGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATGTAGCAGCTCCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGccc                                                                             <  2:927803/150‑1 (MQ=255)
                                                                        aCGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTg                                                                                                                                                   <  2:1353517/79‑1 (MQ=255)
                                                                        aCGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTg                                                                                                                                                   >  1:1353517/1‑79 (MQ=255)
                                                                          gCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATg                                                                          >  1:31981/1‑150 (MQ=255)
                                                                           cccTACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGt                                                                         >  1:626120/1‑150 (MQ=255)
                                                                           cccGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGt                                                                         >  2:1587130/1‑150 (MQ=255)
                                                                           cccGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGt                                                                         <  1:1071521/150‑1 (MQ=255)
                                                                            ccGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATcttc                                                                                               <  2:622298/127‑1 (MQ=255)
                                                                            ccGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATcttc                                                                                               >  1:622298/1‑127 (MQ=255)
                                                                            ccGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTc                                                                        >  2:1714837/1‑150 (MQ=255)
                                                                            ccGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTc                                                                        >  2:469009/1‑150 (MQ=255)
                                                                               aCCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCtt                                                                     <  2:1280077/150‑1 (MQ=255)
                                                                               aCCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCACGGCcgcag                                                                                                                                            <  2:99227/79‑3 (MQ=255)
                                                                               aCCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCACGGCcgcag                                                                                                                                            >  1:99227/1‑77 (MQ=255)
                                                                                   ggggTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTg                                                                 >  2:579746/1‑150 (MQ=255)
                                                                                    gggTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTgc                                                                <  2:1457481/150‑1 (MQ=255)
                                                                                      gTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGAt                                                                                                   >  1:24783/1‑113 (MQ=255)
                                                                                      gTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGAt                                                                                                   <  2:24783/113‑1 (MQ=255)
                                                                                         ccTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGcccc                                                           <  2:919740/150‑1 (MQ=255)
                                                                                          cTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCca                                                          <  1:638709/150‑1 (MQ=255)
                                                                                                     tCCCATCCGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGCAAGGGTCGCTCTTCCCCTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAgg                                               >  2:747054/1‑150 (MQ=255)
                                                                                                     tCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCGCACCGTACAgg                                               <  1:1272052/150‑1 (MQ=255)
                                                                                                        cATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTgg                                                                                                                                         <  2:1233232/57‑1 (MQ=255)
                                                                                                        cATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTgg                                                                                                                                         >  1:1233232/1‑57 (MQ=255)
                                                                                                                 aaGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGc                                   <  2:309423/150‑1 (MQ=255)
                                                                                                                  aGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCTGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCt                                  >  2:931127/1‑150 (MQ=255)
                                                                                                                    gtgtAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATcttc                                                                                               <  1:1042497/87‑1 (MQ=255)
                                                                                                                    gtgtAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATcttc                                                                                               >  2:1042497/1‑87 (MQ=255)
                                                                                                                    gtgtAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCACTCCCGCATCGACGCTGAAAGGGGCGAGCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTc                                >  1:1673789/1‑150 (MQ=255)
                                                                                                                    gtgtAGCAGCTGCACAAACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATAAGAACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTc                                >  2:274717/1‑150 (MQ=255)
                                                                                                                     tgtAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCa                               <  1:579746/150‑1 (MQ=255)
                                                                                                                       tAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGc                             >  2:1661253/1‑150 (MQ=255)
                                                                                                                           cGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCa                         <  1:1727772/150‑1 (MQ=255)
                                                                                                                            gCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAg                        >  2:918699/1‑150 (MQ=255)
                                                                                                                              tGTACAGACTGGCCCTGCGGCCGTGGCCGTGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATACCCATGTCCTTGGTTCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGcc                      >  2:524129/1‑150 (MQ=255)
                                                                                                                                  cAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTc                  <  1:437178/150‑1 (MQ=255)
                                                                                                                                     aCTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCAc               <  1:918699/150‑1 (MQ=255)
                                                                                                                                       tGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGc             >  1:60150/1‑150 (MQ=255)
                                                                                                                                       tGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGc             <  2:1603571/150‑1 (MQ=255)
                                                                                                                                       tGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGAACATATCAATGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGc             >  2:56754/1‑150 (MQ=255)
                                                                                                                                         gCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCTACGccc           >  2:1079258/1‑150 (MQ=255)
                                                                                                                                         gCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGccc           >  1:1105735/1‑150 (MQ=255)
                                                                                                                                           ccTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAg         >  1:1223046/1‑150 (MQ=255)
                                                                                                                                           ccTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAg         >  1:404338/1‑150 (MQ=255)
                                                                                                                                                  ccGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGt  <  2:1105735/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGCGGATGATGGCGTCGACCTTGCGCATTTCTTCGTCGCTGAGCCCTGAATAGGCCTTGGCAAAGGGCACCAACGCCCGACCAGGGGTACCTGCCGGTGCATCCCATACGGCAAAGGTGTAGTCGCTGTACAGACTGGCCCTGCGGCCGTGGCCGCGTTGGGCATAGATCAGCACCGCATCGACGCTGAAAGGGTCGATCTTCCACTTCCACCAGATGCCCATGTCCTTGGTGCGCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGT  >  NC_002947/1262954‑1263249

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