Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,703,316:1 +C 5.6% coding (97/270 nt) PP_5596 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,703,3161.C 5.6% 294.8 / 10.2 71coding (97/270 nt)PP_5596hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (37/30);  new base C (2/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CCACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGCCCAACCAGTCGGTG  >  NC_002947/3703167‑3703465
                                                                                                                                                      |                                                                                                                                                     
ccACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCt                                                                                                                                                        <  1:879391/150‑1 (MQ=255)
 cACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTg                                                                                                                                                      <  1:521113/150‑1 (MQ=255)
 cACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTg                                                                                                                                                      <  1:1040966/150‑1 (MQ=255)
  aacGCCCGATCAATCAAAATGTCATTCATGGAGAAGAGGAGGCCGCCTTCGATGCCTAGCCAACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTgg                                                                                                                                                     >  2:342707/3‑150 (MQ=255)
        cGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTgg                                                                                                                                               >  2:2338980/1‑150 (MQ=255)
        cGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTgg                                                                                                                                               <  2:1096612/150‑1 (MQ=255)
               tcaaAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGCCATTGGGCAAGtt                                                                                                                                        >  2:736936/1‑150 (MQ=255)
                caaAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTtg                                                                                                                                       <  1:1705762/150‑1 (MQ=255)
                 aaaaTGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTtgt                                                                                                                                      >  2:2377520/1‑150 (MQ=255)
                 aaaaTGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTtgt                                                                                                                                      >  2:1893394/1‑150 (MQ=255)
                        cATCCATGGAGAAGAGGAGGCCCCCTGCAATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTTGACAATACCGAGGCCTGGCATTGGGCAAGTTGTGATGCgg                                                                                                                               <  1:2337747/150‑1 (MQ=255)
                         aTCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCggg                                                                                                                              >  1:2074843/1‑150 (MQ=255)
                                gagaAGAGGAGGCCCCCTGCGTTGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCCAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGc                                                                                                                       <  2:1397028/150‑1 (MQ=255)
                                     gaggagGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTgt                                                                                                                                               >  2:35352/1‑120 (MQ=255)
                                     gaggagGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTgt                                                                                                                                               <  1:35352/121‑2 (MQ=255)
                                     gaggagGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGcatcat                                                                                                                  >  1:80664/1‑150 (MQ=255)
                                           gCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCaa                                                                                                            >  1:1550297/1‑150 (MQ=255)
                                              cccTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGAAGTTTCATCATCCGCCTGGAACGCCTGGAGAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGtt                                                                                                         <  1:1105139/150‑1 (MQ=255)
                                                         cTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCgaga                                                                                              >  2:1912925/1‑150 (MQ=255)
                                                         cTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCAGGCCTCGGCATTgtcc                                                                                                                                             >  1:988415/1‑100 (MQ=255)
                                                         cTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCAGGCCTCGGCATTgtcc                                                                                                                                             <  2:988415/103‑4 (MQ=255)
                                                                         cATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAg                                                                              >  1:1548780/1‑150 (MQ=255)
                                                                         cATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGACCTGAGATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCGAGTTCGGTCGCGAGAAAATCAGGAAGGTCAg                                                                              >  2:693807/1‑150 (MQ=255)
                                                                            aaTCGATTACAACCTAAATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCg                                                                           >  2:2394215/1‑150 (MQ=255)
                                                                            aaTCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCg                                                                           >  1:1451245/1‑150 (MQ=255)
                                                                             aTCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCgg                                                                          <  1:1912925/150‑1 (MQ=255)
                                                                              tCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGc                                                                         >  1:1283837/1‑150 (MQ=255)
                                                                              tCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGc                                                                         <  1:1514574/150‑1 (MQ=255)
                                                                                 aTTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGaa                                                                                    >  2:1939665/1‑136 (MQ=255)
                                                                                 aTTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGaa                                                                                    <  1:1939665/136‑1 (MQ=255)
                                                                                      aaaCTGCATGACCAAGAGCGCAGTTTGATCATCCGCCTGGAACGCCTGGACAATGGCGGGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGTTCGCGAGACAATCAGGAAGGTCAGTCGGCAGATCGcc                                                                 >  2:826389/1‑150 (MQ=255)
                                                                                       aaCTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGcatc                                                                                                                    >  1:1801883/1‑98 (MQ=255)
                                                                                       aaCTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGcatc                                                                                                                    <  2:1801883/98‑1 (MQ=255)
                                                                                       aaCTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCg                                                                >  2:2284313/1‑150 (MQ=255)
                                                                                        aCTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTc                                                                                                        <  1:1192280/109‑1 (MQ=255)
                                                                                        aCTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTc                                                                                                        >  2:1192280/1‑109 (MQ=255)
                                                                                         cTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGaa                                                                                    >  1:32401/1‑128 (MQ=255)
                                                                                         cTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGaa                                                                                    <  2:32401/128‑1 (MQ=255)
                                                                                         cTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAg                                                              >  1:399901/1‑150 (MQ=255)
                                                                                          tGCATGACCAACCGCGCAGTTTCCTCATCCGCCTGGGTCGCCTTGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTTATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGc                                                             <  1:2338980/150‑1 (MQ=255)
                                                                                          tGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGc                                                             >  2:1628840/1‑150 (MQ=255)
                                                                                             aTGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCaa                                                                                                                                           >  1:1898534/1‑69 (MQ=255)
                                                                                             aTGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCaa                                                                                                                                           <  2:1898534/69‑1 (MQ=255)
                                                                                              tGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTa                                                         >  2:1454455/1‑150 (MQ=255)
                                                                                                  cAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGc                                                     >  2:2223623/1‑150 (MQ=255)
                                                                                                   aaCCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCa                                                    >  2:802382/1‑150 (MQ=255)
                                                                                                       gcgcAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTAc                                                >  2:1979467/1‑150 (MQ=255)
                                                                                                          cAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTAccgc                                             <  1:1454455/150‑1 (MQ=255)
                                                                                                             tttCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGt                                          <  1:802382/150‑1 (MQ=255)
                                                                                                               tcatcaTCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGt                                          <  2:1187423/148‑1 (MQ=255)
                                                                                                               tcatcaTCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGt                                          >  1:1187423/1‑148 (MQ=255)
                                                                                                                  tcatcCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGt                                                                             >  2:500369/1‑110 (MQ=255)
                                                                                                                  tcatcCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGt                                                                             <  1:500369/110‑1 (MQ=255)
                                                                                                                     tcCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGc                                  >  1:917056/1‑150 (MQ=255)
                                                                                                                     tcCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGc                                  <  1:1979467/150‑1 (MQ=255)
                                                                                                                     tcCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGc                                  >  2:1579022/1‑150 (MQ=255)
                                                                                                                             gAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCg                          <  2:1334759/150‑1 (MQ=255)
                                                                                                                              aaCGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCgg                         <  2:971114/150‑1 (MQ=255)
                                                                                                                              aaCGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCgg                         <  2:1548780/150‑1 (MQ=255)
                                                                                                                              aaCGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCgg                         <  1:1579022/150‑1 (MQ=255)
                                                                                                                               aCGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCggg                        <  1:1628840/150‑1 (MQ=255)
                                                                                                                                cGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCgggg                       <  1:1156646/150‑1 (MQ=255)
                                                                                                                                   cTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTAccgcc                                            >  2:912823/1‑126 (MQ=255)
                                                                                                                                   cTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTAccgcc                                            <  1:912823/126‑1 (MQ=255)
                                                                                                                                   cTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGc                    >  1:591118/1‑150 (MQ=255)
                                                                                                                                      gACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAg                 <  2:917056/150‑1 (MQ=255)
                                                                                                                                       ccAATGCCAGGCCTCGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCg                          <  2:1802909/139‑1 (MQ=255)
                                                                                                                                       ccAATGCCAGGCCTCGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATGTCg                          >  1:1802909/2‑140 (MQ=255)
                                                                                                                                         aaTGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGccc              >  2:1840168/1‑150 (MQ=255)
                                                                                                                                         aaTGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGccc              >  2:1246456/1‑150 (MQ=255)
                                                                                                                                            gCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGCCCAAc           <  1:106232/150‑1 (MQ=255)
                                                                                                                                              cGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGt                                          <  1:1916801/117‑1 (MQ=255)
                                                                                                                                              cGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGt                                          >  2:1916801/1‑117 (MQ=255)
                                                                                                                                              cGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGa                      >  1:1179363/1‑137 (MQ=255)
                                                                                                                                              cGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGa                      <  2:1179363/137‑1 (MQ=255)
                                                                                                                                                     tGGCATTGGGCAAGTTTTGATGCGGGCATTGGCATCATTCCCAATTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGCCCAACCAGTCGGTg  <  1:342707/150‑1 (MQ=255)
                                                                                                                                                     tGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGa                      <  2:147514/130‑1 (MQ=255)
                                                                                                                                                     tGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCACTTCCGCCGTGAGCTGGCGGATTTCGGGGa                      >  1:147514/1‑130 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCACCGCCCGATCAATCAAAATGTCATCCATGGAGAAGAGGAGGCCCCCTGCGATGCCTAGCCCACAGCTGTACATAATCGATTACAAACTGCATGACCAACCGCGCAGTTTCATCATCCGCCTGGAACGCCTGGACAATGCCGAGGCCTGGCATTGGGCAAGTTGTGATGCGGGCATTGGCATCATTCCCAAGTTCGGTCGCGAGAAAATCAGGAAGGTCAGTCGGCCGATGGCCGAGCGTTACGGCATTACCGCCGTGAGCTGGCGGATTTCGGGGAGCAAGCCCAACCAGTCGGTG  >  NC_002947/3703167‑3703465

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: