Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP025268 1,936,791 (G)6→5 pseudogene (1097/1180 nt) purT → formate‑dependent phosphoribosylglycinamide formyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0252681,936,7860G.99.2% 1190.4 / 1.1 251pseudogene (1092/1180 nt)purTformate‑dependent phosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base . (99/149);  total (100/151)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

TGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGAT  >  NZ_CP025268/1936648‑1936914
                                                                                                                                          |                                                                                                                                
tGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTgg                                                                                                                                 >  1:547874/1‑140 (MQ=255)
tGGTCCTGCAGCTTCTGCCGTTATTATGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATGCGTTTATTTGGTAAGCCGGAAATTGATGGAAGCCGTCGTCTgg                                                                                                                                 >  2:4675744/1‑140 (MQ=255)
 ggTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggg                                                                                                                                <  1:3734801/140‑1 (MQ=255)
  gTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTgggg                                                                                                                               <  2:1741201/140‑1 (MQ=255)
   tCCTGCAGCTTCTGCCGTTCTTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggggg                                                                                                                              <  1:524794/140‑1 (MQ=255)
   tCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCATTCAGATTGTCACGTTTAATAATGTGCAGGATGCAGTAGTCGTAGATTGGCAGCTTCGTATATTTGGTAAGCTGGAAATTGAGGGCATCCGTCGTCTggggg                                                                                                                              <  2:697211/140‑1 (MQ=255)
   tCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggggg                                                                                                                              >  1:1485078/1‑140 (MQ=255)
   tCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggggg                                                                                                                              <  2:458334/140‑1 (MQ=255)
   tCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggggg                                                                                                                              <  2:3532476/140‑1 (MQ=255)
   tCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggggg                                                                                                                              <  1:4112792/140‑1 (MQ=255)
   tCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggggg                                                                                                                              <  2:2331107/140‑1 (MQ=255)
    ccTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTgggggt                                                                                                                             >  1:3343172/1‑139 (MQ=255)
     cTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTgggggtg                                                                                                                            <  1:1227192/140‑3 (MQ=255)
        gggCTTCTGCCGTTATTCTGCCACAACTGACCGGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCa                                                                                                                        <  2:2877665/138‑1 (MQ=255)
        cAGCTTCTGCTGTTATTCTGCCACACCTGCCCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCa                                                                                                                        <  2:524597/140‑1 (MQ=255)
        cAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCa                                                                                                                        <  2:3753902/140‑1 (MQ=255)
        cAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCa                                                                                                                        <  2:3137876/140‑1 (MQ=255)
        cAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCa                                                                                                                        <  1:2175802/140‑1 (MQ=255)
        cAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCa                                                                                                                        <  2:5109586/140‑1 (MQ=255)
        cAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCa                                                                                                                        >  1:1708383/1‑140 (MQ=255)
        cAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTCGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTgggggtgtca                                                                                                                         >  2:4814313/1‑135 (MQ=255)
          gCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACt                                                                                                                      <  1:4023352/140‑1 (MQ=255)
          gCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACt                                                                                                                      <  2:1579684/140‑1 (MQ=255)
          gCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACt                                                                                                                      <  2:2832817/140‑1 (MQ=255)
          gCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACt                                                                                                                      <  1:2371402/140‑1 (MQ=255)
          gCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACt                                                                                                                      <  1:4970567/140‑1 (MQ=255)
          ccTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACt                                                                                                                      <  1:3271210/139‑1 (MQ=255)
            ttCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTgg                                                                                                                    <  1:679591/140‑1 (MQ=255)
            ttCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTgg                                                                                                                    <  2:1911082/140‑1 (MQ=255)
            ttCTGCCGTTATTCTGCCACAACTGACCAGGCAGAATGTCACGTTTGCTAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTgg                                                                                                                    <  2:4179247/140‑1 (MQ=255)
             tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGc                                                                                                                   <  1:312219/140‑1 (MQ=255)
             tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGc                                                                                                                   >  1:3897067/1‑140 (MQ=255)
             tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGc                                                                                                                   >  1:644587/1‑140 (MQ=255)
             tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGc                                                                                                                   >  1:3673176/1‑140 (MQ=255)
             tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTCGGGGTGGCACTGGc                                                                                                                   >  2:2512377/1‑140 (MQ=255)
              cTGCCGTTATTCTGCAACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCt                                                                                                                  <  1:2702031/140‑1 (MQ=255)
                gCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTAc                                                                                                                <  1:3084771/140‑1 (MQ=255)
                    ttATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCa                                                                                                            <  2:130502/140‑1 (MQ=255)
                    ttATGCTGCCACAACTGACCAGTCAGAATGTCCCGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGGCGTCGTCTGGGGGTGGCACTGGCTACTGCa                                                                                                            <  2:3566530/140‑1 (MQ=255)
                      aTTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCaga                                                                                                          <  1:3859336/140‑1 (MQ=255)
                      aTTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCaga                                                                                                          <  1:1715037/140‑1 (MQ=255)
                           gCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAgtg                                                                                                     <  1:2350652/140‑1 (MQ=255)
                            ccACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAgtgt                                                                                                    >  1:847303/1‑140 (MQ=255)
                                 aCTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGttgttg                                                                                               >  2:1981175/1‑140 (MQ=255)
                                  cTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGa                                                                                              <  2:3822794/140‑1 (MQ=255)
                                  cTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGa                                                                                              <  1:1776137/140‑1 (MQ=255)
                                    gACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACg                                                                                            >  2:846300/1‑140 (MQ=255)
                                      ccAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGcc                                                                                          >  1:5134449/1‑140 (MQ=255)
                                        aGTCAGAATGTCACGCTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGAGTCGTTTATTTGGTAAGCCGGACATTGATGGCAGCCGTCGTCTGGCGGCGGCACTGGCTACTGCAGAGGGTGATGTTGACGCCAt                                                                                        >  2:3488035/1‑140 (MQ=255)
                                          tccgagtGTCACTTGGGATACTGGGTAGATTGCCGTAGTGCCAGCGTTGCCGATTCGTTTATTTGGTAAGCCGGCAATTGATGGCGGCCGTCGTCTGGGGGTGGCACGGGCTACTGGAGGGAGTGTTGTTGACGCCATTg                                                                                      <  2:1787545/134‑1 (MQ=255)
                                           cAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGa                                                                                     <  2:216069/140‑1 (MQ=255)
                                            aGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGCCGCCATTGaa                                                                                    <  1:2224455/140‑1 (MQ=255)
                                                    aCGTTTGATAATGTGCAGAATGGCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGaa                                                                            >  2:2703193/1‑140 (MQ=255)
                                                    aCGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGaa                                                                            >  1:4204322/1‑140 (MQ=255)
                                                    aCGTTTGATAATGTGCAGAATGCCGTAAGCGCAGATTTGCAGATTCGTATATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGCGTGGGGGTGGCACTGGCTACTAGAGAGAGTGTTGTTGACGCCATTCTGCGCGCGaa                                                                            >  2:3804343/1‑140 (MQ=255)
                                                       tttGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGAACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                                         <  1:1401859/140‑1 (MQ=255)
                                                       tttGATAATGTGCAGAATGACGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                                         <  2:2785582/140‑1 (MQ=255)
                                                          gATAATGTGCAGAATGCCGTACGCGCAGATTTGCAGATTCGTTTAGTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGCGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCAcgc                                                                      >  2:2989005/1‑140 (MQ=255)
                                                           aTAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGACGCAcgcc                                                                     >  1:2074287/1‑140 (MQ=255)
                                                            tAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCAcgccg                                                                    <  1:3432315/140‑1 (MQ=255)
                                                            tAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCAcgccg                                                                    <  1:360826/140‑1 (MQ=255)
                                                                gtgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCgg                                                                <  2:113625/140‑1 (MQ=255)
                                                                 tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                               >  1:5059395/1‑140 (MQ=255)
                                                                 tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                               >  2:2118882/1‑140 (MQ=255)
                                                                 tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                               <  1:2932273/140‑1 (MQ=255)
                                                                 tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                               >  2:3232894/1‑140 (MQ=255)
                                                                 tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                               <  2:1592635/140‑1 (MQ=255)
                                                                 tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                               >  1:1340812/1‑140 (MQ=255)
                                                                 tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                               >  1:1149056/1‑140 (MQ=255)
                                                                 tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                               <  1:4552434/140‑1 (MQ=255)
                                                                 tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGACACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCCAAGCACGCCGCCGGa                                                               <  2:1942891/140‑1 (MQ=255)
                                                                   cAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACa                                                             <  1:4681723/140‑1 (MQ=255)
                                                                   cAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACa                                                             <  1:2586899/140‑1 (MQ=255)
                                                                    aGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggggg                                                                                                                              <  2:4121606/75‑1 (MQ=255)
                                                                      aaTGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGt                                                          <  2:742589/140‑1 (MQ=255)
                                                                      aaTGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGt                                                          >  1:4365068/1‑140 (MQ=255)
                                                                        tGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTaa                                                        <  1:1899401/140‑1 (MQ=255)
                                                                        tGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTaa                                                        >  1:3240862/1‑140 (MQ=255)
                                                                        tGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTaa                                                        <  1:1990222/140‑1 (MQ=255)
                                                                        tGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTaa                                                        >  1:2412484/1‑140 (MQ=255)
                                                                        tGCCGGACGCGCAGGTATGCAGATTCGTGTATATGGTAAGGCGGAAATTGATGGCAGCTGTCGCCTGGTGGTGCCACTGCCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGCAGCACGCCGCCGGACAGGTat                                                        <  2:4306639/140‑2 (MQ=255)
                                                                          ccGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGGGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACggggcaaa                                                      >  2:101865/1‑135 (MQ=255)
                                                                             tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTa                                                   <  1:1521311/140‑1 (MQ=255)
                                                                             tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTa                                                   <  1:3188596/140‑1 (MQ=255)
                                                                             tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTa                                                   <  2:3412039/140‑1 (MQ=255)
                                                                             tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTa                                                   <  2:2011385/140‑1 (MQ=255)
                                                                             tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTa                                                   <  1:355199/140‑1 (MQ=255)
                                                                             tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTa                                                   >  1:723262/1‑140 (MQ=255)
                                                                               ggCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACa                                                 <  1:3926405/140‑1 (MQ=255)
                                                                                                tCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGCGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGCTAAACCCGGCaaa                                >  2:3060704/1‑140 (MQ=255)
                                                                                                            tAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACaa                    >  2:2622417/1‑140 (MQ=255)
                                                                                                            tAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACaa                    >  2:4358945/1‑140 (MQ=255)
                                                                                                            tAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACaa                    >  1:1668124/1‑140 (MQ=255)
                                                                                                            tAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACaa                    >  2:620161/1‑140 (MQ=255)
                                                                                                             aaGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACaaa                   <  1:3654705/140‑1 (MQ=255)
                                                                                                                   gAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGATGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGcc             <  1:2843345/140‑1 (MQ=255)
                                                                                                                   gAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGcc             <  1:3450506/140‑1 (MQ=255)
                                                                                                                   gAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGcc             <  2:1402869/140‑1 (MQ=255)
                                                                                                                   gAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGcc             <  1:1037764/140‑1 (MQ=255)
                                                                                                                          aTGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGCAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAACCCGCTACAAAAATGCCCGATCCt      >  2:3994775/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                
TGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGAT  >  NZ_CP025268/1936648‑1936914

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