Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,509,502 G→C 11.0% G353G (GGG→GGC 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,5020GC11.0% 250.8 / 17.8 82G353G (GGG→GGCurtAurea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base C (6/3);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TACGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGGCGCAGCCATG  >  NC_002947/5509353‑5509651
                                                                                                                                                     |                                                                                                                                                     
tACGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCggg                                                                                                                                                       >  1:845571/1‑150 (MQ=255)
  cGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCCCCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCa                                                                                                                                                     >  2:617648/1‑150 (MQ=255)
  cGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCa                                                                                                                                                     <  2:1337930/150‑1 (MQ=255)
  cGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCa                                                                                                                                                     <  2:1337279/150‑1 (MQ=255)
      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                                                                                                    <  2:1324834/150‑1 (MQ=255)
                                                   cccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaac                                                                                                    >  1:1174024/1‑150 (MQ=255)
                                                    ccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCAt                                                                                                                 >  2:146644/1‑136 (MQ=255)
                                                    ccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCAt                                                                                                                 <  1:146644/136‑1 (MQ=255)
                                                    ccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaacc                                                                                                   <  2:832208/150‑1 (MQ=255)
                                                    ccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaac                                                                                                    >  2:858080/1‑149 (MQ=255)
                                                    ccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaac                                                                                                    <  1:858080/149‑1 (MQ=255)
                                                        ggAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATga                                                                                                               >  2:1399215/1‑133 (MQ=255)
                                                        ggAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATga                                                                                                               <  1:1399215/134‑2 (MQ=255)
                                                          aatccACCTACGTGGGCACCCACATGTGGGCCCAGCCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAACCTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAaccaccacc                                                                                             <  1:30939/147‑1 (MQ=255)
                                                            gCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTg                                                                                           >  2:712281/1‑150 (MQ=255)
                                                               aCCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGcac                                                                                        >  1:697089/1‑150 (MQ=255)
                                                               aCCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGcac                                                                                        >  2:436959/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:369348/150‑1 (MQ=255)
                                                                                 aTGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGc                                                                      >  2:879253/1‑150 (MQ=255)
                                                                                   gtgGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCGGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCga                                                                    <  2:358991/150‑1 (MQ=255)
                                                                                     gggCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCgaga                                                                  <  1:436959/150‑1 (MQ=255)
                                                                                        cccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTc                                                               <  2:1066842/150‑1 (MQ=255)
                                                                                         ccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCa                                                              <  1:879253/150‑1 (MQ=255)
                                                                                         ccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCa                                                              <  2:803659/150‑1 (MQ=255)
                                                                                         ccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCa                                                              >  2:936049/1‑150 (MQ=255)
                                                                                            ggCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGtga                                                                              <  2:698880/131‑1 (MQ=255)
                                                                                            ggCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGtga                                                                              >  1:698880/1‑131 (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                                                    >  2:393368/1‑150 (MQ=255)
                                                                                                   gagaAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGc                                                    <  1:353547/150‑1 (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                                                                                                        <  1:808423/77‑1 (MQ=255)
                                                                                                                        gTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGac                                                                                                        >  2:808423/1‑77 (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                                                    <  2:366197/113‑1 (MQ=255)
                                                                                                                                        aaGCTCTGCCCGGCCAGTGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGc                                                    >  1:366197/1‑113 (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)
                                                                                                                                                     |                                                                                                                                                     
TACGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGGCGCAGCCATG  >  NC_002947/5509353‑5509651

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: