Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,057,920 G→A 6.5% T10M (ACG→ATG)  cobT ← nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,057,9200GA 6.5% 176.6 / 2.9 62T10M (ACG→ATG) cobTnicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base A (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

CCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGA  >  CP009273/2057774‑2058066
                                                                                                                                                  |                                                                                                                                                  
ccACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTa                                                                                                                                                  >  2:315751/1‑149 (MQ=255)
  acaTGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTAtt                                                                                                                                                >  2:511219/1‑149 (MQ=255)
          tacgCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCCCCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAAt                                                                                                                                        <  2:656684/146‑1 (MQ=255)
          tATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGCTGGCTGGGATCGTCTTCAGTAAAt                                                                                                                                        >  2:315086/1‑149 (MQ=255)
          tATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAAt                                                                                                                                        <  2:582296/149‑1 (MQ=255)
               cATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGGAAACCggaa                                                                                                                                   >  1:402029/1‑147 (MQ=255)
                       cccGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTg                                                                                                                           >  1:79594/1‑149 (MQ=255)
                          ggCATTCCTGCCAGTGGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGGTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCAt                                                                                                                        <  1:173327/149‑1 (MQ=255)
                              ttCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCACGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCa                                                                                                                    <  1:550674/149‑1 (MQ=255)
                              ttCCCGCCAGTTGTATGGCAAGCACCTCCTGCAGTCCCCGGCTCCCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTAGACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTGAGTCCATCGGCAACTATTTGCATcaca                                                                                                                    <  1:489548/149‑1 (MQ=255)
                                ccTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATAGGCAAGTATTTGCATCTCACg                                                                                                                  >  2:686012/1‑149 (MQ=255)
                                      aGTTGTTTGGCAAGCACCCCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGGCAATATGCCGTCGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCtt                                                                                                            <  2:435801/149‑1 (MQ=255)
                                       gTTGTATGCCAAGCACCTCCGGCTTTCCCCGGCTACCAACAGGTTTCAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCttt                                                                                                           <  1:206083/149‑1 (MQ=255)
                                         tGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTCACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTtat                                                                                                         <  1:530074/149‑1 (MQ=255)
                                         tGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTtat                                                                                                         <  2:602896/149‑1 (MQ=255)
                                              ggtaaGCCCCTCCACTTTCCCCTGGCTACCAACAGTTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATCGCAg                                                                                                    <  1:446481/146‑1 (MQ=255)
                                               gCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGGTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGGCATAGCGGTAGAACCGATGCCTGGGACCGTATGCAGTAACTCGGCACGGATTTCCATCTCACGTCTCTTTATAGCCGt                                                                                                   >  1:1262/1‑149 (MQ=255)
                                                cAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCATATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTg                                                                                                  >  1:451880/1‑149 (MQ=255)
                                                cAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCATATTCAGTAAATCGGCAAGTATTTGCATCTCACGGCTCTTTATAGCCGTg                                                                                                  >  1:451582/1‑149 (MQ=255)
                                                 aaGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGc                                                                                                 >  1:345722/1‑149 (MQ=255)
                                                   gCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCca                                                                                               >  1:423683/1‑149 (MQ=255)
                                                     aCCTCCTGCTTTCCCATGCTACCCACAGGTGTGAGTAAACCGGCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCcagc                                                                                             <  1:14246/149‑1 (MQ=255)
                                                     aCCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGGATTCAGTAAATCGGCACGTATTTGCATCTCAAGTCCCTTGATCGCAGTGCcaga                                                                                             >  1:749413/1‑148 (MQ=255)
                                                      ccTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCcagca                                                                                            <  2:759701/149‑1 (MQ=255)
                                                        tCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAgg                                                                                          <  2:79594/149‑1 (MQ=255)
                                                        tCCAGCTTTCCCAGGCTACAAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCaggagg                                                                                          >  1:589398/1‑149 (MQ=255)
                                                         ccAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGa                                                                                         >  1:550332/1‑149 (MQ=255)
                                                           aGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAg                                                                                       >  2:401864/1‑149 (MQ=255)
                                                                   ctggCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTc                                                                               >  1:581500/3‑149 (MQ=255)
                                                                   cAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTc                                                                               >  2:380583/1‑149 (MQ=255)
                                                                   cAGCCTACCAAAAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTc                                                                               <  2:427022/149‑1 (MQ=255)
                                                                        tACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTCAATCGGCAAGTATTTGCATCTCACGGCTCTTTATAGCGGTGCCAGCAGGAAGGCCAGTTCACCaa                                                                          >  2:775088/1‑149 (MQ=255)
                                                                         aCCCACAGGTTTGAGGAACCCGGCAATATGCCGTCGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAg                                                                         <  2:60085/149‑1 (MQ=255)
                                                                         aCCAACAGGTTTGAGTAACCCGTCAGTATGACGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAg                                                                         >  2:448804/1‑149 (MQ=255)
                                                                          ccAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCTATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGc                                                                        >  1:642360/1‑149 (MQ=255)
                                                                          ccAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCACGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGc                                                                        <  1:572558/149‑1 (MQ=255)
                                                                          ccAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGc                                                                        <  2:231405/149‑1 (MQ=255)
                                                                          ccAACAGGTTTGAGTAACCCGTCAATATGCCGGTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGc                                                                        <  2:287513/149‑1 (MQ=255)
                                                                             aCAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCCGTTCACCAAGctct                                                                     >  1:289842/1‑149 (MQ=255)
                                                                             aCAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCCCCAAGctct                                                                     >  2:261065/1‑149 (MQ=255)
                                                                             aCAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATGCGGTAAATCGGCAAGTATTTGCATCTCCCGGCTCTGTATGGCAGGGCCAGCAGGGAGACCAGGTCACCAGGctcg                                                                     >  1:679676/1‑148 (MQ=255)
                                                                                 gTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGCAGACCAGTTCACCAAGCTCTATgg                                                                 >  1:72104/1‑149 (MQ=255)
                                                                                 gTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATgg                                                                 <  1:606974/149‑1 (MQ=255)
                                                                                  tttGAGTAACCCGTCAATATGCCGTTTGGAACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCACGCTCTATGgc                                                                >  1:658442/1‑149 (MQ=255)
                                                                                       gTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGcctgccc                                                           >  1:413202/1‑147 (MQ=255)
                                                                                              gTCAATATGCCGTTGTGCAGGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGt                                                    >  2:630261/1‑149 (MQ=255)
                                                                                                      gCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCATATTCAGTAAATCGTCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCACGCTCTATGGCTGCTCCCAGCGTATCGCCCg                                            <  2:451880/149‑1 (MQ=255)
                                                                                                      gCCATTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCATATTCAGTAAATCGCCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCACGCTCTATGGCTGCTCCCAGCGTATCGCCCg                                            <  2:451582/149‑1 (MQ=255)
                                                                                                       ccGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAAGAAGACCAGGTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGt                                           >  2:455125/1‑149 (MQ=255)
                                                                                                           tgtgCACGCGACATAGCGGTAGAATCGCTGGCTGGCACGGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTgc                                       <  2:163427/149‑1 (MQ=255)
                                                                                                           tgtgCACCCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCCATTCAGTAGTGTTGATTGCTGCTCCCAGCGTATCGCCCGTTTgc                                       <  2:425077/149‑1 (MQ=255)
                                                                                                           tgggCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTgc                                       <  2:402029/146‑1 (MQ=255)
                                                                                                            gtgCCCGCGACATAGCGGGAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTCTTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTgcc                                      <  1:761355/149‑1 (MQ=255)
                                                                                                            gtgCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGGTTCGCCCGTTTgcc                                      >  1:255108/1‑149 (MQ=255)
                                                                                                              gCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTgcccg                                    >  1:692391/1‑149 (MQ=255)
                                                                                                                    taCATAGCGGTAGAAGCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGc                              <  1:544305/148‑1 (MQ=255)
                                                                                                                    gACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGc                              <  2:423683/149‑1 (MQ=255)
                                                                                                                        tAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTAc                          >  2:271197/1‑149 (MQ=255)
                                                                                                                                         gCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCg         >  2:191097/1‑149 (MQ=255)
                                                                                                                                           tGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCCATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAg       <  2:642360/149‑1 (MQ=255)
                                                                                                                                             ggATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGa     <  2:613110/149‑1 (MQ=255)
                                                                                                                                                tCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGCACGAGGATGa  >  1:74637/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGCAGTGCCAGCAGGAAGACCAGTTCACCAAGCTCTATGGCTGCTCCCAGCGTATCGCCCGTTTGCCCGCCCAGCGTACGTTTGAGAAGCTGACCGAGGATGA  >  CP009273/2057774‑2058066

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