Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,007,368 C→G 6.1% A369A (GCC→GCG ygeW → putative carbamoyltransferase YgeW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,007,3680CG 6.1% 220.7 / 4.6 65A369A (GCC→GCGygeWputative carbamoyltransferase YgeW
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base G (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

TGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCC  >  NC_000913/3007224‑3007515
                                                                                                                                                |                                                                                                                                                   
tGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGGAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCgat                                                                                                                                                <  2:697522/150‑1 (MQ=255)
       gAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTg                                                                                                                                         <  1:354644/150‑1 (MQ=255)
        aaGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCGGCgatgatataag                                                                                                                                        <  2:749202/150‑6 (MQ=255)
        aaGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCGGCgatgatataag                                                                                                                                        <  1:38453/150‑6 (MQ=255)
        aaGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGt                                                                                                                                        <  1:5112/150‑1 (MQ=255)
          gCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTcc                                                                                                                                      <  1:834461/150‑1 (MQ=255)
                      cACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGcc                                                                                                                          >  2:1258753/1‑150 (MQ=255)
                        ctgcctgcCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCaa                                                                                                                        >  1:1024730/1‑150 (MQ=255)
                         tgcctgccAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCaaa                                                                                                                       >  1:465903/1‑150 (MQ=255)
                          gcctgccAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAc                                                                                                                      >  2:872050/1‑150 (MQ=255)
                             tgccAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAg                                                                                                                   >  2:1071888/1‑150 (MQ=255)
                                cAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTg                                                                                                                <  2:370384/150‑1 (MQ=255)
                                cAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTg                                                                                                                <  2:354298/150‑1 (MQ=255)
                                 aGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGc                                                                                                               >  1:935913/1‑150 (MQ=255)
                                  gCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCa                                                                                                              >  2:442197/1‑150 (MQ=255)
                                     gATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTg                                                                                                           >  1:294274/1‑150 (MQ=255)
                                       tatCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCt                                                                                                         <  2:38741/150‑1 (MQ=255)
                                        atCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTc                                                                                                        <  1:479096/150‑1 (MQ=255)
                                          cAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGa                                                                                                      >  2:1021982/1‑150 (MQ=255)
                                           aGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAg                                                                                                     >  2:36232/1‑150 (MQ=255)
                                           aGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAg                                                                                                     <  2:979166/150‑1 (MQ=255)
                                           aGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAg                                                                                                     >  1:3623/1‑150 (MQ=255)
                                                cgTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAct                                                                                                <  1:585517/150‑1 (MQ=255)
                                                cgTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAct                                                                                                <  1:690186/150‑1 (MQ=255)
                                                      cTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaa                                                                                          <  1:1258753/150‑1 (MQ=255)
                                                      cTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaa                                                                                          >  1:1049589/1‑150 (MQ=255)
                                                      cTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCGGaa                                                                                          >  2:377402/1‑150 (MQ=255)
                                                       tGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaaa                                                                                         >  1:1111999/1‑150 (MQ=255)
                                                       tGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaaa                                                                                         <  2:1049589/150‑1 (MQ=255)
                                                         tAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaaga                                                                                       >  2:501545/1‑150 (MQ=255)
                                                                ggTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCaa                                                                                >  2:689646/1‑150 (MQ=255)
                                                                  tGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGa                                                                              >  1:91462/1‑150 (MQ=255)
                                                                  tGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGa                                                                              >  1:1374235/1‑150 (MQ=255)
                                                                     aGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcg                                                                           >  2:524504/1‑150 (MQ=255)
                                                                        gACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAACCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGt                                                                        >  1:419506/1‑150 (MQ=255)
                                                                          cTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGa                                                                      >  1:42183/1‑150 (MQ=255)
                                                                            gAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGaaa                                                                    >  1:663483/1‑150 (MQ=255)
                                                                               ggCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATaa                                                                 >  2:133729/1‑150 (MQ=255)
                                                                                 cGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAAcc                                                               <  1:905104/150‑1 (MQ=255)
                                                                                       cGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTTCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCgg                                                         <  1:524504/150‑1 (MQ=255)
                                                                                        gAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCATGTACCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCggg                                                        >  2:1143871/1‑150 (MQ=255)
                                                                                           aaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCg                                                     <  2:527930/150‑1 (MQ=255)
                                                                                                 cGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGc                                               <  2:206822/150‑1 (MQ=255)
                                                                                                 cGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGc                                               >  1:310256/1‑150 (MQ=255)
                                                                                                  gTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCt                                              <  1:1021982/150‑1 (MQ=255)
                                                                                                      cGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTc                                          <  1:46439/150‑1 (MQ=255)
                                                                                                          acctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGcc                                      <  1:1071888/150‑1 (MQ=255)
                                                                                                           cctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGccc                                     <  1:36232/150‑1 (MQ=255)
                                                                                                            ctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGcccg                                    <  2:812219/150‑2 (MQ=255)
                                                                                                             tacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCctc                                   >  2:774968/1‑150 (MQ=255)
                                                                                                              acAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCctct                                  <  2:663483/150‑1 (MQ=255)
                                                                                                               cAAATAAGCAAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTAGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTa                                 >  2:81260/1‑150 (MQ=255)
                                                                                                               cAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGcccgagg                                 <  1:872050/150‑5 (MQ=255)
                                                                                                                    aagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTcc                            <  1:501545/150‑1 (MQ=255)
                                                                                                                      gCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCag                          >  2:675624/1‑150 (MQ=255)
                                                                                                                         aGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAgg                       <  1:1161568/150‑1 (MQ=255)
                                                                                                                               aaGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAATTACGCCAAACCAGGTGCACTGTTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaaaa                 >  1:945227/1‑150 (MQ=255)
                                                                                                                                aGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAag                <  2:81444/150‑1 (MQ=255)
                                                                                                                                   cTTATATCATCGCCGCGATGATAATGTCCCGTAAATAAGCCAAACCAGGTGAAATGCTCGAGAAAATGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggat             >  2:1009382/1‑150 (MQ=255)
                                                                                                                                    ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGCTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggata            <  1:418758/150‑1 (MQ=255)
                                                                                                                                     aggatcatcGCGGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatag           >  1:1289802/4‑150 (MQ=255)
                                                                                                                                     aggatcatcGCGGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatag           >  1:667825/4‑150 (MQ=255)
                                                                                                                                       tatCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTTAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatagga         <  1:442197/150‑1 (MQ=255)
                                                                                                                                         tcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAggataggata       <  2:42183/150‑1 (MQ=255)
                                                                                                                                          catcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGatat      >  2:1426963/1‑150 (MQ=255)
                                                                                                                                              gccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGCTCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTcc  >  2:554404/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCC  >  NC_000913/3007224‑3007515

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: