Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,040,0750CG41.1% 78.0 / 25.3 59L200V (CTG→GTG) yedYmembrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base C (4/29);  new base G (23/0);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 5.54e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGG  >  NC_000913/2039929‑2040224
                                                                                                                                                  |                                                                                                                                                      
cGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGaa                                                                                                                                                    <  2:475727/151‑1 (MQ=255)
 gATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGaaa                                                                                                                                                   <  2:352855/151‑1 (MQ=255)
    ggtggtGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAatat                                                                                                                                                <  1:807133/151‑1 (MQ=255)
      tggtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATcc                                                                                                                                              >  2:523906/1‑151 (MQ=255)
       ggtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGgaaaaacc                                                                                                                                              >  1:739092/1‑145 (MQ=255)
        gtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCtt                                                                                                                                            <  2:526965/151‑1 (MQ=255)
                    ggTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGa                                                                                                                                <  1:639227/150‑1 (MQ=255)
                    ggTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAAAATCCTTATGTCGAAgg                                                                                                                                 >  2:427439/1‑149 (MQ=255)
                          ccGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGt                                                                                                                          <  2:437140/150‑1 (MQ=255)
                               gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTcgac                                                                                                                    <  1:686509/151‑1 (MQ=255)
                                cacaAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGGAAAATACTTATTTCGGAGGGGTGCGtcttgccg                                                                                                                   >  1:184503/1‑146 (MQ=255)
                                   aaaTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGCGGAACATCCTTTTCAAGAAGGAGTGCGcctc                                                                                                                       >  1:384074/1‑144 (MQ=255)
                                    aaTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTGGAAAATTCTTATTTGGAAGGGTTGCTTCTCAACGaa                                                                                                                 >  1:875468/1‑149 (MQ=255)
                                       tGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCaa                                                                                                             <  2:243457/150‑1 (MQ=255)
                                             cgcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCCTATTTCGAAGGAATGGGTCTCGACGAAGCAAAGCATcc                                                                                                      >  2:594380/1‑151 (MQ=255)
                                             cgcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATACTTATGTCGAAGGGTTGGGGCTCGACGGAGCAAAGAATcc                                                                                                      >  2:779273/1‑151 (MQ=255)
                                             cgcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTTGAAGGATTGCGTCTCGACGAAGCAATGCATcc                                                                                                      >  2:765485/1‑151 (MQ=255)
                                             cgcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATcc                                                                                                      <  1:743986/151‑1 (MQ=255)
                                                 tGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGc                                                                                                    <  2:955447/149‑1 (MQ=255)
                                                  gCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTc                                                                                                  <  1:55811/150‑1 (MQ=255)
                                                    cGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGCAAAATTCTTATTTCGGAAGGATGCGGCTCGACGGAgaa                                                                                                              >  2:829941/1‑134 (MQ=255)
                                                    cGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATTCTTATGTCGAAAGATTGGGGCTCGACGAAGCAATGCATCCCCTc                                                                                                  >  2:154784/1‑148 (MQ=255)
                                                     gAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCacac                                                                                              <  2:685351/151‑1 (MQ=255)
                                                     gAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCacac                                                                                              <  1:405113/151‑1 (MQ=255)
                                                       aCCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGTAATATCCTTATGTTGGCGGGTTGCGTCTCCACGGAGCTATgt                                                                                                          >  1:126645/1‑136 (MQ=255)
                                                           aCCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACtgatga                                                                                        <  2:106446/151‑1 (MQ=255)
                                                           aCCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACtgatg                                                                                         <  2:816703/150‑1 (MQ=255)
                                                            ccAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGAGGGGTGCAAAATACTTATTTCCAAGGGATGCGTCTCGACGGGGGAAAGCAACCGCCCACACtcatga                                                                                        >  1:169906/1‑150 (MQ=255)
                                                                    gCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGgt                                                                                <  1:937032/150‑1 (MQ=255)
                                                                        aGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTCGAAGGATTGGGTCTCGGCGAAGCAAAGCAATCGCCCCCACCGATGACCGGAgggg                                                                               >  1:341499/1‑145 (MQ=255)
                                                                           aTGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCTTTTGTCGGAGGATTTGGTCTCGACGGAGCAAAGGATTCGCCCACACTGATTACCGTAgggggtt                                                                            >  2:426808/1‑142 (MQ=255)
                                                                           aTGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATgg                                                                        <  2:151792/151‑1 (MQ=255)
                                                                           aTGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATg                                                                         <  2:387003/150‑1 (MQ=255)
                                                                               cGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCaa                                                                     <  2:73739/150‑1 (MQ=255)
                                                                                gCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTCGAAGGATTTCGTCTCGACGGAGCAAAGCATCCCCCCCCACTGATGACCGGAGGGGGTTTTGGCAAgg                                                                   >  2:1030655/1‑151 (MQ=255)
                                                                                gCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAgg                                                                   <  1:681793/151‑1 (MQ=255)
                                                                                        aaCAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATTCCTATGTCGAAGGATTGCGTCTCGACGAAGCAAAGCATCCGCCCCCACTGATGAACCGAGGGGTTTATGGCAAGGCGGTAcc                                                            >  2:839861/1‑150 (MQ=255)
                                                                                             tttATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTTGGGTGGAAAATCATTATGTCCAAGGAATGGGGCTTGGCGAAGCAATGCATCCGCTCACACTGATGATAGCAGGGGTTTTTTGGAAGGGGGCACCGcc                                                         >  1:338921/1‑148 (MQ=255)
                                                                                             tttATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGGATCCGCTCACACTGATGACCGTAGGGGGTTATGGCCAGGCGTTACcccccc                                                       >  2:274031/1‑148 (MQ=255)
                                                                                              ttATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGAAAAAAACTTTTTTTGGAGGGTTGCGGTTGGAGGAAGGAAATCATCCCCTCACACTGATGGACGTAGGGGGTTTTGGGAAGGGGGTAACCCcac                                                       >  2:629123/1‑149 (MQ=255)
                                                                                                  gCACCGGAACAGATGCCGGCCCAGCAGCACCGTTTTATCGGCGGTGGGTGAAAATATCCTTATGTGGAAGGAGTGCGTCGCGACGAAGCAATGCACCCGCTCCCACTGATACCCGTGGGTGTTTGTGGCAAGGCGTTACCGCCACaaaa                                                   <  1:63117/149‑1 (MQ=255)
                                                                                                   cACCGGAACAGATGCCAGGCCAGCAGGCCCGTTTTATCGGGGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCTTCTCGACGAAGAAATGCACCCGCTCACACTGAGGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATg                                                 <  1:511414/150‑1 (MQ=255)
                                                                                                       ggAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGCAAAATCCCTATGTTCAAGGGTTGCGTCTCGACGAAGGAAAGCATCCGCTCACACTGCTGACCGGAGGGGTTTTTGGGAAGGAGGTACCGCCACcaaaaag                                                >  1:216845/1‑144 (MQ=255)
                                                                                                         agcAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGGGGAATAACCATTTGGCGAGGGATTGCGGCTCGACAAAGCAAAGCATCCGCGCACACAGATGACCGGGGGGGGTTATGGGGAGGCGGTACCGCCACCAAATggag                                              >  1:976359/3‑145 (MQ=255)
                                                                                                             gATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTTTTTCGAAGGATTGCGTCTCGACGGAGAAAAGCATCCCCTCACACTGATGACCGTAGGGGGTTTTGGCAAGGCGGTACCGCCCCAAAAAGGCGcggcggg                                       >  1:664381/1‑149 (MQ=255)
                                                                                                              aTGCCAGGCCAGCAGGCCCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGc                                     <  2:623782/151‑1 (MQ=255)
                                                                                                                   aggccagCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATTTCGAAGGATTGTGTCCCCACGGAGCAAAGCCTTCCCCCACACCTATGACCGGAGGGGTTTTTGGCAAAGGGTTTACGCCACCAAAAGGGGCGCCGGTGCGGCt                                 >  2:979277/1‑150 (MQ=255)
                                                                                                                   agaccagAAGGCCCCCTTTATCGGGGGGGGTTTAAAATTTCTTTATGTCGAAGGATTGCGTCTCGCCGAAGAAATGTATCCGCTCACACTGATGTCCGTAGGTTTTTATGGAAAGGGGTAACCGCCACAAAATGGCGCGCCGGTGCGACt                                 <  2:636089/147‑1 (MQ=255)
                                                                                                                    ggccagGAGGGCCGCTTTGTCGGGGGTGGGCGGAAATATCGTTAGGGCGATGGATTGCGCCTGGACGAAGAAATGCACCCGCTCCCGCTGGTGACGGGAGGTGTTTATGGCAGGGCGTGACCGCCACAAAATGGCGCGCCGGTGCGAca                                 <  1:328577/149‑2 (MQ=255)
                                                                                                                    ggccagCAGGCCCGTTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGa                               <  2:718890/151‑1 (MQ=255)
                                                                                                                      ccagCAGGACCGTTTTATCGGGGGTGGGCTGAAATATCTTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTGGGTGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGAtt                             <  1:820398/151‑1 (MQ=255)
                                                                                                                         ggagGACCGTTTTATGGGCGGTGGGTGAAAATATCCTTATGTCGAGGGATTGGGTCTCGAGGAAGCAATGCACCCGCTCACCCGGATGCCCGGGGGTTTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGAtt                             <  1:856615/146‑1 (MQ=255)
                                                                                                                          cagGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATGTCCAAGGGTTGCGGCTCGACGAAGCAATGCATCCCCTCACACCGATTAACGTAGGTGGTTATGGGAAGGCGTTAACCCCCCAAAAAGGCGCGGCGGGGCGACCGATtgt                           >  2:230349/1‑149 (MQ=255)
                                                                                                                            ggACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAGGGATTGCGTCCCGACGAAGAAATGCATCCGCCCACACGGATGACCGTGGGTTTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATtgtg                          <  1:668495/148‑1 (MQ=255)
                                                                                                                                    ttAGCGGCGGGGGGCTAAAATATCCTTATGTCGAAGGATTGCGCCTCGACGAGACAATGCATCCGCTCACACTGATGGCCGTAGGTGTTTATGGAAAGTCTTTGCCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAat                <  1:238025/150‑1 (MQ=255)
                                                                                                                                        cggcggTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGGGCGACTGATTGTGCCGTGGAAATATgg            >  1:948182/1‑150 (MQ=255)
                                                                                                                                         ggcggTGGGCTTAAATATCCTTATGGCGAAGGATTGCTTCTGGACAAAGCAAGGCATCCGCTCCCACTGATGACCGGAGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAatat              <  2:184584/147‑1 (MQ=255)
                                                                                                                                            ggTGGGCTAAAATATCCTTATGTAGAAGAATTGCGTCTCGACGAAGCAATGCTCCCGCTCACTCTGATGACCGGAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTa       <  2:179773/151‑1 (MQ=255)
                                                                                                                                                 gCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAggg  >  2:580945/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGG  >  NC_000913/2039929‑2040224

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