Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,934,883 Δ1 bp 5.1% coding (1752/2130 nt) flhA ← flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,934,8830C. 5.1% 380.8 / 13.6 117coding (1752/2130 nt)flhAflagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base C (56/55);  new base . (3/3);  total (59/58)
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.94e-01

Only 100 of 121 total aligned reads displayed.

ACCGTCTTCCTGACCCTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAG  >  NC_002947/4934736‑4935029
                                                                                                                                                   |                                                                                                                                                   
 ccGTCTTCCTGACCCTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGccacc                                                                                                                                                  >  2:1521006/1‑150 (MQ=255)
  cGTCTTCCTGACCCTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTACTCACGTGTAGCTCCGACTCAAAGCAAAAAGTGCGGTAGGTGTGTGCTCTACACAATCCGACGCGCACCGCCGccacca                                                                                                                                                 >  1:541295/1‑150 (MQ=255)
       tCCTGACCCTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTgcg                                                                                                                                           >  1:1174819/1‑150 (MQ=255)
         cTGACCCTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTgcggc                                                                                                                                         <  1:2035771/150‑1 (MQ=255)
          gtgccgggtgtcatttccttgtcggggctgctcaaggtgctccaagccctgttgtcggaacaggtgccagtgcgcgatattcgcagtattgccgaggccatcgcgaacaacgcgggcaagagtcaagaTACCGCC‑‑‑ACCAGTGCGGCGGTa                                                                                                                                     >  1:3361‑M1/129‑150 (MQ=255)
           gACCCTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGcggcgg                                                                                                                                       >  1:2330729/1‑150 (MQ=255)
             cccTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCATCGCCGCCACCAGTGCGGCGGTa                                                                                                                                     <  1:2231483/150‑1 (MQ=255)
              ccTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGta                                                                                                                                              >  1:1060820/1‑140 (MQ=255)
              ccTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGta                                                                                                                                              <  2:1060820/141‑2 (MQ=255)
              ccTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTAt                                                                                                                                    >  2:1992442/1‑150 (MQ=255)
                 gCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCtt                                                                                                                                 >  1:1817969/1‑150 (MQ=255)
                   ccGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGa                                                                                                                               >  2:1239653/1‑150 (MQ=255)
                        ccTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCtt                                                                                                                          >  1:1085881/1‑150 (MQ=255)
                         cTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTg                                                                                                                         <  1:1694013/150‑1 (MQ=255)
                                       cAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGGGGTAGGTTGACTCTTGCCCGCGTTGTTcgc                                                                                                           <  2:93202/150‑1 (MQ=255)
                                       cAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTCCCCACcttgtcctt                                                                                                           <  2:1800608/150‑5 (MQ=255)
                                         gCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGa                                                                                                         <  2:1085881/150‑1 (MQ=255)
                                          cAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATg                                                                                                       >  1:1798293/1‑150 (MQ=255)
                                             ggtgctgcaagccctgttgtcggaacaggtgccagtgcgcgatattcgcagtattgccgaggccatcgcgaacaacgcgggcaagagtcaagaTACCGCCGCCACCAGTGCGGCGGTATc                                                                                                                                   >  1:32074‑M1/94‑120 (MQ=255)
                                             ggtgctgcaagccctgttgtcggaacaggtgccagtgcgcgatattcgcagtattgccgaggccatcgcgaacaacgcgggcaagagtcaagaTACCGCCGCCACCAGTGCGGCGGTATc                                                                                                                                   <  2:32074‑M1/27‑1 (MQ=255)
                                                       aaCCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATa                                                                                           <  2:1173031/150‑1 (MQ=255)
                                                       aaCCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATa                                                                                           >  1:1216279/1‑150 (MQ=255)
                                                       aaCCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATa                                                                                           <  2:1174819/150‑1 (MQ=255)
                                                        aCCTTGGCTGCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATAc                                                                                          >  1:439330/1‑150 (MQ=255)
                                                          cTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCACCGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTg                                                                                        <  1:1234940/150‑1 (MQ=255)
                                                           ttGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACta                                                                                        <  1:763967/149‑2 (MQ=255)
                                                           ttGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACta                                                                                        >  2:763967/1‑148 (MQ=255)
                                                           ttGGCTCCAGGGTAATCACTGGGAGCTCCGACTCAACGCCGAGAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTGGAGTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGc                                                                                       <  2:851058/150‑1 (MQ=255)
                                                              gCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGaa                                                                                    <  2:968992/150‑1 (MQ=255)
                                                                tCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAata                                                                                  >  2:1337629/1‑150 (MQ=255)
                                                                tCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAata                                                                                  >  1:116262/1‑150 (MQ=255)
                                                                  cAGGGTAATCACTGGTAGCTCCGCATCCACTCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCAACGCCGCCACCAGTGCGGCGCTATCTTCGCTCCGGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATc                                                                                <  1:420827/150‑1 (MQ=255)
                                                                  cAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATc                                                                                <  1:1337629/150‑1 (MQ=255)
                                                                  cAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATc                                                                                <  2:1216279/150‑1 (MQ=255)
                                                                   gaacaggtgccagtgcgcgatattcgcagtattgccgaggccatcgcgaacaacgcgggcaagagtcaagaTACCGCCGCCACCAGTGCGGCGGTATc                                                                                                                                   <  2:180364‑M1/27‑1 (MQ=255)
                                                                   gaacaggtgccagtgcgcgatattcgcagtattgccgaggccatcgcgaacaacgcgggcaagagtcaagaTACCGCCGCCACCAGTGCGGCGGTATc                                                                                                                                   >  1:180364‑M1/72‑98 (MQ=255)
                                                                   aGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTg                                                                                        <  2:132969/141‑1 (MQ=255)
                                                                   aGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTg                                                                                        >  1:132969/1‑141 (MQ=255)
                                                                   aGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATcg                                                                               >  1:567551/1‑150 (MQ=255)
                                                                   aGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATcg                                                                               >  2:2279580/1‑150 (MQ=255)
                                                                   aGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATcg                                                                               >  2:39888/1‑150 (MQ=255)
                                                                     ggTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGGCTTGGCCTCGGCAATACTGCGAATATcgcg                                                                             <  1:1872877/150‑1 (MQ=255)
                                                                      gTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGGATATCgcgc                                                                            >  2:111142/1‑150 (MQ=255)
                                                                      gTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCgcgc                                                                            >  1:2075408/1‑150 (MQ=255)
                                                                            aCTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCgtt                                                                                                                 <  2:1439561/107‑1 (MQ=255)
                                                                            aCTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCgtt                                                                                                                 >  1:1439561/1‑107 (MQ=255)
                                                                              tGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCAc                                                                    <  2:768476/150‑1 (MQ=255)
                                                                              tGGTAGCTCCGACTCAACCCCGACAATCCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCCTGACTCTTGCCAGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCAc                                                                    <  1:49241/150‑1 (MQ=255)
                                                                                gTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCtt                                                                                                                                 <  1:630215/87‑1 (MQ=255)
                                                                                gTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCtt                                                                                                                                 >  2:630215/1‑87 (MQ=255)
                                                                                gTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCt                                                                  >  2:236907/1‑150 (MQ=255)
                                                                                gTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCt                                                                  >  1:1716464/1‑150 (MQ=255)
                                                                                    cTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTc                                                              <  1:1239653/150‑1 (MQ=255)
                                                                                           tCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATcgcg                                                                             <  2:534142/128‑1 (MQ=255)
                                                                                           tCAACGCCGACAATGCTTTGCACGATGGCGCGACACATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGacaaca                                                      <  2:1798293/150‑1 (MQ=255)
                                                                                                gCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAgg                                                    >  2:1277131/1‑148 (MQ=255)
                                                                                                gCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGc                                                  >  2:1554748/1‑150 (MQ=255)
                                                                                                   gACAATGCTTTGCACGATGGCGCGACACAGTCCGACGCTCACCGCCGCCACCAGTGCGGCGGTAGCTCCCCTGAGGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTg                                               <  1:419685/150‑1 (MQ=255)
                                                                                                    aCAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATcgc                                                                              <  2:1612707/118‑1 (MQ=255)
                                                                                                    aCAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATcgc                                                                              >  1:1612707/1‑118 (MQ=255)
                                                                                                      aaTGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTGCCGACAACAGGGCTTgcag                                            >  2:267927/1‑150 (MQ=255)
                                                                                                       aTGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTgcagc                                           >  2:1579671/1‑150 (MQ=255)
                                                                                                          cTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCgg                                                                                                <  2:1740346/94‑1 (MQ=255)
                                                                                                          cTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCAcc                                        >  1:78456/1‑150 (MQ=255)
                                                                                                          cTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCAcc                                        <  1:267927/150‑1 (MQ=255)
                                                                                                           tttGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGAc                                                                                                                              >  2:2234629/1‑63 (MQ=255)
                                                                                                           tttGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGAc                                                                                                                              <  1:2234629/63‑1 (MQ=255)
                                                                                                           tttGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCt                                       >  1:1009403/1‑150 (MQ=255)
                                                                                                              gCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGa                                    >  1:145718/1‑150 (MQ=255)
                                                                                                              gCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGa                                    <  1:2279580/150‑1 (MQ=255)
                                                                                                               cACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTcgc                                                                                                           <  1:956545/78‑1 (MQ=255)
                                                                                                               cACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTcgc                                                                                                           >  2:956545/1‑78 (MQ=255)
                                                                                                               cACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGag                                   <  2:1606463/150‑1 (MQ=255)
                                                                                                                aCGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCa                                                                     >  2:1245423/1‑115 (MQ=255)
                                                                                                                aCGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCa                                                                     <  1:1245423/115‑1 (MQ=255)
                                                                                                                  gATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTTACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGagcag                                >  1:2108375/1‑150 (MQ=255)
                                                                                                                  gATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGa                                                           >  1:835092/1‑123 (MQ=255)
                                                                                                                   aacaacgcgggcaagagtcaagaTACCGCCGCACTGGTGGCGGCGGTGTCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGc                                                  >  1:2282261‑M1/24‑131 (MQ=255)
                                                                                                                   aacaacgcgggcaagagtcaagaTACCGCCGCACTGGTGGCGGCGGTGTCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGc                                                  <  2:2282261‑M1/108‑1 (MQ=255)
                                                                                                                      aacgcgggcaagagtcaagaTACCGCCGCACTGGTGGCGGCGGTGTCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGcccc                            <  1:275269‑M1/130‑1 (MQ=255)
                                                                                                                          gACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACa                        >  2:1415785/1‑150 (MQ=255)
                                                                                                                          gACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACa                        >  2:1642670/1‑150 (MQ=255)
                                                                                                                            cacaATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAg                      <  2:1336028/150‑1 (MQ=255)
                                                                                                                             acaATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAgg                     >  2:1855628/1‑150 (MQ=255)
                                                                                                                             acaATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAgg                     >  1:1725273/1‑150 (MQ=255)
                                                                                                                               aaTCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGaa                   >  2:1605933/1‑150 (MQ=255)
                                                                                                                                 tCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAAt                 <  1:1642670/150‑1 (MQ=255)
                                                                                                                                 tCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAAt                 >  2:551650/1‑150 (MQ=255)
                                                                                                                                 tCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAAt                 <  2:2075408/150‑1 (MQ=255)
                                                                                                                                   cGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGa               <  2:567551/150‑1 (MQ=255)
                                                                                                                                    gACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGc                                                  >  1:1796937/1‑114 (MQ=255)
                                                                                                                                    gACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGc                                                  <  2:1796937/114‑1 (MQ=255)
                                                                                                                                      cgcgCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGGATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGacac            <  1:111142/150‑1 (MQ=255)
                                                                                                                                      aagaCACCGCCGCACTGGTGGCGGCGGTGTCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGaa                                                                                    >  1:1976286‑M1/5‑78 (MQ=255)
                                                                                                                                      aagaCACCGCCGCACTGGTGGCGGCGGTGTCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGaa                                                                                    <  2:1976286‑M1/74‑1 (MQ=255)
                                                                                                                                       gcgcACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACAcc           <  1:1521006/150‑1 (MQ=255)
                                                                                                                                          cACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGacaaca                                                      >  2:2401324/1‑104 (MQ=255)
                                                                                                                                          cACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGacaaca                                                      <  1:2401324/104‑1 (MQ=255)
                                                                                                                                             cgccgcCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGaca             >  2:2337742/1‑142 (MQ=255)
                                                                                                                                                cgccACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAg  <  2:1687433/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCGTCTTCCTGACCCTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAG  >  NC_002947/4934736‑4935029

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