Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,478,807 T→C 5.4% N329D (AAC→GAC) ‡ proA ← gamma‑glutamyl phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,478,8070TC 5.4% 1271.0 / 35.2 393N329D (AAC→GAC) ‡proAgamma‑glutamyl phosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (183/187);  new base C (5/16);  total (190/203)
Fisher's exact test for biased strand distribution p-value = 2.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 399 total aligned reads displayed.

GAAGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATTT  >  NC_002947/5478658‑5478956
                                                                                                                                                     |                                                                                                                                                     
gAAGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGtt                                                                                                                                                       <  1:4782967/150‑1 (MQ=255)
 aaGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTc                                                                                                                                                      >  1:8381291/1‑150 (MQ=255)
 aaGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTc                                                                                                                                                      <  2:7436908/150‑1 (MQ=255)
  aGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCa                                                                                                                                                     <  1:3958450/150‑1 (MQ=255)
      aTCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGcc                                                                                                                                                 >  2:1621758/1‑150 (MQ=255)
      aTCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGcc                                                                                                                                                 >  2:2054384/1‑150 (MQ=255)
       tCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCg                                                                                                                                                >  2:1959160/1‑150 (MQ=255)
       tCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCg                                                                                                                                                <  2:6220120/150‑1 (MQ=255)
         ggCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGGTCGACGGCCTGAACCAGGCCATc                                                                                                                                              <  2:4299460/150‑1 (MQ=255)
         ggCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATCGCCTGGTTCAGGCCGTc                                                                                                                                              <  1:8676659/150‑1 (MQ=255)
         ggCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGACGGCCTGAACCAGGCCATc                                                                                                                                              <  2:7628825/150‑1 (MQ=255)
            aaaGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGAc                                                                                                                                           >  1:1085933/1‑150 (MQ=255)
            aaaGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGAc                                                                                                                                           <  2:2810514/150‑1 (MQ=255)
               gCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCAc                                                                                                                                        >  1:2393503/1‑150 (MQ=255)
                   gTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGa                                                                                                                                    >  1:1350601/1‑150 (MQ=255)
                    tCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGaa                                                                                                                                   <  1:6647387/150‑1 (MQ=255)
                     cGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGGTCGACGGCCTGAACCAGGCCATCGAGCAcatcaac                                                                                                                                  <  1:2580536/150‑7 (MQ=38)
                      ggggTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATc                                                                                                                                 <  1:4492129/150‑1 (MQ=255)
                        ggTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGa                                                                                                                               <  2:6268526/150‑1 (MQ=255)
                         gtgtTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGAc                                                                                                                              <  2:7240707/150‑1 (MQ=255)
                          tgtTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACa                                                                                                                             <  2:5115989/150‑1 (MQ=255)
                           gtTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAg                                                                                                                            >  2:1803629/1‑150 (MQ=255)
                             tGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGa                                                                                                                          <  1:5592235/150‑1 (MQ=255)
                               aGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATc                                                                                                                        <  2:6678453/150‑1 (MQ=255)
                                     aCCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCg                                                                                                                  >  2:4510834/1‑150 (MQ=255)
                                       gaagcaggAAGCGCCCTCTGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAg                                                                                                                >  2:8228637/8‑150 (MQ=255)
                                       cGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAg                                                                                                                <  2:1149045/150‑1 (MQ=255)
                                       cGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAg                                                                                                                <  2:7011350/150‑1 (MQ=255)
                                       cGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAg                                                                                                                <  1:643194/150‑1 (MQ=255)
                                        gACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAgg                                                                                                               <  1:5819555/150‑1 (MQ=255)
                                         aCGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGt                                                                                                              <  2:3693078/150‑1 (MQ=255)
                                             cGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTc                                                                                                          >  2:7115277/1‑150 (MQ=255)
                                             cGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTc                                                                                                          <  2:1716132/150‑1 (MQ=255)
                                                aTCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCgt                                                                                                       <  2:2866298/150‑1 (MQ=255)
                                                   gACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGtgtg                                                                                                    >  1:7876098/1‑150 (MQ=255)
                                                    aCCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGc                                                                                                   >  1:6219923/1‑150 (MQ=255)
                                                      cTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCa                                                                                                 >  2:4017552/1‑150 (MQ=255)
                                                      cTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCa                                                                                                 <  1:525571/150‑1 (MQ=255)
                                                       tCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAg                                                                                                >  1:1778706/1‑150 (MQ=255)
                                                       tCAGCCATGAACTGACGGGCTACACCCTGGTGTTCGCTGATGATCGAGTAGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAg                                                                                                >  2:2171655/1‑150 (MQ=255)
                                                         aGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTc                                                                                              >  1:7520298/1‑150 (MQ=255)
                                                          gCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCg                                                                                             <  1:3979565/150‑1 (MQ=255)
                                                             aTGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCt                                                                                          >  1:487369/1‑150 (MQ=255)
                                                                aaCTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCg                                                                                       <  2:1221718/150‑1 (MQ=255)
                                                                aaCTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCg                                                                                       >  2:3531043/1‑150 (MQ=255)
                                                                 aCTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCgg                                                                                      <  1:7153145/150‑1 (MQ=255)
                                                                   tGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTg                                                                                    <  1:6539503/150‑1 (MQ=255)
                                                                    gACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACa                                                                                                                             >  1:6902232/1‑108 (MQ=255)
                                                                         gCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGt                                                                              >  1:2697807/1‑150 (MQ=255)
                                                                          cTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGcctgac                                                                             <  2:3152129/150‑5 (MQ=255)
                                                                          cTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGtt                                                                             >  2:1102883/1‑150 (MQ=255)
                                                                            tCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTg                                                                           <  1:1633314/150‑1 (MQ=255)
                                                                                ccTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGgcgc                                                                       >  2:6161540/1‑150 (MQ=255)
                                                                                ccTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGgcgc                                                                       >  1:3447682/1‑150 (MQ=255)
                                                                                ccTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGgcgc                                                                       >  2:8148163/1‑150 (MQ=255)
                                                                                 cTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCt                                                                      <  2:4827974/150‑1 (MQ=255)
                                                                                   ggTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtga                                                                    <  2:591247/150‑1 (MQ=255)
                                                                                    gtgtTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtgat                                                                   <  1:5731914/150‑1 (MQ=255)
                                                                                      gtTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtgatga                                                                 <  2:5720336/150‑1 (MQ=255)
                                                                                       gtCGCTGAGGATCGAGTCGGTGTGGGGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTgatgat                                                                <  2:2437963/149‑1 (MQ=255)
                                                                                        tCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCGTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGTTGATa                                                               >  2:261385/1‑150 (MQ=255)
                                                                                          gCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGc                                                             <  2:1306253/150‑1 (MQ=255)
                                                                                          gCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGc                                                             >  2:5339122/1‑150 (MQ=255)
                                                                                           cTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGcc                                                            <  1:5339122/150‑1 (MQ=255)
                                                                                             gatgatCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTg                                                          >  1:6008877/1‑150 (MQ=255)
                                                                                              atgatCGCGTCGGGGGGGTGCGCGCCAGAGTGGTTGATGTGCTCGCTGGCCTGGGTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTgg                                                         <  2:1334806/150‑1 (MQ=255)
                                                                                              atgatCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTgg                                                         >  1:6795836/1‑150 (MQ=255)
                                                                                              atgatCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTgg                                                         <  1:1015694/150‑1 (MQ=255)
                                                                                              atgatCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTgg                                                         <  1:943173/150‑1 (MQ=255)
                                                                                              atgatCGAGTCGGTGTGGTGCGAGCCAGAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTgg                                                         <  2:5196998/150‑1 (MQ=255)
                                                                                                gatCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGt                                                       <  1:4141664/150‑1 (MQ=255)
                                                                                                  tCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTAc                                                     >  1:4296094/1‑150 (MQ=255)
                                                                                                    gAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGc                                                   >  1:2016042/1‑150 (MQ=255)
                                                                                                      gTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTc                                                 <  2:395513/150‑1 (MQ=255)
                                                                                                      gTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTc                                                 <  2:680719/150‑1 (MQ=255)
                                                                                                       tCGGTGTGGTGCGAGCCATAGTAGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCg                                                >  1:3322776/1‑150 (MQ=255)
                                                                                                        cGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGTGCTGATGATAGCCTGGGTACGCTCGc                                               <  2:5771670/150‑1 (MQ=255)
                                                                                                        cGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGc                                               <  2:8434016/150‑1 (MQ=255)
                                                                                                         ggTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCa                                              >  1:8581439/1‑150 (MQ=255)
                                                                                                              ggtgCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACcgc                                         >  2:8592706/1‑150 (MQ=255)
                                                                                                                 gCGAGCCATAGTGGTCGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCa                                      <  1:6135190/150‑1 (MQ=255)
                                                                                                                     gCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTc                                  >  1:3833560/1‑150 (MQ=255)
                                                                                                                     gCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTc                                  <  2:3303200/150‑1 (MQ=255)
                                                                                                                        aTAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCAc                               <  2:6469510/150‑1 (MQ=255)
                                                                                                                         tAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACg                              >  2:905295/1‑150 (MQ=255)
                                                                                                                              gTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCtt                         >  2:4699404/1‑150 (MQ=255)
                                                                                                                               ttGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTtc                        <  1:7267532/150‑1 (MQ=255)
                                                                                                                                    gtgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCAc                   >  1:1313432/1‑150 (MQ=255)
                                                                                                                                        tCGACGGCCTGAACCAGGCCATCGAGCACGCGAATCGACAGGATCGCCTCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAg               >  2:1549241/1‑150 (MQ=255)
                                                                                                                                        tCGACGGCCTGAACCAGGCCATCGAGCACGCGAATCGACAGGATCGCATCGAGGTCTTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAg               >  2:5028245/1‑150 (MQ=255)
                                                                                                                                          gATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAgcg             >  1:3636753/1‑150 (MQ=255)
                                                                                                                                          gATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAgcg             >  2:750921/1‑150 (MQ=255)
                                                                                                                                            tGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAgcggc           <  1:1959160/150‑1 (MQ=255)
                                                                                                                                             ggCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAgcggcg          <  1:2918757/150‑1 (MQ=255)
                                                                                                                                               ccTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCggg        >  1:7928083/1‑150 (MQ=255)
                                                                                                                                               ccTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCggg        >  2:798224/1‑150 (MQ=255)
                                                                                                                                                cTGGTTCAGGCCGTCGACCACGCGAATCGACATGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGc       >  1:5924057/1‑150 (MQ=255)
                                                                                                                                                cTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGc       >  2:2272239/1‑150 (MQ=255)
                                                                                                                                                 tGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGcc      >  2:7199202/1‑150 (MQ=255)
                                                                                                                                                   gTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCAt    <  2:6900523/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATTT  >  NC_002947/5478658‑5478956

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: