Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,175,418 G→A 6.4% L390L (CTG→TTG)  arsB‑I ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,175,4180GA 6.4% 221.5 / 3.5 77L390L (CTG→TTG) arsB‑Iarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base A (2/3);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

ATGAACAAGACTTTCATCAGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGATTGGATCGACAAGGCCCCTACCAGCA  >  NC_002947/2175271‑2175564
                                                                                                                                                   |                                                                                                                                                  
aTGAACAAGACTTTCATCAGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCa                                                                                                                                                  >  2:299466/1‑150 (MQ=255)
 tGAACAAGACTTTCATCAGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCAc                                                                                                                                                 <  1:825433/150‑1 (MQ=255)
  gAACAAGACTTTCATCAGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACg                                                                                                                                                >  2:1006920/1‑150 (MQ=255)
   aaCAAGACTTTCATCAGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGt                                                                                                                                               <  1:1776874/150‑1 (MQ=255)
     cAAGACTTTCATCAGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGc                                                                                                                                             >  2:1319098/1‑150 (MQ=255)
        gACTTTCATCAGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGc                                                                                                                                          >  1:197135/1‑150 (MQ=255)
               atcaGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCa                                                                                                                                   >  1:1086119/1‑150 (MQ=255)
                  aGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGt                                                                                                                                               >  2:906271/1‑135 (MQ=255)
                  aGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGt                                                                                                                                               <  1:906271/135‑1 (MQ=255)
                        ctcCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCa                                                                                                                          <  2:1435886/150‑1 (MQ=255)
                         tcCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTGGCCTCAGATGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAg                                                                                                                         >  2:751042/1‑150 (MQ=255)
                         tcCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAg                                                                                                                         >  1:678946/1‑150 (MQ=255)
                         tcCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAg                                                                                                                         >  1:503376/1‑150 (MQ=255)
                           cGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGc                                                                                                                       >  2:1882001/1‑150 (MQ=255)
                            ggTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCt                                                                                                                      <  1:1882001/150‑1 (MQ=255)
                            ggTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGGTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGAAAAGTGGCCAggat                                                                                                                      >  2:1808673/1‑148 (MQ=255)
                                  aaaGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAACAAAGTGGCCCGGCTGCCGag                                                                                                                >  2:394574/1‑149 (MQ=255)
                                   aaGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAACACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGAtt                                                                                                               >  2:1562786/1‑150 (MQ=255)
                                    aGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTg                                                                                                              <  2:678946/150‑1 (MQ=255)
                                    aGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTg                                                                                                              <  1:1638538/150‑1 (MQ=255)
                                    aGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTg                                                                                                              <  2:1086119/150‑1 (MQ=255)
                                      cTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTggg                                                                                                            <  2:1040988/150‑1 (MQ=255)
                                      cTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTggg                                                                                                            <  1:299466/150‑1 (MQ=255)
                                      cTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTggg                                                                                                            <  2:880196/150‑1 (MQ=255)
                                                tGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAg                                                                                                                                                    >  1:864477/1‑100 (MQ=255)
                                                tGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAg                                                                                                                                                    <  2:864477/100‑1 (MQ=255)
                                                  gCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGCAGTAGCCCACCGTGATGACCCTGCACTTCCTTGGCAAAGCGGTCGTGAATAGCAAAGTGGCCAGGCTGCCTATTGGGGTGATTTTCGGa                                                                                                <  2:1191114/150‑1 (MQ=255)
                                                   ccAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGAc                                                                                               >  2:1892956/1‑150 (MQ=255)
                                                   ccAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGAc                                                                                               >  2:699245/1‑150 (MQ=255)
                                                    cAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGAcc                                                                                              >  2:133070/1‑150 (MQ=255)
                                                    cAATGCGGCCAAGGTGGCAAGAAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGAcc                                                                                              <  2:410320/150‑1 (MQ=255)
                                                     aaTGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAACACGt                                                                                                                                               >  1:663366/1‑100 (MQ=255)
                                                     aaTGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAACACGt                                                                                                                                               <  2:663366/100‑1 (MQ=255)
                                                                 gTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTgg                                                                                                                             >  1:1352975/1‑106 (MQ=255)
                                                                 gTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTgg                                                                                                                             <  2:1352975/106‑1 (MQ=255)
                                                                      aaGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAg                                                                                                                                           >  2:1790310/1‑87 (MQ=255)
                                                                      aaGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGAc                                                                            >  2:1317987/1‑150 (MQ=255)
                                                                      aaGCAAGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAg                                                                                                                                           <  1:1790310/87‑1 (MQ=255)
                                                                          ccgATCTGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCACAGGTGCTGACTATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTg                                                                        <  2:1508345/148‑1 (MQ=255)
                                                                          aGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTCGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGAGTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTg                                                                        <  1:1850393/150‑1 (MQ=255)
                                                                           ggATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTgg                                                                       <  1:1317987/150‑1 (MQ=255)
                                                                                 gCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGa                                                                 >  2:938820/1‑150 (MQ=255)
                                                                                 gCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCAAGTGCAGCCACAGCAAAGTGGCCAGGATGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGa                                                                 >  2:123290/1‑150 (MQ=255)
                                                                                  cAGCGTTAGGACGATGCCGGTACGGGAGTAGGAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTggagatcgg                                                                >  1:1508345/1‑143 (MQ=255)
                                                                                       ttAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCa                                                                                                                                            <  2:1057518/69‑1 (MQ=255)
                                                                                       ttAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCa                                                                                                                                            >  1:1057518/1‑69 (MQ=255)
                                                                                       ttAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATgg                                                           >  2:1201820/1‑150 (MQ=255)
                                                                                       ttAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGGGATGACGATGTTCTGGCGGGCCAGCACGTGCAGACACAGCAAAGTGGCCAGGCAGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGCTGACATTGGCGTAGAccggcg                                                           >  2:503071/1‑146 (MQ=255)
                                                                                           gaggATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGc                                                                                                                                    >  1:1939956/4‑73 (MQ=255)
                                                                                           gACGATGCCGGTCCGGAAGTAGTAGCCCCAGGTGATGGCGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTc                                                       <  1:503071/150‑1 (MQ=255)
                                                                                           gACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGc                                                                                                                                    <  2:1939956/73‑1 (MQ=255)
                                                                                           gACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTc                                                       <  2:197135/150‑1 (MQ=255)
                                                                                           gACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGAATTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGCGCGTGCTGCAGATTGGGGTGATTTTCGGACAGAGGTGACAAACGATGACGTTGGGGTAGACCAGGGCCTc                                                       >  1:1575636/1‑150 (MQ=255)
                                                                                             cGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTcgc                                                     <  1:642823/150‑1 (MQ=255)
                                                                                                  ccGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGCCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGAACg                                                >  1:632044/1‑150 (MQ=255)
                                                                                                   cGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGa                                               <  1:133070/150‑1 (MQ=255)
                                                                                                        cGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCg                                          <  2:413873/150‑1 (MQ=255)
                                                                                                          aaaGTAATAGCCCCAGGTGATTACGATGTTCTTGCGTGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCGGCCAATTGGGGTGATTTTCGGACCGAGGTCACATCCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGt                                        <  2:632044/149‑1 (MQ=255)
                                                                                                           aaGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTg                                       <  2:503376/150‑1 (MQ=255)
                                                                                                                gtagCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCgg                                                                                                 >  1:412409/1‑87 (MQ=255)
                                                                                                                gtagCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCgg                                                                                                 <  2:412409/87‑1 (MQ=255)
                                                                                                                       cAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGAt                           <  1:938820/150‑1 (MQ=255)
                                                                                                                       cAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGAt                           <  1:1481404/150‑1 (MQ=255)
                                                                                                                        aGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGAtt                          <  1:699245/150‑1 (MQ=255)
                                                                                                                        aGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCCGTGCCCTCGGCGGAtc                          <  2:568027/150‑2 (MQ=255)
                                                                                                                          gTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGt                                                                                         >  2:859318/1‑85 (MQ=255)
                                                                                                                          gTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGt                                                                                         <  1:859318/85‑1 (MQ=255)
                                                                                                                             atgaCGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTg                                                                        >  2:612576/1‑99 (MQ=255)
                                                                                                                             atgaCGATGTTCTTGCGCACCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTg                                                                        <  1:612576/99‑1 (MQ=255)
                                                                                                                                aCGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGATTGGATCGAc                  >  2:1510092/1‑150 (MQ=255)
                                                                                                                                  gATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGATTGGATCGACaa                <  2:1785786/150‑1 (MQ=255)
                                                                                                                                     gTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACAGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGATTGGATGGACAAGGc             <  2:1575636/150‑1 (MQ=255)
                                                                                                                                      ttCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGATTGGATCGACAAGGcc            <  2:1606101/150‑1 (MQ=255)
                                                                                                                                        cTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGATTGGATCGACAAGGcccc          >  2:717473/1‑150 (MQ=255)
                                                                                                                                         ttGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTg                                                                        <  1:1659156/87‑1 (MQ=255)
                                                                                                                                         ttGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTg                                                                        >  2:1659156/1‑87 (MQ=255)
                                                                                                                                         ttGCGGGCCAACACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGATTGGATCGACAAGGCCCCt         <  1:1562786/150‑1 (MQ=255)
                                                                                                                                             gggCCAGCACGTCCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATGTTCGTACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGAAGCCGGTGGCCTCGGCGGATTGGATCGACAAGGCCCCTACCa     >  1:636353/1‑150 (MQ=255)
                                                                                                                                                ccAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGATTGGATCGACAAGGCCCCTACCAGca  >  2:283874/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATGAACAAGACTTTCATCAGGTATCTCCGGTTCAAAAGCTCAGGCGGATGGCCAATGCGGCCAAGGTGGCAAGCAGGATCGGCAGCGTTAGGACGATGCCGGTACGGAAGTAGTAGCCCCAGGTGATGACGATGTTCTTGCGGGCCAGCACGTGCAGCCACAGCAAAGTGGCCAGGCTGCCGATTGGGGTGATTTTCGGACCGAGGTCACAACCGATGACGTTGGCGTAGACCATGGCCTCGCGCACGATGCCGGTGGCCTCGGCGGATTGGATCGACAAGGCCCCTACCAGCA  >  NC_002947/2175271‑2175564

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