Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,019,358 T→C intergenic (‑82/+275) ompA ← / ← sulA outer membrane protein 3a/suppressor of lon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,019,3580TC89.7% 615.7 / 53.4 232intergenic (‑82/+275)ompA/sulAouter membrane protein 3a/suppressor of lon
Reads supporting (aligned to +/- strand):  ref base T (18/6);  new base C (104/104);  total (122/110)
Fisher's exact test for biased strand distribution p-value = 2.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 233 total aligned reads displayed.

GTCTTTTTCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTG  >  NC_000913/1019266‑1019449
                                                                                            |                                                                                           
   tttttCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCt                                                                                          >  1:6587852/1‑93 (MQ=255)
   tttttCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCt                                                                                          >  2:6565710/1‑93 (MQ=255)
      ttCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCaa                                                                                       >  2:3891105/1‑93 (MQ=255)
      ttCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCaa                                                                                       >  1:2986095/1‑93 (MQ=255)
       tCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAat                                                                                      <  2:2351413/93‑1 (MQ=255)
         aTTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAatat                                                                                    >  2:4942630/1‑93 (MQ=255)
         aTTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAatat                                                                                    >  2:5441224/1‑93 (MQ=255)
           tttttGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTa                                                                                  <  2:7068737/93‑1 (MQ=255)
            ttttGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTaa                                                                                 <  1:7451848/93‑1 (MQ=255)
             tttGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAAc                                                                                <  2:7494326/93‑1 (MQ=255)
               tGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAActc                                                                              <  2:2342298/93‑1 (MQ=255)
                 cgcCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTa                                                                            <  1:3114041/93‑1 (MQ=255)
                  gcCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTAc                                                                           >  1:2015776/1‑93 (MQ=255)
                  gcCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTAc                                                                           <  2:1385524/93‑1 (MQ=255)
                    cTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCg                                                                         >  1:7172640/1‑93 (MQ=255)
                    cTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCg                                                                         >  1:816949/1‑93 (MQ=255)
                    cTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCg                                                                         >  2:761657/1‑93 (MQ=255)
                     tCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGa                                                                        >  1:2247/1‑93 (MQ=255)
                     tCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGa                                                                        <  2:3350773/93‑1 (MQ=255)
                     tCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGa                                                                        >  1:7387837/1‑93 (MQ=255)
                      cGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGat                                                                       <  2:1856885/93‑1 (MQ=255)
                       gTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGata                                                                      >  1:3892995/1‑93 (MQ=255)
                        ttATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGatat                                                                     >  1:6820796/1‑93 (MQ=255)
                        ttATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGatat                                                                     >  1:5087525/1‑93 (MQ=255)
                         tATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATc                                                                    >  1:3711842/1‑93 (MQ=255)
                         tATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATc                                                                    >  2:3135774/1‑93 (MQ=255)
                            catcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTc                                                                 <  2:178721/93‑1 (MQ=255)
                             atcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCg                                                                <  1:1788709/93‑1 (MQ=255)
                              tcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCgg                                                               >  2:6256853/1‑93 (MQ=255)
                              tcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCgg                                                               >  1:5999938/1‑93 (MQ=255)
                              tcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCgg                                                               >  2:5523758/1‑93 (MQ=255)
                               ccAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGc                                                              <  1:920429/93‑1 (MQ=255)
                               ccAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGc                                                              >  1:7481810/1‑93 (MQ=255)
                                 aaaaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCtt                                                            <  1:1912770/93‑1 (MQ=255)
                                  aaaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCTATATCTTCGGCTTa                                                           <  2:6444484/93‑1 (MQ=255)
                                  aaaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTa                                                           <  1:4921932/93‑1 (MQ=255)
                                   aaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTAt                                                          >  2:6006238/1‑93 (MQ=255)
                                   aaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTAt                                                          >  2:2412982/1‑93 (MQ=255)
                                        gCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGa                                                     <  2:2915173/93‑1 (MQ=255)
                                        gCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGa                                                     >  1:2887151/1‑93 (MQ=255)
                                         ccATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAg                                                    <  2:1604793/93‑1 (MQ=255)
                                          cATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGc                                                   <  1:4044644/93‑1 (MQ=255)
                                            tGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAc                                                 >  2:4201192/1‑93 (MQ=255)
                                            tGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAc                                                 >  1:3843754/1‑93 (MQ=255)
                                            tGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAc                                                 >  2:1986052/1‑93 (MQ=255)
                                              aaTATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAccc                                               <  1:6491078/93‑1 (MQ=255)
                                               atatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAcccc                                              >  2:6419330/1‑93 (MQ=255)
                                               atatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAcccc                                              <  1:5292651/93‑1 (MQ=255)
                                               atatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAcccc                                              >  1:1374134/1‑93 (MQ=255)
                                                tatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCt                                             >  2:5627236/1‑93 (MQ=255)
                                                tatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCt                                             >  1:7317489/1‑93 (MQ=255)
                                                  tctcCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTgg                                           >  2:5026627/1‑93 (MQ=255)
                                                   ctcCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGAACCCCTGGc                                          <  2:1375717/93‑1 (MQ=255)
                                                   ctcCAACGAGATAACACGGGTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGc                                          <  1:547605/93‑1 (MQ=255)
                                                    tcCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCg                                         >  2:5257422/1‑93 (MQ=255)
                                                     ccAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGa                                        <  2:6768956/93‑1 (MQ=255)
                                                     ccAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGa                                        <  1:3687625/93‑1 (MQ=255)
                                                     ccAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGa                                        <  1:3297417/93‑1 (MQ=255)
                                                      cAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGAt                                       <  2:2520942/93‑1 (MQ=255)
                                                      cAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGAt                                       >  2:1052410/1‑93 (MQ=255)
                                                       aaCGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATg                                      >  2:274862/1‑93 (MQ=255)
                                                       aaCGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATg                                      >  2:6335951/1‑93 (MQ=255)
                                                       aaCGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATg                                      <  1:236622/93‑1 (MQ=255)
                                                       aaCGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATg                                      >  2:6995227/1‑93 (MQ=255)
                                                        aCGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGt                                     <  1:5689545/93‑1 (MQ=255)
                                                        aCGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGt                                     >  2:5836799/1‑93 (MQ=255)
                                                           agaTAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTaaa                                  >  2:4688714/1‑93 (MQ=255)
                                                            gaTAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAg                                 <  2:2807873/93‑1 (MQ=255)
                                                             aTAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGt                                <  2:3728516/93‑1 (MQ=255)
                                                                 cacGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTAc                            >  1:5327300/1‑93 (MQ=255)
                                                                  acGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACa                           >  2:5814562/1‑93 (MQ=255)
                                                                   cGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACaa                          <  1:1028794/93‑1 (MQ=255)
                                                                   cGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACaa                          <  2:1248421/93‑1 (MQ=255)
                                                                   cGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACaa                          >  2:7178343/1‑93 (MQ=255)
                                                                   cGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACaa                          >  2:4412427/1‑93 (MQ=255)
                                                                    ggTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAAc                         >  1:5155917/1‑93 (MQ=255)
                                                                    ggTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAAc                         >  2:1103744/1‑93 (MQ=255)
                                                                     gTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACg                        >  1:4304225/1‑93 (MQ=255)
                                                                       tAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTa                      <  2:6593661/93‑1 (MQ=255)
                                                                       tAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTa                      <  1:4172618/93‑1 (MQ=255)
                                                                       tAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTa                      <  1:476094/93‑1 (MQ=255)
                                                                         aaTCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGt                    <  2:897839/93‑1 (MQ=255)
                                                                         aaTCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGt                    <  2:4031247/93‑1 (MQ=255)
                                                                          aTCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGtt                   >  1:6536924/1‑93 (MQ=255)
                                                                          aTCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGtt                   <  2:2830213/93‑1 (MQ=255)
                                                                             cTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGaa                >  1:2653677/1‑93 (MQ=255)
                                                                               tCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGaaaa              <  2:6683643/93‑1 (MQ=255)
                                                                                cACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAAc             <  2:3463330/93‑1 (MQ=255)
                                                                                cACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAAc             <  1:161920/93‑1 (MQ=255)
                                                                                cACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAAc             <  2:2655647/93‑1 (MQ=255)
                                                                                cACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAAc             <  1:6911540/93‑1 (MQ=255)
                                                                                   cGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTa          <  1:443071/93‑1 (MQ=255)
                                                                                   cGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTa          <  1:2749/93‑1 (MQ=255)
                                                                                     ggggATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACa        >  1:4338001/1‑93 (MQ=255)
                                                                                     ggggATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACa        <  1:4634661/93‑1 (MQ=255)
                                                                                     ggggATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACa        <  1:127110/93‑1 (MQ=255)
                                                                                      gggATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACaa       >  2:5917259/1‑93 (MQ=255)
                                                                                       ggATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAg      <  1:2574686/93‑1 (MQ=255)
                                                                                        gATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAgt     <  1:5825813/93‑1 (MQ=255)
                                                                                           cTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGtgtg  >  1:17386/1‑93 (MQ=255)
                                                                                            |                                                                                           
GTCTTTTTCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTG  >  NC_000913/1019266‑1019449

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