Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 954,436 T→C 5.6% L92P (CTG→CCG)  cyoE → protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947954,4360TC 5.6% 2396.0 / 77.2 744L92P (CTG→CCG) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base T (399/301);  new base C (17/25);  total (418/326)
Fisher's exact test for biased strand distribution p-value = 3.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 747 total aligned reads displayed.

GGTGATCGGTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGCCTGT  >  NC_002947/954287‑954585
                                                                                                                                                     |                                                                                                                                                     
    aTCGGTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTc                                                                                                                                                   >  2:2818154/1‑150 (MQ=255)
     tCGGTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCg                                                                                                                                                  <  1:1200510/150‑1 (MQ=255)
     tCGGTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCg                                                                                                                                                  >  2:5281003/1‑150 (MQ=255)
      cGGTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCgg                                                                                                                                                 <  1:6598532/150‑1 (MQ=255)
        gTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCgggg                                                                                                                                               >  1:7991187/1‑150 (MQ=255)
        gTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCgggg                                                                                                                                               >  2:7980502/1‑150 (MQ=255)
        gTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCggg                                                                                                                                                <  2:4530820/149‑1 (MQ=255)
           cTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCCGAGCAGGGTgg                                                                                                                                            <  1:2772292/150‑1 (MQ=255)
           cTTCGCGGGTGGGTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGAGGGAGCGCCCCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGCTCGACGCCACCCTGCTCGGGGTgg                                                                                                                                            <  2:4081404/150‑1 (MQ=255)
            ttCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGc                                                                                                                                           >  2:7519148/1‑150 (MQ=255)
            ttCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGc                                                                                                                                           >  1:3730906/1‑150 (MQ=255)
              cGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCgg                                                                                                                                         <  2:3574685/150‑1 (MQ=255)
               gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCggg                                                                                                                                        <  1:7004110/150‑1 (MQ=255)
               gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCggg                                                                                                                                        <  1:8575506/150‑1 (MQ=255)
               gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCggg                                                                                                                                        <  2:2429963/150‑1 (MQ=255)
               gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCCGAGCAGGGTGGCgta                                                                                                                                        <  1:7601837/150‑3 (MQ=255)
                  ggtggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGttt                                                                                                                                     <  1:1532409/150‑1 (MQ=255)
                  ggtggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGttt                                                                                                                                     <  2:1342652/150‑1 (MQ=255)
                  ggtggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCCGAGCAGGGTGGCgtagat                                                                                                                                     <  1:3430151/150‑6 (MQ=255)
                   gtggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTc                                                                                                                                    <  1:6465351/150‑1 (MQ=255)
                    tggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCa                                                                                                                                   <  1:6246935/150‑1 (MQ=255)
                    tgggTGCGTCCGGATGCGTGGTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGACCCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCa                                                                                                                                   <  2:1215935/150‑1 (MQ=255)
                       ttGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGcc                                                                                                                                >  1:1856827/1‑150 (MQ=255)
                            tCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCctgctg                                                                                                                           <  2:7183445/150‑1 (MQ=255)
                            tCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCctgctg                                                                                                                           >  1:4788480/1‑150 (MQ=255)
                             ccGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGt                                                                                                                          >  1:2564220/1‑150 (MQ=255)
                               ggATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTAt                                                                                                                        >  1:5987827/1‑150 (MQ=255)
                               ggATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTAt                                                                                                                        >  2:6799777/1‑150 (MQ=255)
                                     gtgtTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAg                                                                                                                  <  2:339488/150‑1 (MQ=255)
                                     gtgtTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAg                                                                                                                  >  2:5647708/1‑150 (MQ=255)
                                       gtTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGc                                                                                                                <  1:7842499/150‑1 (MQ=255)
                                        ttCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGcc                                                                                                               >  2:7869729/1‑150 (MQ=255)
                                           aacaacTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAAc                                                                                                            >  2:3898108/1‑150 (MQ=255)
                                           aacaacTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAAc                                                                                                            >  2:8672602/1‑150 (MQ=255)
                                                cTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCt                                                                                                       >  1:4033790/1‑150 (MQ=255)
                                                 tGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTg                                                                                                      >  2:1735823/1‑150 (MQ=255)
                                                  gCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGt                                                                                                     <  1:3274616/150‑1 (MQ=255)
                                                   cATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTc                                                                                                    >  2:2935420/1‑150 (MQ=255)
                                                   cATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTc                                                                                                    >  1:2883416/1‑150 (MQ=255)
                                                     tCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTg                                                                                                  >  1:8486225/1‑150 (MQ=255)
                                                     tCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTg                                                                                                  >  1:527355/1‑150 (MQ=255)
                                                      cGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGc                                                                                                 <  1:5710266/150‑1 (MQ=255)
                                                       gACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCg                                                                                                <  1:5565983/150‑1 (MQ=255)
                                                         ccGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGtt                                                                                              >  2:4639704/1‑150 (MQ=255)
                                                             gACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTgc                                                                                          >  1:3973662/1‑150 (MQ=255)
                                                             gACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTgc                                                                                          <  2:8289819/150‑1 (MQ=255)
                                                               cATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTgcgc                                                                                        <  1:1104760/150‑1 (MQ=255)
                                                                aTCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCa                                                                                       >  2:4825248/1‑150 (MQ=255)
                                                                aTCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCa                                                                                       <  1:898269/150‑1 (MQ=255)
                                                                   gACCACAAGATGGAGCGCACCAAGAACCGTGTCGTGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTg                                                                                    >  1:24202/1‑150 (MQ=255)
                                                                     ccACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGAt                                                                                  <  1:7338479/150‑1 (MQ=255)
                                                                      cacaAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGAtt                                                                                 <  1:5429516/150‑1 (MQ=255)
                                                                       acaAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTg                                                                                >  2:4194352/1‑150 (MQ=255)
                                                                         aaGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGc                                                                              <  2:629800/150‑1 (MQ=255)
                                                                               gAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTgtc                                                                        <  2:5994648/150‑1 (MQ=255)
                                                                                  cgcACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTgtcgtc                                                                     <  1:4475314/150‑1 (MQ=255)
                                                                                   gcACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTGTGTCGTCt                                                                    >  2:5889760/1‑150 (MQ=255)
                                                                                    cACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTa                                                                   <  1:189540/150‑1 (MQ=255)
                                                                                    cACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTa                                                                   <  2:3058493/150‑1 (MQ=255)
                                                                                     aCCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTAc                                                                  >  2:5095079/1‑150 (MQ=255)
                                                                                       cAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGt                                                                >  1:2972740/1‑150 (MQ=255)
                                                                                        aaGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTc                                                               <  1:211759/150‑1 (MQ=255)
                                                                                        aaGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGAGTGGCTTTGTCGTCTACGTc                                                               >  2:5169303/1‑150 (MQ=255)
                                                                                         tgAACCATGACACGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCg                                                              >  2:6506495/2‑150 (MQ=255)
                                                                                             ccGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGttt                                                          >  2:1954452/1‑150 (MQ=255)
                                                                                             ccGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGttt                                                          >  1:8299991/1‑150 (MQ=255)
                                                                                                  tCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACa                                                     <  1:1706760/150‑1 (MQ=255)
                                                                                                     tGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGcc                                                  <  2:1004011/150‑1 (MQ=255)
                                                                                                     tGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGcc                                                  <  1:7763372/150‑1 (MQ=255)
                                                                                                       gTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCtg                                                >  1:3157278/1‑150 (MQ=255)
                                                                                                         tCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCtgtg                                              >  1:4296890/1‑150 (MQ=255)
                                                                                                            gggCGGCCTGTCGCAGCCCCTCGCGCAGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCt                                           <  1:4976742/150‑1 (MQ=255)
                                                                                                             ggcggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTg                                          >  1:4564025/1‑150 (MQ=255)
                                                                                                               cggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGaa                                        <  2:3268614/150‑1 (MQ=255)
                                                                                                               cggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGaa                                        <  2:5494035/150‑1 (MQ=255)
                                                                                                                ggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAg                                       <  2:736288/150‑1 (MQ=255)
                                                                                                                ggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAg                                       <  2:8135573/150‑1 (MQ=255)
                                                                                                                ggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAg                                       <  1:7138782/150‑1 (MQ=255)
                                                                                                                ggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAg                                       <  2:3704744/150‑1 (MQ=255)
                                                                                                                 gcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGc                                      <  1:3571438/150‑1 (MQ=255)
                                                                                                                  cATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCg                                     <  1:2262756/150‑1 (MQ=255)
                                                                                                                       cGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAAt                                <  1:7486198/150‑1 (MQ=255)
                                                                                                                         cTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCg                              <  2:7058423/150‑1 (MQ=255)
                                                                                                                            ccccTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTg                           <  2:157135/150‑1 (MQ=255)
                                                                                                                              ccTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCa                         >  1:1021183/1‑150 (MQ=255)
                                                                                                                               cTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCAc                        <  2:8395945/150‑1 (MQ=255)
                                                                                                                                tCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACg                       >  2:6943439/1‑150 (MQ=255)
                                                                                                                                  gcgcTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGc                     >  1:8174726/1‑150 (MQ=255)
                                                                                                                                     cTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCAcc                  <  1:6110841/150‑1 (MQ=255)
                                                                                                                                      tGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCt                 >  2:2240062/1‑150 (MQ=255)
                                                                                                                                        aTCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTg               <  2:7259359/150‑1 (MQ=255)
                                                                                                                                          acccGCCACCCCGAGCAGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGt             >  1:417004/4‑150 (MQ=255)
                                                                                                                                             cGCCACCCCGAGCAGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCgg          >  1:5860168/1‑150 (MQ=255)
                                                                                                                                              gCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGc         >  1:5147649/1‑150 (MQ=255)
                                                                                                                                                cACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAg       >  1:6861349/1‑150 (MQ=255)
                                                                                                                                                cACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAg       >  2:7592392/1‑150 (MQ=255)
                                                                                                                                                cACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAg       >  2:44314/1‑150 (MQ=255)
                                                                                                                                                 aCCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGc      >  2:6445981/1‑150 (MQ=255)
                                                                                                                                                   ccTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGCCt    >  1:3858994/1‑150 (MQ=255)
                                                                                                                                                     tGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAcccgac  <  2:2596055/150‑4 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGATCGGTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGCCTGT  >  NC_002947/954287‑954585

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