Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,575,825 A→T D319E (GAT→GAA)  ydeN ← sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,575,8250AT99.1% 444.0 / NA 110D319E (GAT→GAA) ydeNsulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (63/46);  minor base G (1/0);  total (64/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCAGGATT  >  NZ_CP009273/1575739‑1575911
                                                                                      |                                                                                      
gCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCa                                                                                     >  1:1553736/1‑90 (MQ=255)
 cACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCaa                                                                                    >  1:2376679/1‑90 (MQ=255)
 cACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCaa                                                                                    <  1:2298033/90‑1 (MQ=255)
  aCCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAAc                                                                                   <  2:1919139/90‑1 (MQ=255)
  aCCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAAc                                                                                   <  2:1421976/90‑1 (MQ=255)
  aCCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAAc                                                                                   <  2:355486/90‑1 (MQ=255)
       tATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAAt                                                                              <  1:915259/90‑1 (MQ=255)
        aTCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATa                                                                             >  1:1673405/1‑90 (MQ=255)
             aGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACgg                                                                        >  2:2559171/1‑90 (MQ=255)
             aGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACgg                                                                        >  2:1065426/1‑90 (MQ=255)
              ggTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGa                                                                       <  1:2194649/90‑1 (MQ=255)
              ggTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGa                                                                       <  2:2867966/90‑1 (MQ=255)
                 aaaGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGc                                                                    <  2:1326993/90‑1 (MQ=255)
                      gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgtagt                                                               >  2:1306327/1‑90 (MQ=255)
                      gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgtagt                                                               >  1:1075833/1‑90 (MQ=255)
                      gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgtagt                                                               >  1:2548942/1‑90 (MQ=255)
                      gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgtagt                                                               >  1:783999/1‑90 (MQ=255)
                      gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgtagt                                                               >  1:393745/1‑90 (MQ=255)
                      gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgtagt                                                               >  2:2387361/1‑90 (MQ=255)
                      gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgtagt                                                               <  2:487025/90‑1 (MQ=255)
                        ataatTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGtagtag                                                             >  1:2016559/1‑90 (MQ=255)
                         taatTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTagtagt                                                            >  1:1217367/1‑90 (MQ=255)
                             tgtattgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATc                                                        <  2:31252/90‑1 (MQ=255)
                              gtattgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATct                                                       >  1:2996444/1‑90 (MQ=255)
                                attgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                     >  2:2765700/1‑90 (MQ=255)
                                attgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                     >  1:1666402/1‑90 (MQ=255)
                                attgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                     >  2:1983604/1‑90 (MQ=255)
                                   gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGa                                                                       <  2:3007157/69‑1 (MQ=255)
                                   gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGa                                                                       <  2:839624/69‑1 (MQ=255)
                                   gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGa                                                                       >  1:839624/1‑69 (MQ=255)
                                   gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGa                                                                       >  1:3007157/1‑69 (MQ=255)
                                   gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGt                                                  >  2:1070398/1‑90 (MQ=255)
                                   gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGt                                                  >  1:2677447/1‑90 (MQ=255)
                                     cATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGttt                                                >  2:2204009/1‑90 (MQ=255)
                                     cATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAAC‑‑AATAAACgg                                                                        >  2:1128007/1‑64 (MQ=255)
                                     cATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAAC‑‑AATAAACgg                                                                        <  1:1128007/64‑1 (MQ=255)
                                        aCTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGAc                                             >  2:633719/1‑90 (MQ=255)
                                        aCTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGAc                                             >  1:908574/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCGGTTGTTCGGGAATGGGTTTTACCCCCTGGTCAACTGAATAAAGGGAGGCATAGTAGTTATCTGCTGTTTGACCAc                                          >  1:1728566/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACa                                                                                  >  2:1241435/1‑50 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACa                                                                                  <  1:1241435/50‑1 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  2:1845300/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  2:1243965/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  2:1306087/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  2:685778/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  1:2219654/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  2:1704997/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  1:674360/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  2:2052018/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  1:1699914/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  1:439135/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  1:1451579/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                          >  1:1288960/1‑90 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGCTCGAGAATGCGTTTTACACCCTGTTCAACAGaa                                                                               <  1:1012472/53‑1 (MQ=255)
                                           gTCCGTTTTTCTTCAGTTGCTCGAGAATGCGTTTTACACCCTGTTCAACAGaa                                                                               >  2:1012472/1‑53 (MQ=255)
                                            tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         <  1:1845300/90‑1 (MQ=255)
                                            tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         <  1:1385644/90‑1 (MQ=255)
                                            tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         <  2:1387897/90‑1 (MQ=255)
                                            tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         <  1:2052018/90‑1 (MQ=255)
                                            tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         <  1:669365/90‑1 (MQ=255)
                                            tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         <  2:783999/90‑1 (MQ=255)
                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGaa                                                                      >  1:1604837/1‑60 (MQ=255)
                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGaa                                                                      <  2:1604837/60‑1 (MQ=255)
                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         >  2:114252/1‑89 (MQ=255)
                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         >  1:2069146/1‑89 (MQ=255)
                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         <  2:2069146/89‑1 (MQ=255)
                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         <  1:114252/89‑1 (MQ=255)
                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                        >  1:2780779/1‑90 (MQ=255)
                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                        >  1:564082/1‑90 (MQ=255)
                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                        >  2:1108820/1‑90 (MQ=255)
                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                        >  2:705144/1‑90 (MQ=255)
                                                  tttCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTAtt                                   <  2:1699914/90‑1 (MQ=255)
                                                     cttcAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTaaa                                >  2:1215095/1‑90 (MQ=255)
                                                     cttcAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTaaa                                >  2:841257/1‑90 (MQ=255)
                                                          gttgttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTg                             >  2:1188182/1‑88 (MQ=255)
                                                          gttgttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTg                             <  1:1188182/88‑1 (MQ=255)
                                                            tgttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTc                         >  1:1841769/1‑90 (MQ=255)
                                                            tgttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTc                         <  1:1753365/90‑1 (MQ=255)
                                                             gttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCt                        <  1:2856332/90‑1 (MQ=255)
                                                             gttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCt                        >  1:2287982/1‑90 (MQ=255)
                                                             gttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCt                        >  1:2831097/1‑90 (MQ=255)
                                                              ttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTg                       <  2:1666402/90‑1 (MQ=255)
                                                               tCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                     >  1:545457/1‑59 (MQ=255)
                                                               tCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                     <  2:545457/59‑1 (MQ=255)
                                                                   gaATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATAtt                  <  2:114510/90‑1 (MQ=255)
                                                                       gCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAATTTGCTTCTGATATTGATc              >  1:587130/1‑90 (MQ=255)
                                                                        cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc             <  1:1264890/90‑1 (MQ=255)
                                                                        cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc             <  1:2014250/90‑1 (MQ=255)
                                                                        cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc             <  2:2996444/90‑1 (MQ=255)
                                                                        cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc             <  2:1217367/90‑1 (MQ=255)
                                                                        cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc             <  2:1075833/90‑1 (MQ=255)
                                                                        cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc             <  1:2559171/90‑1 (MQ=255)
                                                                        cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc             <  1:2677577/90‑1 (MQ=255)
                                                                             tACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGc        <  2:868297/90‑1 (MQ=255)
                                                                              acacCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCa       <  2:2387154/90‑1 (MQ=255)
                                                                                 cccTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         >  2:713497/1‑53 (MQ=255)
                                                                                 cccTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                         <  1:713497/53‑1 (MQ=255)
                                                                                 cccTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCAGGa    >  1:2184865/1‑90 (MQ=255)
                                                                                   cTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCAGGAtt  <  2:2047078/90‑1 (MQ=255)
                                                                                   cTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCAGGAtt  <  1:2765700/90‑1 (MQ=255)
                                                                                      |                                                                                      
GCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCAGGATT  >  NZ_CP009273/1575739‑1575911

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