Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,934,898 A→G 5.2% D579D (GAT→GAC flhA ← flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,934,8980AG 5.2% 360.0 / 5.4 115D579D (GAT→GACflhAflagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base A (53/56);  new base G (3/3);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 118 total aligned reads displayed.

TGACCCTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAGCTCTTCTGCAAGTTTAGG  >  NC_002947/4934746‑4935047
                                                                                                                                                         |                                                                                                                                                     
gtgccgggtgtcatttccttgtcggggctgctcaaggtgctccaagccctgttgtcggaacaggtgccagtgcgcgatattcgcagtattgccgaggccatcgcgaacaacgcgggcaagagtcaagaTACCGCC‑‑‑ACCAGTGCGGCGGTa                                                                                                                                                       >  1:3361‑M1/129‑150 (MQ=255)
   cccTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCATCGCCGCCACCAGTGCGGCGGTa                                                                                                                                                       <  1:2231483/150‑1 (MQ=255)
    ccTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTAt                                                                                                                                                      >  2:1992442/1‑150 (MQ=255)
       gCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCtt                                                                                                                                                   >  1:1817969/1‑150 (MQ=255)
         ccGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGa                                                                                                                                                 >  2:1239653/1‑150 (MQ=255)
         ccGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGa                                                                                                                                                 >  2:1694013/1‑150 (MQ=255)
              ccTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCtt                                                                                                                                            >  1:1085881/1‑150 (MQ=255)
               cTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTg                                                                                                                                           <  1:1694013/150‑1 (MQ=255)
                    cAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCcgc                                                                                                                                      <  1:2196754/150‑1 (MQ=255)
                             cAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGGGGTAGGTTGACTCTTGCCCGCGTTGTTcgc                                                                                                                             <  2:93202/150‑1 (MQ=255)
                               gCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGa                                                                                                                           <  2:1085881/150‑1 (MQ=255)
                                cAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATg                                                                                                                         >  1:1798293/1‑150 (MQ=255)
                                   ggtgctgcaagccctgttgtcggaacaggtgccagtgcgcgatattcgcagtattgccgaggccatcgcgaacaacgcgggcaagagtcaagaTACCGCCGCCACCAGTGCGGCGGTATc                                                                                                                                                     <  2:32074‑M1/27‑1 (MQ=255)
                                   ggtgctgcaagccctgttgtcggaacaggtgccagtgcgcgatattcgcagtattgccgaggccatcgcgaacaacgcgggcaagagtcaagaTACCGCCGCCACCAGTGCGGCGGTATc                                                                                                                                                     >  1:32074‑M1/94‑120 (MQ=255)
                                             aaCCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATa                                                                                                             <  2:1174819/150‑1 (MQ=255)
                                             aaCCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATa                                                                                                             <  2:1173031/150‑1 (MQ=255)
                                             aaCCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATa                                                                                                             >  1:1216279/1‑150 (MQ=255)
                                              aCCTTGGCTGCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATAc                                                                                                            >  1:439330/1‑150 (MQ=255)
                                                cTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCACCGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTg                                                                                                          <  1:1234940/150‑1 (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)
                                                        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                                                                                                 >  2:2279580/1‑150 (MQ=255)
                                                         aGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATcg                                                                                                 >  2:39888/1‑150 (MQ=255)
                                                         aGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATcg                                                                                                 >  1:567551/1‑150 (MQ=255)
                                                          gggTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATCGCCTCGGCAATACTGCCAATATcgc                                                                                                <  2:1356740/150‑1 (MQ=255)
                                                            gTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCgcgc                                                                                              >  1:2075408/1‑150 (MQ=255)
                                                                  aCTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCgtt                                                                                                                                   >  1:1439561/1‑107 (MQ=255)
                                                                    tGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCAc                                                                                      <  2:541295/150‑1 (MQ=255)
                                                                    tGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCAc                                                                                      <  1:788513/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                                                                                    >  1:1716464/1‑150 (MQ=255)
                                                                      gTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCt                                                                                    >  2:236907/1‑150 (MQ=255)
                                                                      gTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCt                                                                                    <  2:1817969/150‑1 (MQ=255)
                                                                                 tCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGATTATcgcg                                                                                               >  1:534142/1‑128 (MQ=255)
                                                                                 tCAACGCCGACAATGCTTTGCACGATGGCGCGACACATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGacaaca                                                                        <  2:1798293/150‑1 (MQ=255)
                                                                                   aaCGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAg                                                                       <  2:116262/150‑1 (MQ=255)
                                                                                      gCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAgg                                                                      <  1:1277131/148‑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                                                                                                >  1:1612707/1‑118 (MQ=255)
                                                                                          aCAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTgc                                                                <  1:39888/150‑1 (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:267927/150‑1 (MQ=255)
                                                                                                cTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCAcc                                                          >  1:78456/1‑150 (MQ=255)
                                                                                                 tttGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGAc                                                                                                                                                >  2:2234629/1‑63 (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)
                                                                                                        gATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGa                                                                             <  2:835092/123‑1 (MQ=255)
                                                                                                         aacaacgcgggcaagagtcaagaTACCGCCGCACTGGTGGCGGCGGTGTCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGc                                                                    <  2:2282261‑M1/108‑1 (MQ=255)
                                                                                                         aacaacgcgggcaagagtcaagaTACCGCCGCACTGGTGGCGGCGGTGTCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGc                                                                    >  1:2282261‑M1/24‑131 (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                                   >  2:551650/1‑150 (MQ=255)
                                                                                                                       tCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAAt                                   <  1:1642670/150‑1 (MQ=255)
                                                                                                                       gagtcaagaCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAAt                                   <  1:1828212‑M1/141‑1 (MQ=255)
                                                                                                                         cGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGa                                 <  2:567551/150‑1 (MQ=255)
                                                                                                                          gACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGc                                                                    <  2:1796937/114‑1 (MQ=255)
                                                                                                                          gACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGc                                                                    >  1:1796937/1‑114 (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)
                                                                                                                             gcgcACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCCCCTTGAGCAGCGCCGACCAGGAACTGACAcc                             >  2:1586073/1‑150 (MQ=255)
                                                                                                                                cACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGacaaca                                                                        <  1:2401324/104‑1 (MQ=255)
                                                                                                                                   cgccgcCACCAGTGCGGCGGTGTCTTGACTCTTGCCCGCGTTGTTCGCGATGGGCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCAc                       >  2:2392985/1‑150 (MQ=255)
                                                                                                                                   cgccgcCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGaca                               >  2:2337742/1‑142 (MQ=255)
                                                                                                                                   cgccgcCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGaca                               <  1:2337742/142‑1 (MQ=255)
                                                                                                                                      cgccACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAg                    <  2:1687433/150‑1 (MQ=255)
                                                                                                                                            caGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAAGACGTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAGCtcttc              <  1:1855628/150‑1 (MQ=255)
                                                                                                                                             aGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAGCtcttct             <  1:236907/150‑1 (MQ=255)
                                                                                                                                                gcggcgGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAGTATCGCGCACTGGCACCTGTTCCGACAACAGGGGTTGCGGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAGCTCTTCTGCa          >  2:275733/1‑150 (MQ=255)
                                                                                                                                                       gTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAggg                                                                     >  1:312692/1‑85 (MQ=255)
                                                                                                                                                       gTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAggg                                                                     <  2:312692/85‑1 (MQ=255)
                                                                                                                                                       gTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAGCTCTTCTGCAAGTTTa    >  2:540884/1‑150 (MQ=255)
                                                                                                                                                         aTCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAGCTCTTCTGCAAGTTTAgg  >  1:472071/1‑150 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
TGACCCTGCCCGGCCCTTTGCAGACTATTCAGCAAAATCTGTTCCAACCTTGGCTCCAGGGTAATCACTGGTAGCTCCGACTCAACGCCGACAATGCTTTGCACGATGGCGCGACACAATCCGACGCGCACCGCCGCCACCAGTGCGGCGGTATCTTGACTCTTGCCCGCGTTGTTCGCGATGGCCTCGGCAATACTGCGAATATCGCGCACTGGCACCTGTTCCGACAACAGGGCTTGCAGCACCTTGAGCAGCCCCGACAAGGAAATGACACCCGGCACCAGCTCTTCTGCAAGTTTAGG  >  NC_002947/4934746‑4935047

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: