Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,005,399 +C coding (1044/1092 nt) ygeW → putative carbamoyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,005,3981.C99.5% 2182.6 / 0.8 549coding (1043/1092 nt)ygeWputative carbamoyl transferase
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (246/299);  total (246/303)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 619 total aligned reads displayed.

TAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATC  >  NC_000913/3005303‑3005488
                                                                                                 |                                                                                           
 aaaGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATc                                                                                            >  1:4907532/1‑96 (MQ=255)
  aagaagGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcc                                                                                           <  1:9438296/96‑2 (MQ=255)
    gaagGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcctg                                                                                         >  2:13624278/1‑93 (MQ=255)
     aagGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcctgt                                                                                        >  2:10737096/1‑92 (MQ=255)
     aagGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcctgt                                                                                        >  1:14135544/1‑92 (MQ=255)
       ggTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTcc                                                                                      <  2:5106169/96‑1 (MQ=255)
       ggTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTcc                                                                                      <  2:1556440/96‑1 (MQ=255)
       ggTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTcc                                                                                      <  2:1452010/96‑1 (MQ=255)
          gAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGt                                                                                   >  2:9348878/1‑96 (MQ=255)
          gAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGt                                                                                   <  2:11951414/96‑1 (MQ=255)
            ggTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaa                                                                                 >  1:13868349/2‑96 (MQ=255)
             gTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaaa                                                                                >  2:13857971/1‑96 (MQ=255)
              tgccTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGGCCCGTAAAt                                                                               <  1:5018260/93‑1 (MQ=255)
              tGACTGAAGGCGTATTCGCAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTACAt                                                                               >  1:10175668/1‑96 (MQ=255)
               gACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATa                                                                              >  2:3813452/1‑96 (MQ=255)
                aCTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATAc                                                                             <  1:13230/96‑1 (MQ=255)
                  tGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGc                                                                           >  1:11357011/1‑96 (MQ=255)
                  tGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGc                                                                           >  1:7811118/1‑96 (MQ=255)
                   gAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGcc                                                                          >  2:8282678/1‑96 (MQ=255)
                   gAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGcc                                                                          <  2:8735594/96‑1 (MQ=255)
                    aaGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCa                                                                         >  2:1950400/1‑96 (MQ=255)
                    aaGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCa                                                                         <  2:1329696/96‑1 (MQ=255)
                    aaGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCa                                                                         >  2:433590/1‑96 (MQ=255)
                     aGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCaa                                                                        <  1:3220327/96‑1 (MQ=255)
                      ggCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCaaa                                                                       <  2:10053549/96‑1 (MQ=255)
                       gCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAc                                                                      <  2:12881184/96‑1 (MQ=255)
                        cGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAcc                                                                     >  1:10886561/1‑96 (MQ=255)
                              cGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGc                                                               <  1:9747862/96‑1 (MQ=255)
                              cGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGc                                                               >  1:6332528/1‑96 (MQ=255)
                               gAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCa                                                              >  2:6550671/1‑96 (MQ=255)
                               gAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCa                                                              <  2:2485878/96‑1 (MQ=255)
                               gAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCa                                                              <  2:3093603/96‑1 (MQ=255)
                                aaaaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCAc                                                             >  2:7743873/1‑96 (MQ=255)
                                 aaaaTACCGTATCGCTACCTACAAAGAAGCCAGGTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACt                                                            >  2:14055231/1‑96 (MQ=255)
                                 aaaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCGGGTGCACt                                                            <  1:12947582/96‑1 (MQ=255)
                                  gaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTg                                                           >  1:11768265/2‑96 (MQ=255)
                                   aaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGc                                                          <  2:4988333/96‑1 (MQ=255)
                                      aCCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCtgg                                                       <  1:9311597/96‑3 (MQ=255)
                                      aCCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCg                                                       >  2:2225331/1‑96 (MQ=255)
                                        cGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAg                                                     >  2:6934765/1‑96 (MQ=255)
                                        cGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACACCAAACCAGGTGCACTGCTCGAg                                                     >  1:7594923/1‑96 (MQ=255)
                                           aTCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCaa                                                  <  1:13639336/96‑1 (MQ=255)
                                           aTCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCaa                                                  >  1:4706633/1‑96 (MQ=255)
                                            tCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAc                                                 <  2:2658199/96‑1 (MQ=255)
                                             cGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAct                                                <  1:10284887/96‑1 (MQ=255)
                                              gCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActg                                               >  1:9858127/1‑96 (MQ=255)
                                               ctacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgc                                              >  2:6369511/1‑96 (MQ=255)
                                               ctacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgc                                              >  1:10023667/1‑96 (MQ=255)
                                               ctacctacAAAGAAGCCAGCTGGAAGCCGTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgc                                              >  2:11401808/1‑96 (MQ=255)
                                                tacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgct                                             <  1:8011268/96‑1 (MQ=255)
                                                tacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgct                                             >  1:3951037/1‑96 (MQ=255)
                                                 acctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATCCGCCAAACCAGGTGCACTGCTCGAGCAActgctg                                            <  2:11139288/96‑1 (MQ=255)
                                                 acctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgctg                                            >  1:2543428/1‑96 (MQ=255)
                                                 acctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgctg                                            <  2:10306839/96‑1 (MQ=255)
                                                   ctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATCCGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaa                                          <  2:2846518/96‑1 (MQ=255)
                                                   ctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaa                                          >  1:863436/1‑96 (MQ=255)
                                                    tacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaaa                                         >  2:8071933/1‑96 (MQ=255)
                                                     acAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaag                                        >  2:4074017/1‑96 (MQ=255)
                                                        aagaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaagaag                                     >  1:2818718/1‑96 (MQ=255)
                                                        aagaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaagaag                                     >  2:485272/1‑96 (MQ=255)
                                                         agaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgc                                    <  2:6153797/96‑1 (MQ=255)
                                                         agaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgc                                    <  1:13052419/96‑1 (MQ=255)
                                                         agaagCCAGCTGGAAGCATTATATCATCGCCGCGATGATCCTGTCCCGTAAATCCGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaagaatc                                    <  1:3419979/96‑3 (MQ=255)
                                                          gacgCCAGCTGGAAGCCTTATATCATCGACGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgcg                                   <  2:13127489/93‑1 (MQ=255)
                                                           aagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgcgc                                  <  1:3712478/96‑1 (MQ=255)
                                                            agCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCa                                 >  1:1876845/1‑96 (MQ=255)
                                                            agCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCa                                 >  1:2240875/1‑96 (MQ=255)
                                                            agCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCTTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCa                                 <  2:12216822/96‑1 (MQ=255)
                                                              ccAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAg                               <  1:4288419/96‑1 (MQ=255)
                                                              ccAGCGGGCAGTCTTATATCATCGCGGGGATGAGCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAg                               <  1:2997120/96‑1 (MQ=255)
                                                                 gCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAc                            <  2:12520514/96‑1 (MQ=255)
                                                                 gCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAc                            <  2:10238739/96‑1 (MQ=255)
                                                                    ggAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgcg                         >  1:8224574/1‑96 (MQ=255)
                                                                      aaGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTg                       <  2:7493766/96‑1 (MQ=255)
                                                                      aaGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTg                       <  1:6545130/96‑1 (MQ=255)
                                                                       aGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGa                      >  1:12132788/1‑96 (MQ=255)
                                                                        gCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGAGaa                     >  2:3555378/1‑96 (MQ=255)
                                                                         ccTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGaaa                    >  2:9090029/1‑96 (MQ=255)
                                                                          cTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAAt                   >  1:12688882/1‑96 (MQ=255)
                                                                           ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATa                  <  2:5501784/96‑1 (MQ=255)
                                                                            tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATaa                 >  1:204962/1‑96 (MQ=255)
                                                                            tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATaa                 <  1:1450473/96‑1 (MQ=255)
                                                                            tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATaa                 >  2:5690034/1‑96 (MQ=255)
                                                                             atatCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAAc                >  1:1049187/1‑96 (MQ=255)
                                                                               atcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAAccc              >  2:10525650/1‑96 (MQ=255)
                                                                                tcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAAcccc             <  2:985034/96‑1 (MQ=255)
                                                                                 catcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCt            <  2:8600551/96‑1 (MQ=255)
                                                                                  atcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTc           <  1:5633322/96‑1 (MQ=255)
                                                                                  atcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTc           >  2:4600526/1‑96 (MQ=255)
                                                                                  atcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTc           <  1:9105723/96‑1 (MQ=255)
                                                                                    cgccgcGATGATCCTGTCCCGTAGATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCgg         <  1:12531107/96‑1 (MQ=255)
                                                                                    cgccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCgg         >  2:3319506/1‑96 (MQ=255)
                                                                                     gccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCggg        <  2:7109466/96‑1 (MQ=255)
                                                                                      ccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAACTAACCCCTCGGGc       <  1:10994203/96‑1 (MQ=255)
                                                                                        gcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCg     >  1:6006332/1‑96 (MQ=255)
                                                                                        gcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCg     >  2:10944274/1‑96 (MQ=255)
                                                                                          gatgatcCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGa  >  1:709038/1‑95 (MQ=255)
                                                                                          gatgatcCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGAt   >  1:537773/1‑96 (MQ=255)
                                                                                          gatgatcCTGTCCCGTAAATACGCCAAACCAGGAGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGAt   >  1:8956231/1‑96 (MQ=255)
                                                                                           atgatcCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATc  <  1:6956756/96‑1 (MQ=255)
                                                                                                 |                                                                                           
TAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATC  >  NC_000913/3005303‑3005488

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