Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,509,493 A→T 10.9% A350A (GCA→GCT 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,4930AT10.9% 238.5 / 17.8 83A350A (GCA→GCTurtAurea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (35/38);  new base T (6/3);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GCCTACGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGGC  >  NC_002947/5509350‑5509642
                                                                                                                                               |                                                                                                                                                     
gCCTACGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGcc                                                                                                                                                 >  2:640416/1‑150 (MQ=255)
  cTACGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCgg                                                                                                                                               >  2:577227/1‑150 (MQ=255)
  cTACGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGAAAAGGTCACCGACCCCCTGGCCgg                                                                                                                                               >  2:1148881/1‑150 (MQ=255)
   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                                                                                           >  1:1174024/1‑150 (MQ=255)
                                                      cccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaac                                                                                           <  2:1324834/150‑1 (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                                                                               >  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: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                                                     <  2:803659/150‑1 (MQ=255)
                                                                                            ccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCa                                                     >  2:936049/1‑150 (MQ=255)
                                                                                            ccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCa                                                     <  1:879253/150‑1 (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:353547/150‑1 (MQ=255)
                                                                                                      gagaAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGc                                           <  1:942724/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                                           <  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  <  2:560394/150‑1 (MQ=255)
                                                                                                                                               aCTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGgc  >  1:92294/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GCCTACGCCAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCAGCCGCTGCGGGC  >  NC_002947/5509350‑5509642

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: