Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,481,686 T→A 10.9% L133Q (CTG→CAG)  yheS → ABC transporter ATP‑binding protein YheS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,481,6860TA10.9% 179.7 / 9.6 63L133Q (CTG→CAG) yheSABC transporter ATP‑binding protein YheS
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base A (2/5);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCATTGACGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTTGCTGCTGCTCGACGAACCGACTAACCACCTCGATCTCG  >  NC_000913/3481537‑3481835
                                                                                                                                                     |                                                                                                                                                     
tGTCATTGACGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCctgct                                                                                                                                                       <  1:894911/150‑1 (MQ=255)
tGTCATTGACGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCctgct                                                                                                                                                       <  1:170678/150‑1 (MQ=255)
 gTCATTGACGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCctgctg                                                                                                                                                      <  1:1038584/150‑1 (MQ=255)
    aTTGACGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCAc                                                                                                                                                   >  2:553929/1‑150 (MQ=255)
       gACGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCAGCAGGCTGGCAGCACggg                                                                                                                                                <  1:1370305/150‑2 (MQ=255)
       gACGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCAGCAGGCTGGCAGCACggg                                                                                                                                                <  2:452887/150‑2 (MQ=255)
        aCGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGcc                                                                                                                                               <  1:478088/150‑1 (MQ=255)
        aCGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCAGCAGGCTGGCAGCACggga                                                                                                                                               <  1:1169754/150‑3 (MQ=255)
         cGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCt                                                                                                                                              >  2:1367497/1‑150 (MQ=255)
           gCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCAGCAGGCTGGCAGCACgggagcg                                                                                                                                            <  2:913962/150‑6 (MQ=255)
            cGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCgg                                                                                                                                           >  1:928851/1‑150 (MQ=255)
             gACCGTTAATATCGTCAACTATAAGCGCAGCTACACGACGCCAAAGAACGTAACGACGGGCACGCCATTTCGACCATTCATTGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGt                                                                                                                                          <  2:1387241/150‑1 (MQ=255)
               ccGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGttt                                                                                                                                        <  2:1017066/150‑1 (MQ=255)
                 gTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCa                                                                                                                                      <  2:202119/150‑1 (MQ=255)
                 gTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCa                                                                                                                                      <  1:1367497/150‑1 (MQ=255)
                  tGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAg                                                                                                                                     <  1:139063/150‑1 (MQ=255)
                  tGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAg                                                                                                                                     <  2:925060/150‑1 (MQ=255)
                  tGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAg                                                                                                                                     >  1:538769/1‑150 (MQ=255)
                  tGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAg                                                                                                                                     >  2:516882/1‑150 (MQ=255)
                   gAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGc                                                                                                                                    <  1:374603/150‑1 (MQ=255)
                   gAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGc                                                                                                                                    <  1:500920/150‑1 (MQ=255)
                     atatCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCaa                                                                                                                                  <  2:495615/150‑1 (MQ=255)
                       atCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATg                                                                                                                                <  1:553929/150‑1 (MQ=255)
                       atCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATg                                                                                                                                >  2:464352/1‑150 (MQ=255)
                         cGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGaa                                                                                                                              >  1:632868/1‑150 (MQ=255)
                           tCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGaaca                                                                                                                            <  2:817242/150‑1 (MQ=255)
                              aCTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCAGCAGGCTGGCAGCACgg                                                                                                                                                 <  2:12953/126‑1 (MQ=255)
                              aCTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCAGCAGGCTGGCAGCACgg                                                                                                                                                 >  1:12953/1‑126 (MQ=255)
                                      cgcAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTCGAGCGcc                                                                                                                 >  1:759436/1‑150 (MQ=255)
                                       gcAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGccc                                                                                                                <  1:464352/150‑1 (MQ=255)
                                        cAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCg                                                                                                               <  2:1353271/150‑1 (MQ=255)
                                        cAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCg                                                                                                               >  2:65400/1‑150 (MQ=255)
                                              cacGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGt                                                                                                         >  1:199300/1‑150 (MQ=255)
                                                   cgCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGAttt                                                                                                    >  2:701292/1‑150 (MQ=255)
                                                      cAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGAttttta                                                                                                 >  2:1233881/1‑149 (MQ=255)
                                                         cgaacgTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCgg                                                                                              <  2:288249/150‑1 (MQ=255)
                                                         cgaacgTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCgg                                                                                              <  1:611699/150‑1 (MQ=255)
                                                                  cgacgGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCg                                                                                     >  2:1005558/1‑150 (MQ=255)
                                                                    acgGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTa                                                                                   >  1:790669/1‑150 (MQ=255)
                                                                      gggCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATg                                                                                 <  2:928851/150‑1 (MQ=255)
                                                                         cACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGt                                                                              >  2:1382056/1‑150 (MQ=255)
                                                                                  gCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCtt                                                                     >  2:918445/1‑150 (MQ=255)
                                                                                   cGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTg                                                                    >  2:973477/1‑150 (MQ=255)
                                                                                   cGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCAGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTg                                                                    >  1:1228637/1‑150 (MQ=255)
                                                                                    gACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGc                                                                   <  1:65400/150‑1 (MQ=255)
                                                                                      ccATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGccc                                                                 >  2:274434/1‑150 (MQ=255)
                                                                                      ccATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGccc                                                                 <  1:810537/150‑1 (MQ=255)
                                                                                       cATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCa                                                                >  1:513003/1‑150 (MQ=255)
                                                                                            aTGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCg                                                                                     >  1:1108250/1‑124 (MQ=255)
                                                                                            aTGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCg                                                                                     <  2:1108250/124‑1 (MQ=255)
                                                                                            aTGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGgcgc                                                           >  2:1316195/1‑150 (MQ=255)
                                                                                                cAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGAt                                                       >  2:641778/1‑150 (MQ=255)
                                                                                                 aaGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGAtt                                                      <  1:1382056/150‑1 (MQ=255)
                                                                                                           cTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAg                                            >  2:86045/1‑150 (MQ=255)
                                                                                                           cTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAg                                            <  1:973477/150‑1 (MQ=255)
                                                                                                               tGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAg                                            <  1:600267/146‑1 (MQ=255)
                                                                                                               tGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAg                                            >  2:600267/1‑146 (MQ=255)
                                                                                                                   gCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTtgctg                                    >  1:664017/1‑150 (MQ=255)
                                                                                                                    cATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTtgctgc                                   <  1:516882/150‑1 (MQ=255)
                                                                                                                             ttCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTTGCTGCTGCTcgacg                          >  2:1317065/1‑150 (MQ=255)
                                                                                                                                     cGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTTGCTGCTGCTCGACGAACCGACt                  <  1:93606/150‑1 (MQ=255)
                                                                                                                                      gTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTTGCTGCTGCTCGACGAACCGACTa                 <  2:513003/150‑1 (MQ=255)
                                                                                                                                         ctgcCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTTGCTGCTGCTCGACGAACCGACTAacc              <  1:1317065/150‑1 (MQ=255)
                                                                                                                                          tgcCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAg                                            >  1:308609/1‑119 (MQ=255)
                                                                                                                                          tgcCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAg                                            <  2:308609/119‑1 (MQ=255)
                                                                                                                                            ccAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTTGCTGCTGCTCGACGAACCGACTAaccacc           >  2:1250834/1‑150 (MQ=255)
                                                                                                                                                   gctgcACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTTGCTGCTGCTCGACGAACCGACTAACCACctcgatct    >  1:212250/1‑150 (MQ=255)
                                                                                                                                                     tgcACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTTGCTGCTGCTCGACGAACCGACTAACCACctcgatctcg  <  2:538769/150‑1 (MQ=255)
                                                                                                                                                     tgcACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTTGCTGCTGCTCGACGAACCGACTAACCACctcgatctcg  >  2:242942/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTCATTGACGGCGACCGTGAATATCGTCAACTAGAAGCGCAGCTACACGACGCCAACGAACGTAACGACGGGCACGCCATTGCGACCATTCATGGCAAGCTGGATGCTATTGACGCATGGAGTATTCGCTCCCGTGCTGCCAGCCTGCTGCACGGCCTCGGTTTCAGCAATGAACAACTGGAGCGCCCGGTAAGTGATTTTTCCGGGGGCTGGCGTATGCGTCTTAACCTTGCCCAGGCGCTGATTTGCCGTTCAGACTTGCTGCTGCTCGACGAACCGACTAACCACCTCGATCTCG  >  NC_000913/3481537‑3481835

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