Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 33,732 T→C D972D (GAT→GAC carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091333,7320TC81.1% 560.8 / 138.3 280D972D (GAT→GACcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base C (134/93);  total (164/116)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-01

Only 100 of 282 total aligned reads displayed.

CACGGTCGTGCGCTGCTTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTCCGCACATTCAGGA  >  NC_000913/33640‑33824
                                                                                            |                                                                                            
cacGGTCGTGCGCTGCTTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGAt                                                                                              <  1:5900663/93‑1 (MQ=255)
       gTGCGCTGCTTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGAccc                                                                                       <  1:7334419/93‑1 (MQ=255)
       gTGCGCTGCTTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGAccc                                                                                       <  1:529484/93‑1 (MQ=255)
        tGCGCTGCTTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCca                                                                                      <  1:3550668/93‑1 (MQ=255)
        tGCGCTGCTTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCca                                                                                      <  2:3239788/93‑1 (MQ=255)
           gctgctTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCcacgg                                                                                   >  2:7527656/1‑93 (MQ=255)
           gctgctTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCcacgg                                                                                   >  1:4323017/1‑93 (MQ=255)
            ctgctTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCcacggc                                                                                  >  2:6985615/1‑93 (MQ=255)
            ctgctTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCcacggc                                                                                  >  2:5388692/1‑93 (MQ=255)
            ctgctTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCcacggc                                                                                  >  2:2915014/1‑93 (MQ=255)
             tgctTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGACCcacggca                                                                                 >  2:5048319/1‑93 (MQ=255)
             tgctTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCcacggca                                                                                 <  2:7194156/93‑1 (MQ=255)
             tgctTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCcacggca                                                                                 >  2:6764974/1‑93 (MQ=255)
              gctTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGACCcacggcac                                                                                <  2:5192027/93‑1 (MQ=255)
               ctTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCcacggcacg                                                                               >  1:3455540/1‑93 (MQ=255)
                tttCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCcacggcacgg                                                                              >  2:1204268/1‑93 (MQ=255)
                tttCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCcacggcacgg                                                                              >  1:1256941/1‑93 (MQ=255)
                   ccGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGa                                                                           <  1:4790004/93‑1 (MQ=255)
                    cGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGAt                                                                          <  1:3633272/93‑1 (MQ=255)
                     gTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGAtt                                                                         <  1:291218/93‑1 (MQ=255)
                       gcgcGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATtgt                                                                       <  1:5145899/93‑1 (MQ=255)
                          cgAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCt                                                                    <  2:5726197/93‑1 (MQ=255)
                           gAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTg                                                                   >  1:6948058/1‑93 (MQ=255)
                           gAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTg                                                                   <  1:6399113/93‑1 (MQ=255)
                            aaGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTgg                                                                  >  2:811424/1‑93 (MQ=255)
                             aGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTggg                                                                 <  1:6268496/93‑1 (MQ=255)
                              ggCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGc                                                                >  1:905392/1‑93 (MQ=255)
                              ggCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGc                                                                >  1:2319691/1‑93 (MQ=255)
                               gCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCg                                                               <  2:2161373/93‑1 (MQ=255)
                                 gATAAAGAACGCGTGGTGGGCCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGaa                                                             <  1:6602526/93‑1 (MQ=255)
                                 gATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGaa                                                             >  2:4548876/1‑93 (MQ=255)
                                  aTAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAg                                                            >  2:1612909/1‑93 (MQ=255)
                                   ttttgAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGc                                                           >  1:3937233/5‑93 (MQ=255)
                                    aaaGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCa                                                          >  2:4089530/1‑93 (MQ=255)
                                     aaGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAg                                                         <  1:7454249/93‑1 (MQ=255)
                                     aaGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAg                                                         <  2:7085493/93‑1 (MQ=255)
                                      aGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAgg                                                        >  2:6043671/1‑93 (MQ=255)
                                         aCGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTAt                                                     >  1:769549/1‑93 (MQ=255)
                                         aCGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTAt                                                     >  1:617529/1‑93 (MQ=255)
                                         aCGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTAt                                                     >  1:4456335/1‑93 (MQ=255)
                                          cgcgTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATc                                                    >  2:6866865/1‑93 (MQ=255)
                                           gcgTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCa                                                   >  1:1116181/1‑93 (MQ=255)
                                           gcgTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCa                                                   >  1:6784700/1‑93 (MQ=255)
                                           gcgTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCa                                                   >  2:2605714/1‑93 (MQ=255)
                                           gcgTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCa                                                   <  2:118704/93‑1 (MQ=255)
                                             gtggtgGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAAc                                                 >  2:3472424/1‑93 (MQ=255)
                                             gtggtgGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAAc                                                 <  1:147312/93‑1 (MQ=255)
                                              tggtggACCTGGCGGCAAAACTGCTGAAAGAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAAcc                                                <  1:5969846/93‑1 (MQ=255)
                                               ggtggACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGGATCAAccc                                               >  1:1397210/1‑93 (MQ=255)
                                                gtggACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCcg                                              >  1:6431821/1‑93 (MQ=255)
                                                 tggACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCcgc                                             >  1:1613915/1‑93 (MQ=255)
                                                   gACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGt                                           >  2:6205258/1‑93 (MQ=255)
                                                    aCCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTc                                          <  1:5546131/93‑1 (MQ=255)
                                                    aCCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTc                                          >  1:2138063/1‑93 (MQ=255)
                                                      cTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTg                                        >  1:1093648/1‑93 (MQ=255)
                                                      cTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTg                                        >  1:2198196/1‑93 (MQ=255)
                                                       tGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGCAGCAGGTATCAACCCGCGTCTgg                                       >  1:3325522/1‑93 (MQ=255)
                                                       tGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTgg                                       >  2:3807310/1‑93 (MQ=255)
                                                        ggcggcAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGt                                      <  1:2542022/93‑1 (MQ=255)
                                                        ggcggcAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGt                                      >  2:4452924/1‑93 (MQ=255)
                                                           ggcAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTaaa                                   >  1:1953211/1‑93 (MQ=255)
                                                            gcAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAAc                                  <  2:1884048/93‑1 (MQ=255)
                                                            gcAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAAc                                  <  2:2462868/93‑1 (MQ=255)
                                                             cAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACa                                 >  2:267556/1‑93 (MQ=255)
                                                             cAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACa                                 >  2:7175056/1‑93 (MQ=255)
                                                             cAAAACTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACa                                 <  2:6495549/93‑1 (MQ=255)
                                                               aaaCTGCTGAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAg                               >  1:5443885/1‑93 (MQ=255)
                                                               aaaCTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAg                               >  1:3026927/1‑93 (MQ=255)
                                                                 aCTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGt                             >  2:4960667/1‑93 (MQ=255)
                                                                 aCTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGt                             <  2:5877721/93‑1 (MQ=255)
                                                                 aCTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGt                             <  2:7311742/93‑1 (MQ=255)
                                                                 aCTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGt                             <  1:4559498/93‑1 (MQ=255)
                                                                 aCTGCTGAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGt                             <  1:4408398/93‑1 (MQ=255)
                                                                   tgctgAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGc                           >  2:5264294/1‑93 (MQ=255)
                                                                   tgctgAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGc                           <  1:4494964/93‑1 (MQ=255)
                                                                    gctgAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCa                          <  1:6472099/93‑1 (MQ=255)
                                                                     ctgAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCAt                         >  2:6378134/1‑93 (MQ=255)
                                                                      tgAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATg                        >  2:3735973/1‑93 (MQ=255)
                                                                       gAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGa                       <  2:1689620/93‑1 (MQ=255)
                                                                       gAAACAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGa                       >  1:3723218/1‑93 (MQ=255)
                                                                        aaaCAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGaa                      <  1:5852642/93‑1 (MQ=255)
                                                                         aaCAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAg                     <  1:6914895/93‑1 (MQ=255)
                                                                          aCAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAgg                    >  2:2783776/1‑93 (MQ=255)
                                                                          aCAGGGCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAgg                    >  1:7222235/1‑93 (MQ=255)
                                                                           cAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGc                   <  2:3549527/93‑1 (MQ=255)
                                                                             gggCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCg                 >  1:7072957/1‑93 (MQ=255)
                                                                             gggCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCg                 >  1:4610204/1‑93 (MQ=255)
                                                                             gggCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCg                 >  1:2442626/1‑93 (MQ=255)
                                                                              ggCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGt                <  2:4487692/93‑1 (MQ=255)
                                                                              ggCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGt                >  2:6239102/1‑93 (MQ=255)
                                                                               gCTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTc               >  2:5221918/1‑93 (MQ=255)
                                                                                cTTCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTcc              >  2:1025513/1‑93 (MQ=255)
                                                                                 ttCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTCCg             >  2:1990720/1‑93 (MQ=255)
                                                                                  tCGAGCTGGACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTCCGc            >  2:6307660/1‑93 (MQ=255)
                                                                                   cGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTCCGca           >  1:1678472/1‑93 (MQ=255)
                                                                                     aGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTCCGcaca         >  2:3063360/1‑93 (MQ=255)
                                                                                       cTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTCCGCACAtt       <  2:2326158/93‑1 (MQ=255)
                                                                                         ggACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTCCGCACATTCa     >  2:54756/1‑93 (MQ=255)
                                                                                         ggACGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTCCGCACATTCa     <  1:2521974/93‑1 (MQ=255)
                                                                                            tGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTCCGCACATTCAGGa  <  2:5823019/93‑1 (MQ=255)
                                                                                            |                                                                                            
CACGGTCGTGCGCTGCTTTCCGTGCGCGAAGGCGATAAAGAACGCGTGGTGGACCTGGCGGCAAAACTGCTGAAACAGGGCTTCGAGCTGGATGCGACCCACGGCACGGCGATTGTGCTGGGCGAAGCAGGTATCAACCCGCGTCTGGTAAACAAGGTGCATGAAGGCCGTCCGCACATTCAGGA  >  NC_000913/33640‑33824

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