Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,226,250 C→T T808M (ACG→ATG)  metH → cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,226,2500CT98.6% 250.1 / ‑5.9 71T808M (ACG→ATG) metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (37/33);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCACGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAGCAGAACTACAGCGGCCC  >  NC_000913/4226114‑4226385
                                                                                                                                        |                                                                                                                                       
aaaaaTATCGTTGGTGTGGTGCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGcc                                                                                                                                      <  2:1047999/140‑1 (MQ=255)
        cGTTGGTGTGGTGCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTgg                                                                                                                              >  1:372323/1‑140 (MQ=255)
         gTTGGTGTGGTGCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGa                                                                                                                             >  2:1113084/1‑140 (MQ=255)
           tggtgtggtgCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACg                                                                                                                           <  2:372323/140‑1 (MQ=255)
            ggtgtggtgCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACga                                                                                                                          <  1:247301/140‑1 (MQ=255)
             gtgtggtgCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGc                                                                                                                                       >  1:326136/1‑126 (MQ=255)
             gtgtggtgCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACgag                                                                                                                         <  2:365868/140‑1 (MQ=255)
             gtgtggtgCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACgag                                                                                                                         >  2:1012462/1‑140 (MQ=255)
             gtgtggtgCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACgag                                                                                                                         >  1:294272/1‑140 (MQ=255)
                     ctgcAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTaa                                                                                                                 >  1:778469/1‑140 (MQ=255)
                        cAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACg                                                                                                               >  2:373775/1‑139 (MQ=255)
                        cAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACg                                                                                                               <  1:373775/139‑1 (MQ=255)
                        cAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGt                                                                                                              >  2:481168/1‑140 (MQ=255)
                         aaTGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTg                                                                                                             <  1:481168/140‑1 (MQ=255)
                            gTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCg                                                                                                          >  2:161561/1‑140 (MQ=255)
                             tAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAAc                                                                                                                >  1:220796/1‑133 (MQ=255)
                             tAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAAc                                                                                                                <  2:220796/133‑1 (MQ=255)
                             tAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGc                                                                                                           <  2:608119/138‑1 (MQ=255)
                             tAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGc                                                                                                           >  1:608119/1‑138 (MQ=255)
                             tAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGa                                                                                                         >  2:575222/1‑140 (MQ=255)
                             tAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGa                                                                                                         <  2:659569/140‑1 (MQ=255)
                              aacaacTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGaa                                                                                                        <  2:631351/140‑1 (MQ=255)
                               acaacTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGaaa                                                                                                       <  1:724402/140‑1 (MQ=255)
                                  acTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGTCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAagag                                                                                                    >  1:51232/1‑140 (MQ=255)
                                  acTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAagag                                                                                                    >  2:939479/1‑140 (MQ=255)
                                     aCGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATg                                                                                                 <  1:1113084/140‑1 (MQ=255)
                                     aCGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATg                                                                                                 >  2:427542/1‑140 (MQ=255)
                                     aCGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATg                                                                                                 >  2:734967/1‑140 (MQ=255)
                                     aCGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCACGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATg                                                                                                 >  2:511142/1‑140 (MQ=255)
                                           ttGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGt                                                                                           <  1:1012462/140‑1 (MQ=255)
                                             gTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCa                                                                                         >  1:951942/1‑140 (MQ=255)
                                                gATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAggg                                                                                      <  1:427542/140‑1 (MQ=255)
                                                     cGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCa                                                                                 <  1:575222/140‑1 (MQ=255)
                                                     cGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCa                                                                                 >  1:297995/1‑140 (MQ=255)
                                                      ggCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCAc                                                                                <  1:939479/140‑1 (MQ=255)
                                                        cGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTCTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAggg                                                                                      >  2:377235/1‑132 (MQ=255)
                                                        cGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTCTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAggg                                                                                      <  1:377235/132‑1 (MQ=255)
                                                          ttATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTAtt                                                                            <  1:586239/140‑1 (MQ=255)
                                                          ttATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTAtt                                                                            <  2:135583/140‑1 (MQ=255)
                                                           tATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTc                                                                           <  1:604471/140‑1 (MQ=255)
                                                                tGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTAtt                                                                            <  1:501518/134‑1 (MQ=255)
                                                                tGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTAtt                                                                            >  2:501518/1‑134 (MQ=255)
                                                                    tGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTGAGGGCTTCACTATTCCGTTACTGa                                                                  >  1:1102972/1‑140 (MQ=255)
                                                                    tGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGa                                                                  >  2:403545/1‑140 (MQ=255)
                                                                    tGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGa                                                                  >  1:141342/1‑140 (MQ=255)
                                                                     gCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTTACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGAt                                                                 >  1:6020/1‑140 (MQ=255)
                                                                         aaaaaaTTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTggc                                                             <  1:161561/140‑1 (MQ=255)
                                                                                 ctcCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGcgac                                                     <  2:1120928/140‑1 (MQ=255)
                                                                                        ccGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTTACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCa                                              <  2:6020/140‑1 (MQ=255)
                                                                                           cTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTGAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCaaaa                                           <  2:1102972/140‑1 (MQ=255)
                                                                                                  agTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGcaca                                    <  2:778469/140‑1 (MQ=255)
                                                                                                   gTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCacac                                   >  2:728082/1‑140 (MQ=255)
                                                                                                   gTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCacac                                   >  2:1045142/1‑140 (MQ=255)
                                                                                                     gAATGCTGATCTGATTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACAcgg                                 >  1:643186/1‑140 (MQ=255)
                                                                                                             atctgatTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGaaa                         >  2:1071700/1‑140 (MQ=255)
                                                                                                                tgatTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATc                      >  1:808547/1‑140 (MQ=255)
                                                                                                                tgatTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATc                      >  1:164473/1‑140 (MQ=255)
                                                                                                                 gatTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCg                     <  2:761735/140‑1 (MQ=255)
                                                                                                                  atTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGa                    >  2:929751/1‑140 (MQ=255)
                                                                                                                  atTGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGa                    >  2:89722/1‑140 (MQ=255)
                                                                                                                   ttGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAg                   <  1:1045142/140‑1 (MQ=255)
                                                                                                                   ttGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAg                   <  2:141342/140‑1 (MQ=255)
                                                                                                                   ttGGCCTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAg                   <  2:919806/140‑1 (MQ=255)
                                                                                                                       ccTTTCGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAGCAGa               >  1:220891/1‑140 (MQ=255)
                                                                                                                         tttCGGGGCTTATCATGCCGTCGCTAGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAGCAGAAc             <  2:246908/140‑1 (MQ=255)
                                                                                                                            cGGGGCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAGCAGAACTAc          >  2:167946/1‑140 (MQ=255)
                                                                                                                                gCTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAGCAGAACTACAGCg      <  1:734967/140‑1 (MQ=255)
                                                                                                                                 cTTATCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAGCAGAACTACAGCgg     >  2:378989/1‑140 (MQ=255)
                                                                                                                                    aTCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAGCAGAACTACAGCGGccc  >  1:513770/1‑140 (MQ=255)
                                                                                                                                    aTCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAGCAGAACTACAGCGGccc  >  1:1089077/1‑140 (MQ=255)
                                                                                                                                    aTCATGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAGCAGAACTACAGCGGccc  <  2:297995/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGAAAAAATTCTCCGTACCGCTAAAGAAGTGAATGCTGATCTGATTGGCCTTTCGGGGCTTATCACGCCGTCGCTGGACGAGATGGTTAACGTGGCGAAAGAGATGGAGCGTCAGGGCTTCACTATTCCGTTACTGATTGGCGGCGCGACGACCTCAAAAGCGCACACGGCGGTGAAAATCGAGCAGAACTACAGCGGCCC  >  NC_000913/4226114‑4226385

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