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,3580TC99.5% 647.7 / ‑4.5 190intergenic (‑82/+275)ompA/sulAouter membrane protein 3a/suppressor of lon
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (86/103);  total (87/103)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 194 total aligned reads displayed.

CTTTTTCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTG  >  NC_000913/1019268‑1019447
                                                                                          |                                                                                         
cTTTTTCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGc                                                                                         <  2:2221876/93‑1 (MQ=255)
 tttttCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCt                                                                                        >  1:679330/1‑93 (MQ=255)
 tttttCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCt                                                                                        <  1:6507721/93‑1 (MQ=255)
  ttttCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTTACCGGGGGATCCGCTc                                                                                       <  2:8363208/93‑1 (MQ=255)
  ttttCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTc                                                                                       >  1:7388132/1‑93 (MQ=255)
    ttCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCaa                                                                                     <  2:842809/93‑1 (MQ=255)
     tCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAat                                                                                    <  2:1611029/93‑1 (MQ=255)
      cATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAata                                                                                   <  1:5420159/93‑1 (MQ=255)
      cATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAata                                                                                   <  1:4408756/93‑1 (MQ=255)
       aTTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAatat                                                                                  <  2:1401136/93‑1 (MQ=255)
          ttttGCGACTCGTTATCATCCAAAATACGCCATGAATATCCCCAACGAGATAACACTGTTAAATCCTTCACCGGGGGATCCGCTCAATAttta                                                                               <  2:7203895/93‑3 (MQ=255)
             tGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAActc                                                                            <  2:7721333/93‑1 (MQ=255)
             tGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAActc                                                                            <  2:750165/93‑1 (MQ=255)
             tGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAActc                                                                            <  1:772856/93‑1 (MQ=255)
              gcgcCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCTGGGGATCCGCTCAATATTAActct                                                                           >  2:8055636/1‑93 (MQ=255)
              gcgcCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAActct                                                                           >  2:4474676/1‑93 (MQ=255)
                gcCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTAc                                                                         >  2:3678998/1‑93 (MQ=255)
                gcCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTAc                                                                         >  2:4874666/1‑93 (MQ=255)
                 ccTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTAcc                                                                        >  2:7159797/1‑93 (MQ=255)
                   tCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGa                                                                      >  1:1814574/1‑93 (MQ=255)
                   tCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGa                                                                      >  1:5083258/1‑93 (MQ=255)
                   tCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGa                                                                      >  2:5558579/1‑93 (MQ=255)
                    cGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGat                                                                     <  2:8478535/93‑1 (MQ=255)
                      ttATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGatat                                                                   >  1:8028331/1‑93 (MQ=255)
                      ttATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGatat                                                                   >  1:2935042/1‑93 (MQ=255)
                      ttATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGatat                                                                   >  1:692401/1‑93 (MQ=255)
                       tATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGatatt                                                                  >  2:1331534/1‑92 (MQ=255)
                        atcatcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCt                                                                 >  2:8282373/1‑93 (MQ=255)
                        atcatcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCt                                                                 <  2:2820329/93‑1 (MQ=255)
                        atcatcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACACTACCGATATCt                                                                 <  2:632797/93‑1 (MQ=255)
                          catcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTc                                                               <  2:6794134/93‑1 (MQ=255)
                           atcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCg                                                              <  1:2232251/93‑1 (MQ=255)
                            tcCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCgg                                                             >  2:3190864/1‑93 (MQ=255)
                             ccAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGc                                                            <  1:1085207/93‑1 (MQ=255)
                              cAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCt                                                           <  1:3654848/93‑1 (MQ=255)
                              cAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCt                                                           <  2:5978111/93‑1 (MQ=255)
                               aaaaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCtt                                                          >  2:3437578/1‑93 (MQ=255)
                               aaaaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCtt                                                          <  1:937751/93‑1 (MQ=255)
                                aaaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTa                                                         <  2:3273144/93‑1 (MQ=255)
                                 aaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTAt                                                        <  1:5276822/93‑1 (MQ=255)
                                 aaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTAt                                                        <  1:2175653/93‑1 (MQ=255)
                                 aaTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTAt                                                        <  2:1164420/93‑1 (MQ=255)
                                  aTACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATg                                                       <  2:1088274/93‑1 (MQ=255)
                                   tACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGc                                                      <  2:5129698/93‑1 (MQ=255)
                                   aaCGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGc                                                      <  1:688039/92‑1 (MQ=255)
                                    aCGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGcc                                                     <  1:5426911/93‑1 (MQ=255)
                                    aCGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGcc                                                     <  1:3136399/93‑1 (MQ=255)
                                      gCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGa                                                   >  1:7148439/1‑93 (MQ=255)
                                          tGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAc                                               >  1:4593823/1‑93 (MQ=255)
                                          tGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAc                                               <  2:2502752/93‑1 (MQ=255)
                                           gAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAcc                                              <  2:8153239/93‑1 (MQ=255)
                                           gAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAcc                                              <  1:4174541/93‑1 (MQ=255)
                                            aaTATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAccc                                             >  2:3976877/1‑93 (MQ=255)
                                            aaTATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAccc                                             <  1:6482881/93‑1 (MQ=255)
                                            aaTATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAccc                                             <  1:6931526/93‑1 (MQ=255)
                                            aaTATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAccc                                             >  1:3662371/1‑93 (MQ=255)
                                             atatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAcccc                                            >  2:5453385/1‑93 (MQ=255)
                                             atatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAcccc                                            <  2:2008381/93‑1 (MQ=255)
                                             atatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAcccc                                            <  1:2712471/93‑1 (MQ=255)
                                             atatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCAcccc                                            <  2:2917352/93‑1 (MQ=255)
                                              tatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCt                                           <  2:1227000/93‑1 (MQ=255)
                                              tatCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCt                                           <  1:8252465/93‑1 (MQ=255)
                                                    cAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGAt                                     >  2:1935030/1‑93 (MQ=255)
                                                       cGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGTGATGTa                                  >  1:489791/1‑93 (MQ=255)
                                                       cGAGATAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTa                                  <  2:6068950/93‑1 (MQ=255)
                                                        gagaTAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTaa                                 <  2:3339216/93‑1 (MQ=255)
                                                         agaTAACACGGTTAAATCCTTCACCTGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTaaa                                >  1:7020069/1‑93 (MQ=255)
                                                         agaTAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTaaa                                <  2:8716230/93‑1 (MQ=255)
                                                            tAACACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTc                             >  1:1473068/1‑93 (MQ=255)
                                                             aaCACGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCt                            <  1:3519411/93‑1 (MQ=255)
                                                              acacGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTa                           <  1:2695207/93‑1 (MQ=255)
                                                               cacGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTAc                          <  1:386215/93‑1 (MQ=255)
                                                               cacGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTAc                          <  2:3695459/93‑1 (MQ=255)
                                                                acGGTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACa                         >  1:157005/1‑93 (MQ=255)
                                                                  ggTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAAc                       <  1:7914601/93‑1 (MQ=255)
                                                                  ggTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAAc                       >  2:6755817/1‑93 (MQ=255)
                                                                  ggTTAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAAc                       >  2:567517/1‑93 (MQ=255)
                                                                     tAAATCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTa                    <  2:6286755/93‑1 (MQ=255)
                                                                       aaTCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGt                  <  1:1161459/93‑1 (MQ=255)
                                                                       aaTCCTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGt                  <  1:8702565/93‑1 (MQ=255)
                                                                          ccTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGa               >  1:3941419/1‑93 (MQ=255)
                                                                          ccTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGa               >  1:4780341/1‑93 (MQ=255)
                                                                          ccTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGa               >  1:5700513/1‑93 (MQ=255)
                                                                          ccTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGACGTAAAGTCTACAACGTAGTTGa               >  2:4367070/1‑93 (MQ=255)
                                                                           cTTCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGaa              >  2:2462680/1‑93 (MQ=255)
                                                                            ttCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGaaa             >  2:790242/1‑93 (MQ=255)
                                                                            ttCACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGaaa             >  1:749335/1‑93 (MQ=255)
                                                                              cACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAAc           >  2:884607/1‑93 (MQ=255)
                                                                              cACCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAAc           <  2:6452486/93‑1 (MQ=255)
                                                                               aCCGGGGGATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACt          >  1:5857692/1‑93 (MQ=255)
                                                                                 cGGGGTATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCCCCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTa        <  2:7926189/93‑1 (MQ=255)
                                                                                  gggggATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTAc       >  2:7461383/1‑93 (MQ=255)
                                                                                   ggggATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACa      >  1:2838896/1‑93 (MQ=255)
                                                                                   ggggATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCCAGCACCCCTGGCGATTTAAAGTCTACAACGTAGTTGAAAACTTACa      <  2:4313149/93‑1 (MQ=255)
                                                                                    gggATCCGCTCATTATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACaa     <  2:5053768/93‑1 (MQ=255)
                                                                                    gggATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACaa     >  1:592195/1‑93 (MQ=255)
                                                                                     ggATCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAg    <  1:2602882/93‑1 (MQ=255)
                                                                                       aTCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAgtg  <  1:8364278/93‑1 (MQ=255)
                                                                                       aTCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAgtg  <  1:847846/93‑1 (MQ=255)
                                                                                       aTCCGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACACCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAgtg  >  1:6080241/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTTTTTCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTG  >  NC_000913/1019268‑1019447

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