Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,805,970 G→A 6.3% S38S (AGC→AGT ydiY ← acid‑inducible putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,805,9700GA 6.3% 198.7 / 2.6 63S38S (AGC→AGTydiYacid‑inducible putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (34/25);  new base A (2/2);  total (36/27)
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.06e-01

TATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCA  >  NC_000913/1805825‑1806117
                                                                                                                                                 |                                                                                                                                                   
tATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTg                                                                                                                                                  >  4:189757/1‑149 (MQ=255)
     cTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAg                                                                                                                                              >  5:337202/1‑148 (MQ=255)
        gCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGat                                                                                                                                           <  1:290500/148‑1 (MQ=255)
              ccccACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCAt                                                                                                                                    <  6:19650/149‑1 (MQ=255)
                 cacaGCGCCCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGAGCATCCATGa                                                                                                                                  >  4:400661/1‑148 (MQ=255)
                    aGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAg                                                                                                                               >  4:322027/1‑148 (MQ=255)
                       gACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACa                                                                                                                                >  2:386250/1‑144 (MQ=255)
                       gACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACa                                                                                                                                <  1:386250/144‑1 (MQ=255)
                        aCCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTaaa                                                                                                                           >  1:117172/1‑148 (MQ=255)
                         ccAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTg                                                                                                                                                  >  5:87551/1‑124 (MQ=255)
                         ccAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTg                                                                                                                                                  <  6:87551/124‑1 (MQ=255)
                                       gtCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCa                                                                                                           <  2:122016/149‑1 (MQ=255)
                                        tCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAg                                                                                                          >  3:261435/1‑149 (MQ=255)
                                          cATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGc                                                                                                         <  3:189757/148‑1 (MQ=255)
                                            tACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCACTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTg                                                                                                       >  8:171002/1‑148 (MQ=255)
                                              ccAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCa                                                                                                     <  6:337202/148‑1 (MQ=255)
                                               cAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCAtt                                                                                                   >  2:219271/1‑149 (MQ=255)
                                               cAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCAtt                                                                                                   >  4:355821/1‑149 (MQ=255)
                                               cAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGAGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGAAGCTGCAtt                                                                                                   >  5:331571/1‑149 (MQ=255)
                                                  gTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAg                                                                                                >  5:21146/1‑149 (MQ=255)
                                                              tCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACa                                                                                     <  4:208580/148‑1 (MQ=255)
                                                                  ccGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGccc                                                                                <  1:231230/149‑1 (MQ=255)
                                                                  ccGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCACTGGCAGGATCATCCATGACAGTAAAAACGGa                                                                                                                     >  5:209462/1‑112 (MQ=255)
                                                                  ccGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCACTGGCAGGATCATCCATGACAGTAAAAACGGa                                                                                                                     <  6:209462/112‑1 (MQ=255)
                                                                   cGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCa                                                                                                 <  2:383223/131‑1 (MQ=255)
                                                                   cGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCa                                                                                                 >  1:383223/1‑131 (MQ=255)
                                                                         tGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGcc                                                                          >  1:236393/1‑148 (MQ=255)
                                                                           aggagCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAAc                                                                       <  3:355821/149‑1 (MQ=255)
                                                                             gagCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCTCCCCGCCAACa                                                                      >  3:210383/1‑148 (MQ=255)
                                                                             gagCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACa                                                                      >  2:357141/1‑148 (MQ=255)
                                                                              agCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACAt                                                                     >  8:420890/1‑148 (MQ=255)
                                                                                  tttCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGc                                                                                                            >  7:240975/1‑105 (MQ=255)
                                                                                  tttCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGc                                                                                                            <  8:240975/105‑1 (MQ=255)
                                                                                   ttCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAg                                                                                                >  2:253117/1‑116 (MQ=255)
                                                                                   ttCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAg                                                                                                >  2:259529/1‑116 (MQ=255)
                                                                                   ttCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAg                                                                                                <  1:259529/116‑1 (MQ=255)
                                                                                   ttCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAg                                                                                                <  1:253117/116‑1 (MQ=255)
                                                                                     cGTTTTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCCGGATCATCCATGACAGTAAAACCGGAATCATCGGCCGCTGCATTCAGTGACGCAAACCTGCCCCCCGCCAACTTAACGATgg                                                             >  6:405472/1‑149 (MQ=255)
                                                                                       tATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCa                                                           <  2:236393/149‑1 (MQ=255)
                                                                                         ttGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAg                                                          >  7:254396/1‑148 (MQ=255)
                                                                                            ccTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGc                                                                                  <  8:394375/121‑1 (MQ=255)
                                                                                             cTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAg                                                                                                          <  2:192512/96‑1 (MQ=255)
                                                                                             cTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACt                                                     <  6:21146/149‑1 (MQ=255)
                                                                                               gATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGt                                                   >  5:303166/1‑149 (MQ=255)
                                                                                               gATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGt                                                   >  7:56871/1‑149 (MQ=255)
                                                                                                     gCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCa                                              >  7:146042/1‑148 (MQ=255)
                                                                                                       aaGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAAGCATCGGCAGCTGCATGCAGTGCCGCAACCATGCCCCCCGCCAACAGAACGAGGGCAGGAACTGTCTTCaa                                             >  7:164931/1‑147 (MQ=255)
                                                                                                       aaGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGAGCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCAGTCAGTGCCGCAAACATGCCCCCCGCCAACATAACGCTGGCAGGCACTGTCTTCaa                                             >  7:158927/1‑147 (MQ=255)
                                                                                                             tCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCACTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCa                                     <  7:171002/149‑1 (MQ=255)
                                                                                                                   gTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTAt                                >  7:7340/1‑148 (MQ=255)
                                                                                                                      cAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGCCAGTAAAAACGGAAACATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCaa                             <  8:174401/148‑1 (MQ=255)
                                                                                                                      cAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCaa                             <  8:254396/148‑1 (MQ=255)
                                                                                                                        gATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAagag                          >  8:81924/1‑149 (MQ=255)
                                                                                                                          ttACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATTCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTc                        >  5:93364/1‑149 (MQ=255)
                                                                                                                              cTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAAGTGTCTTCAAAAGCTTCATTTAGCAAGAGTCCGTa                    >  2:298314/1‑149 (MQ=255)
                                                                                                                                  gAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCACCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTacaa                 >  1:96448/1‑148 (MQ=255)
                                                                                                                                   aaaGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGGATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTacaac                >  5:90754/1‑148 (MQ=255)
                                                                                                                                   aaaGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTacaac                >  5:79687/1‑148 (MQ=255)
                                                                                                                                       gTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACaaaaa           <  4:261435/149‑1 (MQ=255)
                                                                                                                                            ttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAaga       <  6:331571/148‑1 (MQ=255)
                                                                                                                                              aGCGCTGGCAGGATCATCCATGACAGTAAAAACGGGATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAgaga     <  6:90754/148‑1 (MQ=255)
                                                                                                                                              aGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAgaga     <  1:357141/148‑1 (MQ=255)
                                                                                                                                              aGCGCGGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCATTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAgaga     <  2:404247/148‑1 (MQ=255)
                                                                                                                                               gcgcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGAGGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGAcc   <  2:117172/149‑1 (MQ=255)
                                                                                                                                                cgcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCa  <  8:56871/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCA  >  NC_000913/1805825‑1806117

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