Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,009,141 G→A 5.2% Q38Q (CAG→CAA) ‡ PP_3533 → ornithine cyclodeaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,009,1410GA 5.2% 1221.6 / 33.9 383Q38Q (CAG→CAA) ‡PP_3533ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base G (174/188);  new base A (10/10);  total (185/198)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 396 total aligned reads displayed.

CTGCTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGGCCTT  >  NC_002947/4008992‑4009290
                                                                                                                                                     |                                                                                                                                                     
ctgctgGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAg                                                                                                                                                       <  2:8125211/150‑1 (MQ=255)
  gctgGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGc                                                                                                                                                     <  1:148551/150‑1 (MQ=255)
  gctgGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGc                                                                                                                                                     <  2:7962437/150‑1 (MQ=255)
     gggCAAGCGCACACTCACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCCTGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCCAACGCTGGCAGGCGtt                                                                                                                                                  <  2:8229066/150‑1 (MQ=255)
     gggCAAGCGCACAATACCAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGtt                                                                                                                                                  <  2:2635375/150‑1 (MQ=255)
     gggCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGtt                                                                                                                                                  <  1:2486773/150‑1 (MQ=255)
      ggCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGttt                                                                                                                                                 <  1:4728479/150‑1 (MQ=255)
       gCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTg                                                                                                                                                <  2:8636796/150‑1 (MQ=255)
        cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCATGCGTTTGa                                                                                                                                               >  1:6682979/1‑150 (MQ=255)
        cAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGa                                                                                                                                               >  2:1700861/1‑150 (MQ=255)
          aGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACa                                                                                                                                             >  2:6061142/1‑150 (MQ=255)
          aGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACa                                                                                                                                             <  2:7091058/150‑1 (MQ=255)
           gcgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACaa                                                                                                                                            <  1:7832306/150‑1 (MQ=255)
           gcgcACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCCTGCCAGCGTTTGaa                                                                                                                                              <  1:1770310/148‑2 (MQ=255)
                      cAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgcgcg                                                                                                                                 <  2:780996/150‑1 (MQ=255)
                      cAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCcgcgcgcg                                                                                                                                 >  1:6493774/1‑150 (MQ=255)
                       aaGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGAGCCGTTTATTGGCGATCTTGCCGCTGCCCTGCGGGACGAGTTCAAACGCTGGCAGGCGTTTGACAAGTCAGCGCGCGt                                                                                                                                >  2:3826476/1‑150 (MQ=255)
                            gCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgcca                                                                                                                           >  2:6482754/1‑150 (MQ=255)
                            gCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgcca                                                                                                                           <  2:5174913/150‑1 (MQ=255)
                             cAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccag                                                                                                                          <  1:4095249/150‑1 (MQ=255)
                             cAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccag                                                                                                                          <  1:5441123/150‑1 (MQ=255)
                              aaGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccagc                                                                                                                         <  1:2726499/150‑1 (MQ=255)
                               agagCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTgccagcc                                                                                                                        <  2:4263102/150‑1 (MQ=255)
                                    aTGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTTTTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCg                                                                                                                   >  2:4428074/1‑150 (MQ=255)
                                       aCGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGaa                                                                                                                <  1:4328337/150‑1 (MQ=255)
                                       aCGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGaa                                                                                                                >  1:5520162/1‑150 (MQ=255)
                                          tATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTg                                                                                                             >  1:3632451/1‑150 (MQ=255)
                                           aTTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTgg                                                                                                            <  2:4529945/150‑1 (MQ=255)
                                           aTTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTgg                                                                                                            <  1:2656552/150‑1 (MQ=255)
                                           aTTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTgg                                                                                                            <  1:5470514/150‑1 (MQ=255)
                                             ttCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGc                                                                                                          >  1:8478514/1‑150 (MQ=255)
                                               cATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGt                                                                                                        <  2:3265529/150‑1 (MQ=255)
                                               cATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGt                                                                                                        <  1:4236422/150‑1 (MQ=255)
                                                aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTTGAAGGCGTTTGACAAGTACGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTg                                                                                                       >  1:4412546/1‑150 (MQ=255)
                                                aTTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTg                                                                                                       >  2:4827043/1‑150 (MQ=255)
                                                 ttGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGa                                                                                                      <  2:385269/150‑1 (MQ=255)
                                                   gATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATc                                                                                                    <  1:1036538/150‑1 (MQ=255)
                                                    aTGTTCCAACCATGTCGGACTTGGTGCATGACCTTGGCGTAGCGCCGTTTAGTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCg                                                                                                   >  2:3675511/1‑150 (MQ=255)
                                                      gTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGaa                                                                                                 <  2:3958110/150‑1 (MQ=255)
                                                       ttCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAAc                                                                                                >  1:890034/1‑150 (MQ=255)
                                                           aaCCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGAt                                                                                            <  1:3355243/150‑1 (MQ=255)
                                                           aaCCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGAt                                                                                            >  1:7980652/1‑150 (MQ=255)
                                                           aaCCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGAt                                                                                            >  1:8077439/1‑150 (MQ=255)
                                                             ccATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGc                                                                                          <  2:7381612/150‑1 (MQ=255)
                                                                    ggACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGc                                                                                   >  1:6987586/1‑150 (MQ=255)
                                                                    ggACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGc                                                                                   >  2:3041413/1‑150 (MQ=255)
                                                                          ggTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaa                                                                             >  2:3558759/1‑150 (MQ=255)
                                                                            tGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaaaa                                                                           <  2:3227369/150‑1 (MQ=255)
                                                                            tGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACaaaa                                                                           <  2:4183024/150‑1 (MQ=255)
                                                                              cATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGc                                                                         >  1:1196629/1‑150 (MQ=255)
                                                                                tGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCg                                                                       >  2:5060196/1‑150 (MQ=255)
                                                                                  aCATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCt                                                                     <  2:7654032/150‑1 (MQ=255)
                                                                                     ttGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATg                                                                  >  2:2485064/1‑150 (MQ=255)
                                                                                     ttGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATg                                                                  >  1:973118/1‑150 (MQ=255)
                                                                                       ggCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGcc                                                                >  2:5496634/1‑150 (MQ=255)
                                                                                        gCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCt                                                               >  2:7438690/1‑150 (MQ=255)
                                                                                         cGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCTGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCtt                                                              <  1:2839861/150‑1 (MQ=255)
                                                                                         cGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCtt                                                              <  1:1570756/150‑1 (MQ=255)
                                                                                             gcgcCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTTGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCaaa                                                          <  1:5496617/150‑1 (MQ=255)
                                                                                             gcgcCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCaaa                                                          <  1:7703479/150‑1 (MQ=255)
                                                                                                ccGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATAc                                                       <  2:1511351/150‑1 (MQ=255)
                                                                                                  gTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGt                                                     >  2:5352832/1‑150 (MQ=255)
                                                                                                      aTTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAAc                                                 <  1:3150776/150‑1 (MQ=255)
                                                                                                       ttGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACg                                                <  1:1982907/150‑1 (MQ=255)
                                                                                                       ttGGCGAGCTGGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACg                                                <  2:7922674/150‑1 (MQ=255)
                                                                                                         ggCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCTGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGc                                              <  2:2839597/150‑1 (MQ=255)
                                                                                                          gCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGcc                                             >  1:5258784/1‑150 (MQ=255)
                                                                                                            gAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCAc                                           >  1:8269684/1‑150 (MQ=255)
                                                                                                                ttGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCgg                                       >  1:2317581/1‑150 (MQ=255)
                                                                                                                ttGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCgg                                       >  1:7355688/1‑150 (MQ=255)
                                                                                                                ttGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCgg                                       >  2:1827900/1‑150 (MQ=255)
                                                                                                                 tGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGc                                      <  2:7235113/150‑1 (MQ=255)
                                                                                                                 ggCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTTGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGc                                      <  2:4669047/149‑1 (MQ=255)
                                                                                                                   ccGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCa                                    >  1:7122959/1‑150 (MQ=255)
                                                                                                                    cGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCaa                                   >  2:3545964/1‑150 (MQ=255)
                                                                                                                    cGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCaa                                   >  1:6807956/1‑150 (MQ=255)
                                                                                                                           cTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCa                            <  2:8163651/150‑1 (MQ=255)
                                                                                                                           cTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCa                            <  1:5403302/150‑1 (MQ=255)
                                                                                                                            tGCGGGACGATTTCAAGCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCAc                           <  1:4184524/150‑1 (MQ=255)
                                                                                                                            tGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCTGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCAc                           >  2:3427217/1‑150 (MQ=255)
                                                                                                                              cGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGc                         <  2:5723215/150‑1 (MQ=255)
                                                                                                                              cGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGc                         <  2:1397285/150‑1 (MQ=255)
                                                                                                                               gggACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCa                        >  2:7644870/1‑150 (MQ=255)
                                                                                                                                ggACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCaa                       <  1:6718657/150‑1 (MQ=255)
                                                                                                                                 gACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAAc                      <  2:8484392/150‑1 (MQ=255)
                                                                                                                                   cGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCTGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCt                    <  2:4756267/150‑1 (MQ=255)
                                                                                                                                    gATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTg                   >  2:1637581/1‑150 (MQ=255)
                                                                                                                                     ttgtCAAACGCCTGCCAGCGTTTGAAAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGc                  >  2:61286/4‑150 (MQ=255)
                                                                                                                                     aTTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGc                  >  2:2337340/1‑150 (MQ=255)
                                                                                                                                      tttCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCa                 >  1:7411441/1‑150 (MQ=255)
                                                                                                                                          aaaCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACg             <  2:4026799/150‑1 (MQ=255)
                                                                                                                                          aaaCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACg             >  1:6009585/1‑150 (MQ=255)
                                                                                                                                          aaaCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACg             >  1:5904620/1‑150 (MQ=255)
                                                                                                                                           aaCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACgg            <  1:3899790/150‑1 (MQ=255)
                                                                                                                                            aCGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGt           >  1:836146/1‑150 (MQ=255)
                                                                                                                                             cGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTg          <  2:3126206/150‑1 (MQ=255)
                                                                                                                                               cTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGAt        <  2:6333141/150‑1 (MQ=255)
                                                                                                                                                tGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATg       <  1:6269296/150‑1 (MQ=255)
                                                                                                                                                cggCAGGCGTATGACAAGTCCGCCCGCGCTGCCCGCCACTCGGAAGTGGGCGTCATCGAACTGACGCCGGTCGCCGACAACCGCCGCGCGGCCTTCCAAGACCTCAACGGCCACCCGGCCAATCCTGCACGCAACCTGCATACGGCGATg       <  2:411743/149‑1 (MQ=255)
                                                                                                                                                   cAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGGcc    >  1:6382113/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCTGGGCAAGCGCACAATAACAAGGTGCAAGAGCATGACGTATTTCATTGATGTTCCAACCATGTCGGACTTGGTGCATGACATTGGCGTAGCGCCGTTTATTGGCGAGCTTGCCGCTGCCCTGCGGGACGATTTCAAACGCTGGCAGGCGTTTGACAAGTCCGCGCGCGTTGCCAGCCACTCGGAAGTGGGCGTGATCGAACTGATGCCGGTCGCCGACAAAAGCCGCTATGCCTTCAAATACGTCAACGGCCACCCGGCCAATACTGCACGCAACCTGCATACGGTGATGGCCTT  >  NC_002947/4008992‑4009290

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