Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,042,365 C→A 8.4% L137I (CTC→ATC) ‡ PP_3563 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,042,3650CA 8.4% 284.9 / 9.6 92L137I (CTC→ATC) ‡PP_3563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (41/43);  new base A (4/4);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGG  >  NC_002947/4042219‑4042508
                                                                                                                                                  |                                                                                                                                               
tCAGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGTTGGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCa                                                                                                                                              >  1:552088/1‑150 (MQ=255)
tCAGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGCCACGGCCGCTGCTCa                                                                                                                                              >  2:912928/1‑150 (MQ=255)
  aGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGc                                                                                                                                            <  1:566756/150‑1 (MQ=255)
  aGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGc                                                                                                                                            <  2:1487039/150‑1 (MQ=255)
   gCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCg                                                                                                                                           >  1:839143/1‑150 (MQ=255)
   gCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCg                                                                                                                                           <  2:1018740/150‑1 (MQ=255)
   gCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCg                                                                                                                                           >  2:1322136/1‑150 (MQ=255)
   gCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCg                                                                                                                                           >  1:1533220/1‑150 (MQ=255)
      tGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgcc                                                                                                                                        >  1:1274031/1‑150 (MQ=255)
            ggCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTgg                                                                                                                                  <  2:1234903/150‑1 (MQ=255)
            ggCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTgg                                                                                                                                  <  1:155633/150‑1 (MQ=255)
            ggCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTgg                                                                                                                                  >  2:719357/1‑150 (MQ=255)
             gCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTggg                                                                                                                                 <  1:2228169/150‑1 (MQ=255)
              ccACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGc                                                                                                                                <  2:630489/150‑1 (MQ=255)
                 cgcgACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGAc                                                                                                                             >  1:597096/1‑150 (MQ=255)
                 cgcgACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATGCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGAc                                                                                                                             >  2:16824/1‑150 (MQ=255)
                   cgACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTg                                                                                                                           <  1:854613/150‑1 (MQ=255)
                      cgcacgcGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCa                                                                                                                        <  2:839143/150‑1 (MQ=255)
                      cgcacgcGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCa                                                                                                                        <  2:1533220/150‑1 (MQ=255)
                         tcgcgCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCa                                                                                                                        >  2:584073/2‑147 (MQ=255)
                         tcgcgCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCa                                                                                                                        <  1:584073/146‑1 (MQ=255)
                               cgccATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAggg                                                                                                               >  2:2358127/1‑150 (MQ=255)
                                gccATGCATGTCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGa                                                                                                              <  1:189719/150‑1 (MQ=255)
                                  catgcatgCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACa                                                                                                            <  1:1285887/150‑1 (MQ=255)
                                   atgcatgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAg                                                                                                           >  2:2182507/1‑150 (MQ=255)
                                    tgcatgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGa                                                                                                          >  2:1750308/1‑150 (MQ=255)
                                     gcatgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCTGCCTGGGCGACTGGCAGGCCGAGGGACAGaa                                                                                                         >  2:2369151/1‑150 (MQ=255)
                                     gcatgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGaa                                                                                                         <  1:856832/150‑1 (MQ=255)
                                         gcaagcGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGtg                                                                                                     <  2:552088/145‑1 (MQ=255)
                                          ccATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCTGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGtgc                                                                                                    <  2:597096/150‑1 (MQ=255)
                                           cATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGtgct                                                                                                   >  1:199398/1‑150 (MQ=255)
                                            aTCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGtgctg                                                                                                  >  2:2293219/1‑150 (MQ=255)
                                               gAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCag                                                                                               <  2:1802310/150‑1 (MQ=255)
                                               gAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCag                                                                                               >  2:1456737/1‑150 (MQ=255)
                                                 ggAAGGATTCTGTCAGGTGGTGCTCGACCCGCAATTTGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCGCATCAGCAGCGGCGGCCTGGGCGACTGGCAGGCCTAGGGACAGAAGGTGCTGCagga                                                                                             <  2:338623/150‑1 (MQ=255)
                                                 ggAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGc                                                                                                                                >  2:1427773/1‑115 (MQ=255)
                                                 ggAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGc                                                                                                                                <  1:1427773/115‑1 (MQ=255)
                                                 ggAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCagga                                                                                             <  1:376306/150‑1 (MQ=255)
                                                  gAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggag                                                                                            <  1:2162854/150‑1 (MQ=255)
                                                  gAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggag                                                                                            <  1:943397/150‑1 (MQ=255)
                                                    aaGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGt                                                                                          <  2:46670/150‑1 (MQ=255)
                                                             cgaGGTGGTGCTCGACCCGCCATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGGCTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAAcgccgc                                                                                 <  2:654056/148‑1 (MQ=255)
                                                             ccAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCgctgct                                                                                                                                                <  2:75513/87‑1 (MQ=255)
                                                             ccAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCgctgct                                                                                                                                                >  1:75513/1‑87 (MQ=255)
                                                               aGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCcgcgc                                                                               <  1:2293219/150‑1 (MQ=255)
                                                                  tggtgCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCa                                                                            >  1:1928344/1‑150 (MQ=255)
                                                                     tgCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCcgcgc                                                                               <  2:944862/144‑1 (MQ=255)
                                                                     tgCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCcgcgc                                                                               >  1:944862/1‑144 (MQ=255)
                                                                       cTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccgt                                                                                                                                      >  1:867917/1‑86 (MQ=255)
                                                                       cTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccgt                                                                                                                                      <  2:867917/87‑2 (MQ=255)
                                                                          gACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCa                                                                    <  1:78984/150‑1 (MQ=255)
                                                                             ccGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACgg                                                                 >  1:1617612/1‑150 (MQ=255)
                                                                             ccGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACgg                                                                 >  1:2073086/1‑150 (MQ=255)
                                                                              cGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccgtgc                                                                                                                                    >  1:1235963/1‑79 (MQ=38)
                                                                              cGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccgtgc                                                                                                                                    <  2:1235963/82‑4 (MQ=39)
                                                                                    tcgtcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCaa                                                              >  2:2294895/1‑146 (MQ=255)
                                                                                    tcgtcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCaa                                                              <  1:2294895/146‑1 (MQ=255)
                                                                                      gtcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATgcggc                                                        >  1:1751057/1‑150 (MQ=255)
                                                                                       tcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGAc                                                                                                                             <  1:757879/80‑1 (MQ=255)
                                                                                       tcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGAc                                                                                                                             >  2:757879/1‑80 (MQ=255)
                                                                                         gATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCggcggc                                                     <  2:199398/150‑1 (MQ=255)
                                                                                         gATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCggcggc                                                     <  1:2279390/150‑1 (MQ=255)
                                                                                             tCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgcc                                                                                                                                        <  1:1119585/63‑1 (MQ=255)
                                                                                             tCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgcc                                                                                                                                        >  2:1119585/1‑63 (MQ=255)
                                                                                             tCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGt                                                                        <  1:641691/127‑1 (MQ=255)
                                                                                             tCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGt                                                                        >  2:641691/1‑127 (MQ=255)
                                                                                                gCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGt                                              >  1:726708/1‑150 (MQ=255)
                                                                                                 cGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGtt                                             <  1:1456737/150‑1 (MQ=255)
                                                                                                  gTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCagga                                                                                             <  2:386991/101‑1 (MQ=255)
                                                                                                  gTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCagga                                                                                             >  1:386991/1‑101 (MQ=255)
                                                                                                     aGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgcc                                         <  2:1751057/150‑1 (MQ=255)
                                                                                                      gCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgccg                                        >  1:2275701/1‑150 (MQ=255)
                                                                                                      gCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccg                                                                                                                                       <  2:503482/55‑1 (MQ=38)
                                                                                                      gCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccg                                                                                                                                       >  1:503482/1‑55 (MQ=38)
                                                                                                       ccAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGttt                                            <  1:72493/145‑1 (MQ=255)
                                                                                                       ccAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGttt                                            >  2:72493/1‑145 (MQ=255)
                                                                                                        cAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgccgcc                                      >  1:1131797/1‑150 (MQ=255)
                                                                                                        cAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgccgcc                                      >  1:1701484/1‑150 (MQ=255)
                                                                                                           gaggGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgcc                                                                                                                                        >  2:1654713/1‑49 (MQ=255)
                                                                                                           gaggGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgcc                                                                                                                                        <  1:1654713/49‑1 (MQ=255)
                                                                                                                         gcgcATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTg                     >  1:2036783/1‑150 (MQ=255)
                                                                                                                         gcgcATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGAAGGCCGAGGGAAAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTg                     <  1:16824/150‑1 (MQ=255)
                                                                                                                                 tGCTGGCACGGCCGCTGCTCAGCGGCTGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGAAAGGTGCCACGGCAATGCGGGGGCCAGCGTTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAg             <  1:2369151/150‑1 (MQ=255)
                                                                                                                                   cccaggcGGCCGCTGAGCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACATAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTg           >  1:552564/7‑150 (MQ=255)
                                                                                                                                   cTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCGGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTg           <  2:726708/150‑1 (MQ=255)
                                                                                                                                      gCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTg        >  2:922454/1‑150 (MQ=255)
                                                                                                                                        aCGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCACTGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGc      >  1:183431/1‑150 (MQ=255)
                                                                                                                                        aCGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCa                             <  2:2228712/127‑1 (MQ=255)
                                                                                                                                        aCGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCa                             >  1:2228712/1‑127 (MQ=255)
                                                                                                                                         cGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGTACAGATGGTGCTGCAGAGGGTCATCGCAGCGCGCAAGGTGCCACAGCAATGCCGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTAGAGCACAGTGCCGGcc     >  2:84541/1‑150 (MQ=255)
                                                                                                                                           gCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGt                                                                        >  2:211339/1‑81 (MQ=255)
                                                                                                                                           gCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGt                                                                        <  1:211339/81‑1 (MQ=255)
                                                                                                                                            ccGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCggg  >  1:222727/1‑150 (MQ=255)
                                                                                                                                            ccGCGGCTCAGCGGCCCCCTGGGCGCCTGCCCGGCCGAGGAATAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCggg  <  1:1812901/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TCAGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGG  >  NC_002947/4042219‑4042508

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