Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,392,826 T→C Q18Q (CAA→CAG uspE ← universal stress protein UspE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,392,8260TC98.9% 346.7 / ‑4.0 93Q18Q (CAA→CAGuspEuniversal stress protein UspE
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (34/58);  total (34/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

TTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGC  >  NZ_CP009273/1392739‑1392909
                                                                                       |                                                                                   
ttCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATgctggctg                                                                                   >  2:2315950/1‑90 (MQ=255)
  cGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGctggctggc                                                                                 <  2:159949/90‑2 (MQ=255)
  cGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCtggctggt                                                                                 <  2:325843/90‑1 (MQ=255)
    tATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGgtcg                                                                               >  1:223516/1‑90 (MQ=255)
       gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTcc                                                                            <  1:1087731/90‑1 (MQ=255)
       gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                            <  1:1134747/90‑1 (MQ=255)
       gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                            <  2:2620353/90‑1 (MQ=255)
       gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                            >  2:1573530/1‑90 (MQ=255)
       gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                            <  1:708176/90‑1 (MQ=255)
       gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                            <  2:2888989/90‑1 (MQ=255)
        agaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCt                                                                           <  2:1818672/90‑1 (MQ=255)
         gaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTg                                                                          >  1:2835614/1‑90 (MQ=255)
             tCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTa                                                                      >  2:2757235/1‑90 (MQ=255)
              cATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAg                                                                     >  1:1739412/1‑90 (MQ=255)
                tAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGa                                                                   <  2:2986701/90‑1 (MQ=255)
                 aGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGAt                                                                  >  2:278737/1‑90 (MQ=255)
                  gATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATc                                                                 <  1:2710175/90‑1 (MQ=255)
                   aTCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCg                                                                >  2:1927552/1‑90 (MQ=255)
                    tCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGa                                                               >  2:746412/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  2:1414966/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  1:2569724/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  2:770733/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             <  2:1465488/90‑1 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  2:2491631/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             <  1:1463430/90‑1 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  1:1679408/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  1:2183848/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  2:740363/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  2:399288/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  1:1830237/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  2:637330/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  2:1843204/1‑90 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGACGGCCTGGTTAGGATGGATa                                                             >  1:509327/1‑90 (MQ=255)
                           aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                                        <  2:1938451/90‑1 (MQ=255)
                           aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                                        <  1:399288/90‑1 (MQ=255)
                           aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                                        <  2:1679408/90‑1 (MQ=255)
                           aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                                        <  1:2491631/90‑1 (MQ=255)
                           aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                                        <  1:2430473/90‑1 (MQ=255)
                           aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                                        <  1:1843204/90‑1 (MQ=255)
                           aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                                        <  2:677896/90‑1 (MQ=255)
                               gCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGc                                                    <  1:770733/90‑1 (MQ=255)
                                cTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCa                                                   <  1:2315950/90‑1 (MQ=255)
                                cTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCa                                                   <  1:1414966/90‑1 (MQ=255)
                                cTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCa                                                   <  2:1230121/90‑1 (MQ=255)
                                cTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCa                                                   <  2:1830237/90‑1 (MQ=255)
                                 tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCAt                                                  <  1:805971/90‑1 (MQ=255)
                                 tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCAt                                                  <  1:995055/90‑1 (MQ=255)
                                   tAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                                                <  2:767344/90‑1 (MQ=255)
                                   tAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                                                <  1:185181/90‑1 (MQ=255)
                                     aTTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTc                                              <  1:2542313/90‑1 (MQ=255)
                                     aTTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTc                                              <  1:2757235/90‑1 (MQ=255)
                                        ttGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGa                                           >  1:1486354/1‑90 (MQ=255)
                                           ccaccaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATAc                                        <  1:1772421/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  2:767005/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  2:1739412/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  2:2835614/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  1:1113105/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  2:1695504/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  1:1578952/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  1:1796480/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  1:746412/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  2:1486354/90‑1 (MQ=255)
                                                 aTCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGcc                                  <  1:278737/90‑1 (MQ=255)
                                                    cGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATa                               <  2:550857/90‑1 (MQ=255)
                                                     gTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATAc                              >  1:1609025/1‑90 (MQ=255)
                                                        gATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGc                                                    <  2:494588/65‑1 (MQ=255)
                                                        gATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGc                                                    >  1:494588/1‑65 (MQ=255)
                                                        gATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAgg                           <  1:1850529/90‑1 (MQ=255)
                                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc                        <  2:1145954/90‑1 (MQ=255)
                                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc                        <  2:757178/90‑1 (MQ=255)
                                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc                        <  2:609503/90‑1 (MQ=255)
                                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc                        <  1:2192374/90‑1 (MQ=255)
                                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc                        <  2:2199693/90‑1 (MQ=255)
                                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc                        <  1:968587/90‑1 (MQ=255)
                                                                  aaCAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGc                                   >  2:2348645/1‑72 (MQ=255)
                                                                  aaCAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGc                                   <  1:2348645/72‑1 (MQ=255)
                                                                  aaCAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc                        <  1:2057068/83‑1 (MQ=255)
                                                                  aaCAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc                        >  2:2057068/1‑83 (MQ=255)
                                                                  aaCAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTac                 >  1:1425605/1‑90 (MQ=255)
                                                                   aCAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  1:2041284/68‑1 (MQ=255)
                                                                   aCAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      >  2:2041284/1‑68 (MQ=255)
                                                                     aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCa                                                   <  1:1402440/53‑1 (MQ=255)
                                                                     aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCa                                                   >  2:1402440/1‑53 (MQ=255)
                                                                     aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                                                <  2:813145/56‑1 (MQ=255)
                                                                     aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                                                >  1:813145/1‑56 (MQ=255)
                                                                     aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      >  1:1270120/1‑66 (MQ=255)
                                                                     aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  2:1270120/66‑1 (MQ=255)
                                                                     aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTacaac              >  2:1730395/1‑90 (MQ=255)
                                                                      gCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                                                >  2:711524/1‑55 (MQ=255)
                                                                      gCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                                                <  1:711524/55‑1 (MQ=255)
                                                                      gCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTacaaca             >  1:1525126/1‑90 (MQ=255)
                                                                      gCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTacaaca             >  2:1986698/1‑90 (MQ=255)
                                                                            cgcAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTc       <  2:2351672/90‑1 (MQ=255)
                                                                                aTGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACg   <  2:910434/90‑1 (MQ=255)
                                                                                 tGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGc  <  1:1927552/90‑1 (MQ=255)
                                                                                 tGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGc  <  2:866716/90‑1 (MQ=255)
                                                                                       |                                                                                   
TTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGC  >  NZ_CP009273/1392739‑1392909

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