Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,616,560 G→C 6.4% S156T (AGC→ACC) ‡ bepA → beta‑barrel assembly‑enhancing protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,616,5600GC 6.4% 240.4 / 4.8 73S156T (AGC→ACC) ‡bepAbeta‑barrel assembly‑enhancing protease
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base C (2/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

TGCTTTCTTTGGCGGCAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACC  >  NC_000913/2616414‑2616709
                                                                                                                                                  |                                                                                                                                                     
tGCTTTCTTTGGCGGCAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCagcgcagcgc                                                                                                                                                    <  1:985298/150‑1 (MQ=255)
tGCTTTCTTTCGCGGAAACGTGGTGATGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCAACCTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAgcagcgcagccc                                                                                                                                                    >  2:371601/1‑148 (MQ=255)
  ctttctttGGCGGCAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAgcgcgc                                                                                                                                                  >  2:1160061/1‑150 (MQ=255)
      ctttGGCGGCAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCt                                                                                                                                              >  2:79434/1‑150 (MQ=255)
        ttGGCGGCAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGa                                                                                                                                            <  1:598263/150‑1 (MQ=255)
           gcggcAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGGTCAGCGGCGCGCTGCGCTGCTGATCt                                                                                                                                         <  1:1435651/150‑1 (MQ=255)
           gcggcAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCt                                                                                                                                         >  2:983499/1‑150 (MQ=255)
            cggcAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTg                                                                                                                                        >  2:14225/1‑150 (MQ=255)
             ggcAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTgg                                                                                                                                       >  2:1031606/1‑150 (MQ=255)
              gcAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTggg                                                                                                                                      <  1:80234/150‑1 (MQ=255)
               cAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCGTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGt                                                                                                                                     >  2:1158495/1‑150 (MQ=255)
                aaCGTGGTGCGGCACGCTGCCCTGGTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTc                                                                                                                                    <  1:209586/150‑1 (MQ=255)
                 aCGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGGTCAGCGGCGCGCTGCGCTGCTGATCttccatc                                                                                                                                   <  1:1299266/150‑7 (MQ=255)
                      gtgCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGcgcg                                                                                                                              <  1:970386/150‑1 (MQ=255)
                       tgctgcACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGcgcc                                                                                                                                                 <  1:1162637/130‑1 (MQ=255)
                       tgctgcACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGcgcc                                                                                                                                                 >  2:1162637/1‑130 (MQ=255)
                       tgctgcACTCAGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGAACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATAAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGt                                                                                                                             >  2:409241/1‑150 (MQ=255)
                        gctgcACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGtt                                                                                                                            >  1:238430/1‑150 (MQ=255)
                         ctgcACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTa                                                                                                                           >  2:1014864/1‑150 (MQ=255)
                               tctGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCt                                                                                                                     >  2:825536/1‑150 (MQ=255)
                                  gCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTAtt                                                                                                                  >  2:1371070/1‑150 (MQ=255)
                                  gCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTAtt                                                                                                                  >  2:715354/1‑150 (MQ=255)
                                   cccTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTAttt                                                                                                                 <  1:485096/150‑1 (MQ=255)
                                          ccGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGc                                                                                                          <  1:1067754/150‑1 (MQ=255)
                                          ccGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGc                                                                                                          <  1:983499/150‑1 (MQ=255)
                                             ttAATCCGATAACGAAAGTCAACTGGCGTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGATGATCAGCATCGCAGCTCGCCGCTGACCTGGTTCGGCGCGTTAGGTTCTATTTTACTGGCGAt                                                                                                       <  2:459148/150‑1 (MQ=255)
                                               aTTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATgg                                                                                                     <  2:678948/150‑1 (MQ=255)
                                                ttCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATg                                                                                                      <  1:87893/148‑1 (MQ=255)
                                                ttCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATg                                                                                                      >  2:87893/1‑148 (MQ=255)
                                                  ccGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCa                                                                                                  >  2:30341/1‑150 (MQ=255)
                                                   cGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAg                                                                                                 >  1:147544/1‑150 (MQ=255)
                                                        aCGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTAtt                                                                                                                  >  2:1048840/1‑128 (MQ=255)
                                                        aCGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTAtt                                                                                                                  <  1:1048840/128‑1 (MQ=255)
                                                        aCGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGc                                                                                            >  1:401049/1‑150 (MQ=255)
                                                         cGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCa                                                                                           >  2:1207515/1‑150 (MQ=255)
                                                         cGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCa                                                                                           <  1:482085/150‑1 (MQ=255)
                                                         cGAAAGTCAACTAGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTcc                                                                                              <  2:1187753/147‑1 (MQ=255)
                                                         cGAAAGTCAACTAGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTcc                                                                                              >  1:1187753/1‑147 (MQ=255)
                                                          gAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAg                                                                                          >  2:1520712/1‑150 (MQ=255)
                                                          caaaGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACCTCACCCAACGTCACCTGGCGCGAGCGATGGAACATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaata                                                                                          <  1:427804/149‑7 (MQ=255)
                                                           aaaGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAgg                                                                                         >  1:1054026/1‑150 (MQ=255)
                                                             aGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTAGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCg                                                                                       >  2:1002587/1‑150 (MQ=255)
                                                              gTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgg                                                                                      >  2:443736/1‑150 (MQ=255)
                                                              gTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgg                                                                                      <  1:814872/150‑1 (MQ=255)
                                                                cAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgggg                                                                                    >  2:940043/1‑150 (MQ=255)
                                                                cAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgggg                                                                                    <  1:1288114/150‑1 (MQ=255)
                                                                     ggCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggc                                                                               <  1:1160061/150‑1 (MQ=255)
                                                                     ggCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggc                                                                               <  2:521785/150‑1 (MQ=255)
                                                                      gCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcg                                                                              >  1:1044790/1‑150 (MQ=255)
                                                                       cTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTGACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcgg                                                                             <  1:130807/150‑1 (MQ=255)
                                                                       cTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcgg                                                                             >  1:1364474/1‑150 (MQ=255)
                                                                       cTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATggcgg                                                                             >  1:1372079/1‑150 (MQ=255)
                                                                         tCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCCCCTGGCGCGACCAATGCAAGATCAGCAGCGCAGCCCGCCGCAGCCCTGGCTCGCCGCGTAAGGATCAATTTTACTGGCGATGCCCAGACCCCAGCCGGGGATGccgccg                                                                           >  2:150188/1‑150 (MQ=255)
                                                                            gTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTg                                                                        <  1:957507/150‑1 (MQ=255)
                                                                              tATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGAc                                                                      <  2:140422/150‑1 (MQ=255)
                                                                                 ggCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCgg                                                                   >  2:450436/1‑150 (MQ=255)
                                                                                  gcgcACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGt                                                                  >  1:678882/1‑150 (MQ=255)
                                                                                   cgcACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTa                                                                 >  2:1330062/1‑150 (MQ=255)
                                                                                   cgcACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCGGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTa                                                                 <  1:1210820/150‑1 (MQ=255)
                                                                                    gcACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgggg                                                                                    <  2:1274586/130‑1 (MQ=255)
                                                                                    gcACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgggg                                                                                    >  1:1274586/1‑130 (MQ=255)
                                                                                      aCGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTacac                                                              <  2:92726/150‑1 (MQ=255)
                                                                                             ctcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCggg                                                       <  2:10862/150‑1 (MQ=255)
                                                                                               cccACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGaa                                                     >  2:1259327/1‑150 (MQ=255)
                                                                                                       aCCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAg                                             <  1:825536/150‑1 (MQ=255)
                                                                                                           aaCGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGga                                         <  1:1371070/150‑1 (MQ=255)
                                                                                                                cACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATc                                    >  2:737460/1‑150 (MQ=255)
                                                                                                                cACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATc                                    <  1:1330062/150‑1 (MQ=255)
                                                                                                                  ccTGGCGCGAGCGATGGAAGATCAGCAGCGCAACGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAg                                  >  1:1029567/1‑150 (MQ=255)
                                                                                                                   cTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGt                                 <  1:14225/150‑1 (MQ=255)
                                                                                                                          gagcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCAccc                          <  1:1123054/150‑1 (MQ=255)
                                                                                                                            gcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCAGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAg                        <  2:147544/150‑1 (MQ=255)
                                                                                                                                     ggtCAGCGGCGCGCTGCGCTGCTGATCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGaa               >  1:945739/3‑150 (MQ=255)
                                                                                                                                      gtCAGCGGCGCGCTGCGCTGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGACCGCGTCAGGGGATGATCAGTTTCCCCCAGCAAAATGAAc              >  1:359162/2‑150 (MQ=255)
                                                                                                                                        cagcagCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAg            >  1:735099/1‑150 (MQ=255)
                                                                                                                                        cagcagCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAg            >  1:29591/1‑150 (MQ=255)
                                                                                                                                          gcagcgcagcGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGa          <  2:1044790/150‑1 (MQ=255)
                                                                                                                                             gcgcagcgcGCCGCTGACCTGGGTCGGCGCGTTAGGTTATATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACACAGCAAAATGAACAGGAAGc       >  2:95424/1‑150 (MQ=255)
                                                                                                                                             gcgcagcgcGCCGCTGACCTGGGTCGGCGCGTTAGGTTATATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTCGCGCGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGCACAGGAAGc       >  2:89910/1‑150 (MQ=255)
                                                                                                                                              cgcagcgcGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCg      >  2:1131563/1‑150 (MQ=255)
                                                                                                                                               gcagcgcGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCgg     >  2:221160/1‑150 (MQ=255)
                                                                                                                                                 agcgcGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGAGCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGAc   <  2:1372079/150‑1 (MQ=255)
                                                                                                                                                  gcgcgcCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGAcc  >  2:1093277/1‑150 (MQ=255)
                                                                                                                                                  gcgcgcCGCTGACCTAGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGAcc  <  1:1002587/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGCTTTCTTTGGCGGCAACGTGGTGCTGCACTCTGCCCTGTTCCGTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACC  >  NC_000913/2616414‑2616709

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