Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,895,074 C→T 6.4% W10* (TGG→TAG)  yegH ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,895,0740CT 6.4% 233.1 / 5.0 78W10* (TGG→TAG) yegHputative transporter
Reads supporting (aligned to +/- strand):  ref base C (34/39);  new base T (2/3);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGGAGAACGAT  >  NC_002947/3894932‑3895222
                                                                                                                                              |                                                                                                                                                    
cATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGt                                                                                                                                               <  1:27422/150‑1 (MQ=255)
 atgatgAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTGGCTAGGcctgt                                                                                                                                              <  1:732537/150‑4 (MQ=255)
  tgatgaGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTgg                                                                                                                                             >  2:519909/1‑150 (MQ=255)
  tgatgaGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTgg                                                                                                                                             <  2:1271109/150‑1 (MQ=255)
           ccAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTa                                                                                                                                    <  1:141497/150‑1 (MQ=255)
               gCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGcca                                                                                                                                <  2:1208804/150‑1 (MQ=255)
                cTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATAAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATGAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccat                                                                                                                               >  2:1180416/1‑150 (MQ=255)
                 tCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccatt                                                                                                                              >  1:865609/1‑150 (MQ=255)
                  cAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattc                                                                                                                             <  2:720823/150‑1 (MQ=255)
                  cAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattc                                                                                                                             <  2:1323653/150‑1 (MQ=255)
                    ggCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTGGCTAGGcctgt                                                                                                                                              >  2:293486/1‑128 (MQ=255)
                    ggCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTGGCTAGGcctgt                                                                                                                                              <  1:293486/131‑4 (MQ=255)
                       cAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTa                                                                                                                        <  1:1631463/150‑1 (MQ=255)
                       cAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTa                                                                                                                        >  1:167813/1‑150 (MQ=255)
                       cAATGACCCGCGCGCGGTCGAGATGATGGGGCTGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTa                                                                                                                        <  2:1614297/150‑1 (MQ=255)
                         aTGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGt                                                                                                                      <  2:167813/150‑1 (MQ=255)
                           gACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTgg                                                                                                                    <  1:88436/150‑1 (MQ=255)
                                 cgcgcgGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTc                                                                                                              <  2:54461/150‑1 (MQ=255)
                                  gcgcgGTCGCGCTGATGGGGAGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATACGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCa                                                                                                             >  1:269317/1‑150 (MQ=255)
                                   cgcgGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCAc                                                                                                            <  1:591072/150‑1 (MQ=255)
                                              tGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATGAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGt                                                                                                 <  1:1180416/150‑1 (MQ=255)
                                                aTGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTca                                                                                               >  2:109751/1‑150 (MQ=255)
                                                  ggggCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGt                                                                                                                      <  1:1552728/125‑1 (MQ=255)
                                                  ggggCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGt                                                                                                                      >  2:1552728/1‑125 (MQ=255)
                                                       ggcAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattc                                                                                                                             >  1:320261/1‑113 (MQ=255)
                                                       ggcAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGccattc                                                                                                                             <  2:320261/113‑1 (MQ=255)
                                                          aaCTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACggg                                                                                     <  1:111269/150‑1 (MQ=255)
                                                                  ggCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGAt                                                                             <  1:519909/150‑1 (MQ=255)
                                                                   gCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGACTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATg                                                                            >  1:845626/1‑150 (MQ=255)
                                                                        gATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAg                                                                                                                                                    >  2:293700/1‑73 (MQ=255)
                                                                        gATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAg                                                                                                                                                    <  1:293700/73‑1 (MQ=255)
                                                                           ggCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                                                    <  1:4076/150‑1 (MQ=255)
                                                                              gATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTgg                                                                                                        <  1:1424154/111‑1 (MQ=255)
                                                                              gATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTgg                                                                                                        >  2:1424154/1‑111 (MQ=255)
                                                                                     acCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGACAGAAAGGCCGAAGAc                                                          >  1:1486288/1‑150 (MQ=255)
                                                                                     acCAGGTTGTCGATAACCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGAc                                                          >  1:1365944/1‑150 (MQ=255)
                                                                                      ccAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACt                                                         >  2:1377328/1‑150 (MQ=255)
                                                                                       cAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTg                                                        >  1:1555740/1‑150 (MQ=255)
                                                                                             gTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGc                                                  >  1:719986/1‑150 (MQ=255)
                                                                                             gTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGc                                                  >  2:1422903/1‑150 (MQ=255)
                                                                                              tCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCAAGACCGGCAGTACACTGGGGGGcc                                                 >  2:1583507/1‑149 (MQ=255)
                                                                                                gATGCCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGGCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTcc                                               >  2:884008/1‑150 (MQ=255)
                                                                                                 aTACCCAGCACCAGCTCGAGGGCGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACgggg                                                                                    <  1:373863/112‑1 (MQ=255)
                                                                                                 aTACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACgggg                                                                                    >  2:373863/1‑112 (MQ=255)
                                                                                                 aTACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                                                    <  2:788332/128‑1 (MQ=255)
                                                                                                 aTACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                                                    >  1:788332/1‑128 (MQ=255)
                                                                                                  tACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGGGGTATGGCCAGAAAGGCCGAAGACTGGGGGGCTCcgc                                             >  2:1620635/1‑150 (MQ=255)
                                                                                                     ccAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                                                    <  1:649868/124‑1 (MQ=255)
                                                                                                     ccAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGaaa                                                                    >  2:649868/1‑124 (MQ=255)
                                                                                                     ccAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGaa                                          >  1:1127233/1‑150 (MQ=255)
                                                                                                          aCCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACggg                                                                                     <  2:1237931/102‑1 (MQ=255)
                                                                                                          aCCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACggg                                                                                     >  1:1237931/1‑102 (MQ=255)
                                                                                                            cAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCCTTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTc                                   <  1:884008/150‑1 (MQ=255)
                                                                                                             aGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCa                                  <  1:1422903/150‑1 (MQ=255)
                                                                                                               cTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGTCGAAGACTGGGGGGCTCCGCGAAGGTGCTCata                                <  1:1620635/150‑1 (MQ=255)
                                                                                                                tCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCatat                               >  2:1462139/1‑150 (MQ=255)
                                                                                                                   aGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGa                                                           <  2:409828/119‑1 (MQ=255)
                                                                                                                   aGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGa                                                           >  1:409828/1‑119 (MQ=255)
                                                                                                                    ggACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGAc                           <  1:1462139/150‑1 (MQ=255)
                                                                                                                       cGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCtt                        <  2:1555740/150‑1 (MQ=255)
                                                                                                                        gATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTg                       <  2:538164/150‑1 (MQ=255)
                                                                                                                          tAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTc                                                                                                <  1:535697/75‑1 (MQ=255)
                                                                                                                          tAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTc                                                                                                >  2:535697/1‑75 (MQ=255)
                                                                                                                          tAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCata                                >  1:923514/1‑139 (MQ=255)
                                                                                                                          tAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCata                                <  2:923514/139‑1 (MQ=255)
                                                                                                                            aGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATaa                   >  2:390989/1‑150 (MQ=255)
                                                                                                                              cGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACa                 >  2:938850/1‑150 (MQ=255)
                                                                                                                                tCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGaa                                          >  1:888670/1‑123 (MQ=255)
                                                                                                                                tCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGaa                                          <  2:888670/123‑1 (MQ=255)
                                                                                                                                   cacggCCTGGCTAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGaa                                                                                 <  2:310032/78‑1 (MQ=255)
                                                                                                                                   cacggCCTGGCTAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGaa                                                                                 >  1:310032/4‑81 (MQ=255)
                                                                                                                                   aCAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAAttt            >  2:302377/1‑149 (MQ=255)
                                                                                                                                    cAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTgg           <  1:1377328/150‑1 (MQ=255)
                                                                                                                                      ggCCTAGCCAGGCCGTGGGGTCGACTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGga         <  2:845626/150‑1 (MQ=255)
                                                                                                                                          tAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCtt                        >  2:19735/1‑131 (MQ=255)
                                                                                                                                          tAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCtt                        <  1:19735/131‑1 (MQ=255)
                                                                                                                                          tAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGGAGAAc     <  1:938850/150‑1 (MQ=255)
                                                                                                                                          tAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGGAGAAc     <  1:109751/150‑1 (MQ=255)
                                                                                                                                             ccAGGCCGTGGGGTCGGCTAGCCATTCCATGAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGGAGAACGAt  >  1:300287/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CATGATGAGCGCCAGGCTCAGGCCAATGACCCGCGCGCGGTCGCGCTGATGGGGCGGCAACTTGTCGGCAAGGATGGCGATGAACACCAGGTTGTCGATACCCAGCACCAGCTCGAGGACGATAAGCGTCAACAGGCCTAGCCAGGCCGTGGGGTCGGCTAGCCATTCCATTAGTGGGTACTCACGTGGGACACGTCACAAGGACGGGGGAAGGATGGCCAGAAAGGCCGAAGACTGGGGGGCTCCGCGAAGGTGCTCATATAGACCTTGATAACAAATTGGGAGAACGAT  >  NC_002947/3894932‑3895222

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