Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,509,506 A→T 10.9% S355C (AGC→TGC)  urtA → urea ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,509,5060AT10.9% 234.7 / 18.8 83S355C (AGC→TGC) urtAurea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (32/41);  new base T (6/3);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

AAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGGCGCAGCCATGGAGC  >  NC_002947/5509359‑5509655
                                                                                                                                                   |                                                                                                                                                     
aaggccaaggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGc                                                                                                                                                     >  2:422070/1‑150 (MQ=255)
   gccaaggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCTCTGCCCGGCCAGTGCTTc                                                                                                                                                  <  2:388676/150‑1 (MQ=255)
      aaggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAg                                                                                                                                               <  1:1195059/150‑1 (MQ=255)
      aaggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAg                                                                                                                                               <  1:422070/150‑1 (MQ=255)
                cAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCgg                                                                                                                                     <  2:190976/150‑1 (MQ=255)
                    tGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCtt                                                                                                                                 >  2:1297243/1‑150 (MQ=255)
                             ggCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCTCTGCCCGGCCAGTGCtt                                                                                                                                                   <  2:775600/123‑1 (MQ=255)
                             ggCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCTCTGCCCGGCCAGTGCtt                                                                                                                                                   >  1:775600/1‑123 (MQ=255)
                                cGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCAt                                                                                                                     >  1:266414/1‑150 (MQ=255)
                                 gTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATg                                                                                                                    >  1:792242/1‑150 (MQ=255)
                                       aaCGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAg                                                                                                              >  1:369348/1‑150 (MQ=255)
                                        aCGACCCCATGGAAGCCACCTACGTGGGCATCCATATCTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCCCTCAGTTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGa                                                                                                             <  1:1384679/150‑1 (MQ=255)
                                             cccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaac                                                                                                        >  1:1174024/1‑150 (MQ=255)
                                             cccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaac                                                                                                        <  2:1324834/150‑1 (MQ=255)
                                              ccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCAt                                                                                                                     <  1:146644/136‑1 (MQ=255)
                                              ccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCAt                                                                                                                     >  2:146644/1‑136 (MQ=255)
                                              ccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaacc                                                                                                       <  2:832208/150‑1 (MQ=255)
                                              ccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaac                                                                                                        >  2:858080/1‑149 (MQ=255)
                                              ccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaac                                                                                                        <  1:858080/149‑1 (MQ=255)
                                                  ggAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATga                                                                                                                   <  1:1399215/134‑2 (MQ=255)
                                                  ggAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATga                                                                                                                   >  2:1399215/1‑133 (MQ=255)
                                                    aatccACCTACGTGGGCACCCACATGTGGGCCCAGCCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAACCTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAaccaccacc                                                                                                 <  1:30939/147‑1 (MQ=255)
                                                      gCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTg                                                                                               >  2:712281/1‑150 (MQ=255)
                                                         aCCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGcac                                                                                            >  2:436959/1‑150 (MQ=255)
                                                         aCCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGcac                                                                                            >  1:697089/1‑150 (MQ=255)
                                                                 gggCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGt                                                                                    <  1:617648/150‑1 (MQ=255)
                                                                   gCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGtga                                                                                  <  1:1290421/150‑1 (MQ=255)
                                                                      tCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGtgatga                                                                               >  1:1066842/1‑150 (MQ=255)
                                                                         acaTGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAg                                                                                                                                               >  2:500210/1‑83 (MQ=255)
                                                                         acaTGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAg                                                                                                                                               <  1:500210/83‑1 (MQ=255)
                                                                         acaTGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCg                                                                            >  2:1357095/1‑150 (MQ=255)
                                                                          caTGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCgg                                                                           <  1:712281/150‑1 (MQ=255)
                                                                           aTGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGc                                                                          >  2:879253/1‑150 (MQ=255)
                                                                           aTGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGc                                                                          <  2:369348/150‑1 (MQ=255)
                                                                             gtgGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCGGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCga                                                                        <  2:358991/150‑1 (MQ=255)
                                                                               gggCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCgaga                                                                      <  1:436959/150‑1 (MQ=255)
                                                                                  cccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTc                                                                   <  2:1066842/150‑1 (MQ=255)
                                                                                   ccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCa                                                                  >  2:936049/1‑150 (MQ=255)
                                                                                   ccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCa                                                                  <  1:879253/150‑1 (MQ=255)
                                                                                   ccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCa                                                                  <  2:803659/150‑1 (MQ=255)
                                                                                      ggCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGtga                                                                                  >  1:698880/1‑131 (MQ=255)
                                                                                      ggCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGtga                                                                                  <  2:698880/131‑1 (MQ=255)
                                                                                      ggCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCgg                                                                           <  2:1185286/138‑1 (MQ=255)
                                                                                      ggCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCgg                                                                           >  1:1185286/1‑138 (MQ=255)
                                                                                          gTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGAc                                                           >  1:1334801/1‑150 (MQ=255)
                                                                                             gagaAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGc                                                        <  1:942724/150‑1 (MQ=255)
                                                                                             gagaAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGc                                                        <  1:353547/150‑1 (MQ=255)
                                                                                             gagaAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGc                                                        >  2:393368/1‑150 (MQ=255)
                                                                                              agaCAGCCAAGTCCACTGACGTTGACAACGTCCGCGAAGCACTGGCCGGGCAGAGCGTCGAGGCGCCGTCGGGCTTCACCCTCACCATGGACAAGACCCACCACCACCTGCACCACCCGGTGATGATCGGCGAGGCTCAGGATGACGGcc                                                       >  2:68114/1‑150 (MQ=255)
                                                                                                    ccAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCa                                                 >  1:879543/1‑150 (MQ=255)
                                                                                                      aaGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGc                                               <  2:1174024/150‑1 (MQ=255)
                                                                                                       aGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCg                                              <  2:879543/150‑1 (MQ=255)
                                                                                                         tCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCgtg                                            <  1:640416/150‑1 (MQ=255)
                                                                                                            aCTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAACGCTTAGAAGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGAGGACGGCCAGTTCAGCgtggtg                                         >  1:958005/1‑150 (MQ=255)
                                                                                                              tGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGtgtg                                       >  1:560394/1‑150 (MQ=255)
                                                                                                                aCGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGCCCAACCACCACCTGCACAAGCCGGtg                                                                                   >  1:601843/1‑104 (MQ=255)
                                                                                                                aCGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGtg                                                                                   <  2:601843/104‑1 (MQ=255)
                                                                                                                  gTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGac                                                                                                            >  2:808423/1‑77 (MQ=255)
                                                                                                                  gTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGac                                                                                                            <  1:808423/77‑1 (MQ=255)
                                                                                                                   ttGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGaaa                                  <  2:266414/150‑1 (MQ=255)
                                                                                                                      aCAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCg                               >  2:1313670/1‑150 (MQ=255)
                                                                                                                       cAAGGCCCGCGAAGCACTGGCCGGGCAGAGCGTCGAGTCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACATGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGa                              <  1:224027/150‑1 (MQ=255)
                                                                                                                        aaGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGag                             >  2:982295/1‑150 (MQ=255)
                                                                                                                        aaGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGag                             <  2:1334801/150‑1 (MQ=255)
                                                                                                                        aaGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGACGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGag                             >  2:1261351/1‑150 (MQ=255)
                                                                                                                        aaGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCAACTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGag                             <  2:697089/150‑1 (MQ=255)
                                                                                                                            gcctcgAAGCTCTGCCCGGCCAGTGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGagcagc                         >  1:1272461/5‑150 (MQ=255)
                                                                                                                              ctcgAAGCTCTGCCCGGCCAGTGCTTCGCGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGTAGCCg                       >  2:568193/3‑150 (MQ=255)
                                                                                                                              ctcgAAGCTCTGCCCGGCCAGTGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCg                       >  2:281991/3‑150 (MQ=255)
                                                                                                                                  aaGCTCTGCCCGGCCAGTGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGc                                                        >  1:366197/1‑113 (MQ=255)
                                                                                                                                  aaGCTCTGCCCGGCCAGTGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGc                                                        <  2:366197/113‑1 (MQ=255)
                                                                                                                                  aaGCTCTGCCCGGCCAGTGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGc                   >  2:544253/1‑150 (MQ=255)
                                                                                                                                   aGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACCGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCg                  <  1:1357095/150‑1 (MQ=255)
                                                                                                                                    gCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCgg                 <  2:958005/150‑1 (MQ=255)
                                                                                                                                     cACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCga                                                                        >  2:1334554/1‑94 (MQ=255)
                                                                                                                                     cACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCga                                                                        <  1:1334554/94‑1 (MQ=255)
                                                                                                                                      aCTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGgc               >  1:92294/1‑150 (MQ=255)
                                                                                                                                      aCTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGgc               <  2:560394/150‑1 (MQ=255)
                                                                                                                                         ggCCGGGCAGAGCGTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGGCGCa            <  1:509194/150‑1 (MQ=255)
                                                                                                                                               gCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCga                                                                        >  2:834465/1‑84 (MQ=255)
                                                                                                                                               gCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCga                                                                        <  1:834465/84‑1 (MQ=255)
                                                                                                                                               gCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGGCGCAGCCATg      <  2:92294/150‑1 (MQ=255)
                                                                                                                                                cAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGGCGCAGCCATgg     <  1:577227/150‑1 (MQ=255)
                                                                                                                                                   agCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGGCGCAGCCATGGAGc  >  1:526800/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGGCGCAGCCATGGAGC  >  NC_002947/5509359‑5509655

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: