Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 232,840 T→A V96E (GTG→GAG)  dnaQ → DNA polymerase III subunit epsilon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273232,8400TA97.5% 314.9 / ‑0.9 81V96E (GTG→GAG) dnaQDNA polymerase III subunit epsilon
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (38/41);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

GTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTT  >  NZ_CP009273/232753‑232923
                                                                                         |                                                                                   
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  1:2551905/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  1:2412000/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  1:758668/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  1:2371821/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  1:1976414/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  2:894727/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   >  1:1669185/1‑90 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  2:2390055/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   >  2:2360834/1‑90 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  1:2864751/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  1:1209741/90‑1 (MQ=255)
  aTTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATc                                                                                 <  1:2585533/90‑1 (MQ=255)
  aTTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATc                                                                                 <  2:2057443/90‑1 (MQ=255)
     gCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCAt                                                                              >  2:1294550/1‑90 (MQ=255)
     gCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCAt                                                                              >  2:2100930/1‑90 (MQ=255)
      ccGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATa                                                                             <  1:924836/90‑1 (MQ=255)
        gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   >  1:351927/1‑82 (MQ=255)
        gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  2:351927/82‑1 (MQ=255)
        gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAAc                                                                           >  2:1271859/1‑90 (MQ=255)
        gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAAc                                                                           >  2:1963420/1‑90 (MQ=255)
        gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAAc                                                                           >  2:1181421/1‑90 (MQ=255)
          tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   >  2:973374/1‑80 (MQ=255)
          tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  1:973374/80‑1 (MQ=255)
          tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGc                                                                         >  2:2883919/1‑90 (MQ=255)
          tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGc                                                                         >  1:2949698/1‑90 (MQ=255)
          tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGc                                                                         >  2:2450340/1‑90 (MQ=255)
            aaTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGTGATCCATAACGCAg                                                                       <  2:2303088/90‑1 (MQ=255)
            aaTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAg                                                                       <  2:1605438/90‑1 (MQ=255)
                 ttGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTc                                                                  <  1:1895721/90‑1 (MQ=255)
                       gATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATc                                                            <  1:511829/90‑1 (MQ=255)
                       gATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATc                                                            <  1:1422234/90‑1 (MQ=255)
                       gATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATc                                                            <  1:2678093/90‑1 (MQ=255)
                       gATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATc                                                            <  1:1181421/90‑1 (MQ=255)
                       gATAAGCCCACCTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATc                                                            >  2:2957246/1‑90 (MQ=255)
                               cACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTAt                                                    <  1:675927/90‑1 (MQ=255)
                               cACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTAt                                                    <  1:1072804/90‑1 (MQ=255)
                                 cGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATgg                                                  <  1:162396/90‑1 (MQ=255)
                                    ttGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACt                                               >  2:2770336/1‑90 (MQ=255)
                                     tGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   >  1:394053/1‑53 (MQ=255)
                                     tGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  2:394053/53‑1 (MQ=255)
                                      gCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTAc                                             >  2:1882285/1‑90 (MQ=255)
                                      gCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTAc                                             >  2:1065567/1‑90 (MQ=255)
                                               gCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCg                                    >  2:2440696/1‑90 (MQ=255)
                                                ccGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGt                                   <  2:2949698/90‑1 (MQ=255)
                                                   aTGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   <  1:2442648/39‑1 (MQ=255)
                                                   aTGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGgaga                                                                                   >  2:2442648/1‑39 (MQ=255)
                                                      aGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTa                             <  1:2409100/90‑1 (MQ=255)
                                                       gTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTaa                            >  1:1306142/1‑90 (MQ=255)
                                                         tCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgc                          <  1:2692744/90‑1 (MQ=255)
                                                         tCATGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgc                          >  1:2414246/1‑90 (MQ=255)
                                                           aTGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgcgc                        <  2:2290055/90‑1 (MQ=255)
                                                           aTGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgcgc                        <  1:1294550/90‑1 (MQ=255)
                                                           aTGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgcgc                        <  1:2844508/90‑1 (MQ=255)
                                                           aTGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgcgc                        <  1:2450340/90‑1 (MQ=255)
                                                           aTGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgcgc                        <  2:2759925/90‑1 (MQ=255)
                                                            tGGACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGcgcg                       >  1:1318015/1‑90 (MQ=255)
                                                             ggACTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGa                      <  2:778629/90‑1 (MQ=255)
                                                               aCTATATTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGata                    <  1:1176884/90‑1 (MQ=255)
                                                                 tataTTCGCGGCGCGGAGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGa                      <  2:1331263/88‑1 (MQ=255)
                                                                 tataTTCGCGGCGCGGAGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGa                      >  1:1331263/1‑88 (MQ=255)
                                                                 tataTTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTc                                                                  <  2:1297583/42‑1 (MQ=255)
                                                                 tataTTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTc                                                                  >  1:1297583/1‑42 (MQ=255)
                                                                 tataTTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCgg                                                          >  2:961289/1‑50 (MQ=255)
                                                                 tataTTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCgg                                                          <  1:961289/50‑1 (MQ=255)
                                                                 tataTTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATAtt                  >  2:387084/1‑90 (MQ=255)
                                                                 tataTTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATAtt                  >  2:212574/1‑90 (MQ=255)
                                                                 tataTTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATAtt                  >  1:2120159/1‑90 (MQ=255)
                                                                 tataTTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATAtt                  >  2:1710450/1‑90 (MQ=255)
                                                                 tataTTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATAtt                  >  1:832975/1‑90 (MQ=255)
                                                                 tataTTCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCAGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATAtt                  >  1:2400574/1‑90 (MQ=255)
                                                                     ttCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGa              <  1:2100930/90‑1 (MQ=255)
                                                                     ttCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGa              >  2:2058345/1‑90 (MQ=255)
                                                                     ttCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGa              <  2:133998/90‑1 (MQ=255)
                                                                     ttCGCGGCGCG‑‑GAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGa              <  2:2167874/90‑1 (MQ=255)
                                                                       cgcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAg            <  1:1072862/90‑1 (MQ=255)
                                                                       cgcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAg            <  1:2883919/90‑1 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATc                                                            >  2:995561/1‑41 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATc                                                            >  2:331548/1‑41 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATc                                                            <  1:995561/41‑1 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATc                                                            <  1:331548/41‑1 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGa                      <  1:133039/79‑1 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGa                      >  2:133039/1‑79 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGa              <  1:330393/87‑1 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGa              >  2:330393/1‑87 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa           >  1:566563/1‑90 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa           >  1:306903/1‑90 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa           >  2:411563/1‑90 (MQ=255)
                                                                        gcggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa           >  2:1777681/1‑90 (MQ=255)
                                                                         cggcgcg‑‑gAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGAc          >  1:2963506/1‑90 (MQ=255)
                                                                          ggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAg                           >  2:1417821/1‑72 (MQ=255)
                                                                          ggcgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAg                           <  1:1417821/72‑1 (MQ=255)
                                                                            cgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGata                    >  2:2291838/1‑77 (MQ=255)
                                                                            cgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGata                    <  1:2291838/77‑1 (MQ=255)
                                                                            cgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCaa       >  2:2838076/1‑90 (MQ=255)
                                                                            cgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCaa       >  2:1202212/1‑90 (MQ=255)
                                                                            cgcg‑‑gAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCaa       <  2:2120159/90‑1 (MQ=255)
                                                                                   aGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACtt  <  1:354564/90‑1 (MQ=255)
                                                                                   aGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACtt  <  2:2659519/90‑1 (MQ=255)
                                                                                   aGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACtt  <  1:2440696/90‑1 (MQ=255)
                                                                                   aGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACtt  <  1:1682063/90‑1 (MQ=255)
                                                                                         |                                                                                   
GTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCG‑‑GAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAATACTT  >  NZ_CP009273/232753‑232923

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: