Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,919,491 A→G D262G (GAC→GGC)  queF → NADPH‑dependent 7‑cyano‑7‑deazaguanine reductase QueF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,919,4910AG95.4% 372.7 / 6.2 108D262G (GAC→GGC) queFNADPH‑dependent 7‑cyano‑7‑deazaguanine reductase QueF
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (52/51);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 108 total aligned reads displayed.

TCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTG  >  NZ_CP009273/2919365‑2919623
                                                                                                                              |                                                                                                                                    
tctcATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCg                                                                                                                            <  2:438287/137‑1 (MQ=255)
tctcATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCg                                                                                                                            >  1:438287/1‑137 (MQ=255)
tctcATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTg                                                                                                                          >  2:489271/1‑139 (MQ=255)
    aTTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTgg                                                                                                                         <  2:1101642/136‑1 (MQ=255)
    aTTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgcg                                                                                                                       <  1:1219142/138‑1 (MQ=255)
    aTTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgcg                                                                                                                       <  1:1279789/138‑1 (MQ=255)
     ttCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgcg                                                                                                                       >  2:613858/1‑137 (MQ=255)
     ttCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgcg                                                                                                                       >  1:575517/1‑137 (MQ=255)
          tcatcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTaa                                                                                                                 >  2:1046712/1‑138 (MQ=255)
          tcatcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTaa                                                                                                                 >  2:840961/1‑138 (MQ=255)
            atcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAg                                                                                                              <  2:1000897/139‑1 (MQ=255)
            atcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAg                                                                                                              <  1:762235/139‑1 (MQ=255)
                 aaCGAGTTCCACGAACAGCGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGAGGCGGCCCTGGCTTTAACCCGTGGCGCATCAATAGCGttt                                                                                                         >  2:343738/1‑139 (MQ=255)
                           aCGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATc                                                                                               <  1:1015695/139‑1 (MQ=255)
                           aCGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATc                                                                                               <  2:162407/139‑1 (MQ=255)
                           aCGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATc                                                                                               <  1:514499/139‑1 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgcg                                                                                                                       <  2:88263/107‑1 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgcg                                                                                                                       >  1:88263/1‑107 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgc                                                                                                                        >  1:903062/1‑106 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgc                                                                                                                        <  2:903062/106‑1 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTc                                                                                                    >  2:545311/1‑126 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTc                                                                                                    <  1:545311/126‑1 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCa                                                                                                 >  1:66048/1‑129 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCa                                                                                                 <  2:66048/129‑1 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACcaca                                                                                        >  2:1039887/1‑138 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACcac                                                                                         >  2:1105702/1‑137 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACcac                                                                                         >  1:759487/1‑137 (MQ=255)
                                   tGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACcac                                                                                         >  1:261951/1‑137 (MQ=255)
                                     cGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAg                                                                                      <  2:575517/138‑1 (MQ=255)
                                      gTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACa                                                                                                                                    >  2:742178/1‑91 (MQ=255)
                                      gTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACa                                                                                                                                    <  1:742178/91‑1 (MQ=255)
                                       tGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACt                                                                                   >  1:519785/1‑139 (MQ=255)
                                           aCGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGtt                                                                               <  1:219856/139‑1 (MQ=255)
                                           aCGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGtt                                                                               >  1:1194020/1‑139 (MQ=255)
                                             gCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgcg                                                                                                                       <  2:713123/97‑1 (MQ=255)
                                             gCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgcg                                                                                                                       >  1:713123/1‑97 (MQ=255)
                                             gCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACcac                                                                                         <  2:907633/127‑1 (MQ=255)
                                             gCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGACTGGGCTTTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACcac                                                                                         >  1:907633/1‑127 (MQ=255)
                                               aTCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTgg                                                                                                                         <  1:908577/93‑1 (MQ=255)
                                               aTCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTgg                                                                                                                         >  2:908577/1‑93 (MQ=255)
                                               aTCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCgg                                                                            >  1:907204/1‑138 (MQ=255)
                                                 cTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGgcg                                                                                                                       >  2:562371/1‑93 (MQ=255)
                                                 cTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGgcg                                                                                                                       <  1:562371/93‑1 (MQ=255)
                                                    tAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACca                                                                                          <  1:648856/119‑1 (MQ=255)
                                                    tAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACca                                                                                          >  2:648856/1‑119 (MQ=255)
                                                           cTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAAtttttt                                                               >  2:470921/1‑139 (MQ=255)
                                                           cTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAAtttttt                                                               >  1:878775/1‑139 (MQ=255)
                                                           cTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAAtttttt                                                               >  1:1190567/1‑139 (MQ=255)
                                                           cTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAAtttttt                                                               >  1:904321/1‑139 (MQ=255)
                                                             gTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTtct                                                             >  1:878129/1‑139 (MQ=255)
                                                              ttACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgcg                                                                                                                       <  2:727044/80‑1 (MQ=255)
                                                              ttACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGgcg                                                                                                                       >  1:727044/1‑80 (MQ=255)
                                                              ttACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTtct                                                             <  1:178998/138‑1 (MQ=255)
                                                              ttACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTtct                                                             >  2:857449/1‑138 (MQ=255)
                                                              ttACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTtct                                                             >  2:1103191/1‑138 (MQ=255)
                                                              ttACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTtc                                                              <  2:600039/137‑1 (MQ=255)
                                                              ttACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAAACTGGTTCGGCAATAAATTTTTtct                                                             >  1:44849/1‑138 (MQ=255)
                                                                        gccagccaGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTTTGGGCATTAACCCGCTGCGCAGTAATAACGATTTTCGCCCATTGCCCACAAGACCTGGTCGGCAATAAATTTTTTCTCAATTTTTCGt                                                  >  2:559178/1‑139 (MQ=255)
                                                                        gccagccaGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTTCgc                                                  >  2:1237359/1‑138 (MQ=255)
                                                                        gccagccaGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGt                                                  >  2:1173303/1‑139 (MQ=255)
                                                                        gccagccaGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGt                                                  >  2:1240756/1‑139 (MQ=255)
                                                                         ccagccagAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGc                                                    <  1:72835/136‑1 (MQ=255)
                                                                                        gAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAg                                   <  2:406527/138‑1 (MQ=255)
                                                                                        gAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAg                                   <  2:907204/138‑1 (MQ=255)
                                                                                         aGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGa                                  <  2:1194020/138‑1 (MQ=255)
                                                                                         aGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGa                                  <  1:798736/138‑1 (MQ=255)
                                                                                         aGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGa                                  <  2:283907/138‑1 (MQ=255)
                                                                                                 cGCACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCCTCGACCACAAGACTGGTTCGGCAATAAGTTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTtgtg                          <  1:500076/138‑1 (MQ=255)
                                                                                                   cACGTTATACCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGaaa                       <  2:594759/139‑1 (MQ=255)
                                                                                                        tataCCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGtca                  >  1:918938/1‑139 (MQ=255)
                                                                                                          taCCCGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGtcat                 <  2:759487/138‑1 (MQ=255)
                                                                                                            cccGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGtcatca               <  1:898795/138‑1 (MQ=255)
                                                                                                             ccGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcagg             <  1:1173303/139‑1 (MQ=255)
                                                                                                             ccGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcagg             <  2:1190567/139‑1 (MQ=255)
                                                                                                             ccGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcagg             <  1:1046712/139‑1 (MQ=255)
                                                                                                             ccGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcagg             <  1:1039887/139‑1 (MQ=255)
                                                                                                             ccGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcagg             <  2:904321/139‑1 (MQ=255)
                                                                                                             ccGTCGTGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcagg             <  2:878775/139‑1 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATc                                                                                               >  2:296245/1‑54 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATc                                                                                               <  1:296245/54‑1 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATc                                                                                               <  1:1272062/54‑1 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGAtt                                          >  1:1035209/1‑107 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGAtt                                          <  2:1035209/107‑1 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAg                                >  1:130319/1‑117 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAg                                <  2:130319/117‑1 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcaggcag          >  1:819316/1‑139 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcaggcag          >  2:924869/1‑139 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcaggcag          >  1:501503/1‑139 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcaggcag          >  2:349822/1‑139 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGTGGCGCAGGAATAACGATTTTGTCCCATCGACCACAAGACTGGGTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGGAGGTTGTGAAAGGGCATcaggcag          >  2:134072/1‑139 (MQ=255)
                                                                                                                tcgtGGCGGTCTGGGCATTAACCCGGGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCCATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATcaggcag          >  1:24323/1‑139 (MQ=255)
                                                                                                                  gtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAAttt                                                       >  1:88348/1‑92 (MQ=255)
                                                                                                                  gtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAAttt                                                       <  2:88348/92‑1 (MQ=255)
                                                                                                                  gtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAggg        <  1:857449/139‑1 (MQ=255)
                                                                                                                  gtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAggg        <  1:616815/139‑1 (MQ=255)
                                                                                                                  gtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAggg        <  2:1170871/139‑1 (MQ=255)
                                                                                                                  gtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAggg        >  1:115099/1‑139 (MQ=255)
                                                                                                                  gtGGCGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAggg        <  1:1105702/139‑1 (MQ=255)
                                                                                                                      cGGTCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTAt    <  2:819316/139‑1 (MQ=255)
                                                                                                                         tCTGGGCATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTg  >  2:257157/1‑138 (MQ=255)
                                                                                                                              |                                                                                                                                    
TCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTG  >  NZ_CP009273/2919365‑2919623

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: