Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,807,897 A→G T256A (ACC→GCC)  proX → glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,807,8970AG98.3% 208.2 / ‑4.3 60T256A (ACC→GCC) proXglycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (30/29);  total (31/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 = 3.03e-01

ACACCTGGACGCCGTACTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCACCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTATGCATGACGGCAAAGCCTCAGAAGGCGATATTC  >  NC_000913/2807751‑2808041
                                                                                                                                                  |                                                                                                                                                
acacCTGGACGCCGTACTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCAt                                                                                                                                              <  1:158841/151‑1 (MQ=255)
acacCTGGACGCCGTACTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCAt                                                                                                                                              >  1:88386/1‑151 (MQ=255)
      ggACGCCGTACTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCata                                                                                                                                         <  1:724840/150‑1 (MQ=255)
      ggACGCCGTACTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCat                                                                                                                                          <  1:348204/149‑1 (MQ=255)
          gCCGTACTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCACCATGCATATCGtt                                                                                                                                    >  1:292415/1‑151 (MQ=255)
             gTACTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCAt                                                                                                                                              <  2:815099/138‑1 (MQ=255)
             gTACTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCAt                                                                                                                                              >  1:815043/1‑138 (MQ=255)
             gTACTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGcc                                                                                                                                 >  2:616663/1‑151 (MQ=255)
               aCTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCca                                                                                                                                <  1:711545/150‑1 (MQ=255)
                 tGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCcaac                                                                                                                              >  1:11833/1‑150 (MQ=255)
                  gggTGAGTAACGAACTGAAGCCGTGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCcaac                                                                                                                              <  1:771946/149‑1 (MQ=255)
                  gggTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCcaaca                                                                                                                             <  1:771420/150‑1 (MQ=255)
                  gggTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCcaac                                                                                                                              <  2:646045/149‑1 (MQ=255)
                    gTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAg                                                                                                                          <  2:592734/151‑1 (MQ=255)
                           aCGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGc                                                                                                                   >  1:664002/1‑151 (MQ=255)
                           aCGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGc                                                                                                                   <  2:612395/151‑1 (MQ=255)
                               aCTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCg                                                                                                                 >  2:329911/1‑149 (MQ=255)
                                  gAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGaa                                                                                                             >  2:2328/1‑150 (MQ=255)
                                    aGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGaaaaa                                                                                                          <  2:393606/151‑1 (MQ=255)
                                      ccGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAAcc                                                                                                        >  2:427693/1‑151 (MQ=255)
                                       cGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAAccc                                                                                                       <  2:264097/151‑1 (MQ=255)
                                       cGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAAccc                                                                                                       >  2:98287/1‑151 (MQ=255)
                                         ggCAAAGATGTCGTCTGGTGGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCgg                                                                                                     <  2:641402/151‑1 (MQ=255)
                                             aaGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGgcagc                                                                                                 <  2:603387/151‑1 (MQ=255)
                                              aGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGgcagca                                                                                                <  2:246534/151‑1 (MQ=255)
                                                        aggTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGt                                                                                      >  1:625473/2‑151 (MQ=255)
                                                         ggTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGtt                                                                                     >  2:551304/1‑151 (MQ=255)
                                                          gTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGttt                                                                                    >  2:154811/1‑151 (MQ=255)
                                                            tGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTg                                                                                   <  2:81341/150‑1 (MQ=255)
                                                              cAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGcc                                                                                 >  1:693532/1‑150 (MQ=255)
                                                                 gTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTa                                                                             <  1:816070/151‑1 (MQ=255)
                                                                  tGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTAt                                                                            >  2:146428/1‑151 (MQ=255)
                                                                   gCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGcc                                                                                                                                 <  1:1138/97‑1 (MQ=255)
                                                                   gCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGcc                                                                                                                                 >  2:1138/1‑97 (MQ=255)
                                                                    ccGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTTGGCCGCGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATg                                                                           >  1:439576/1‑150 (MQ=255)
                                                                     cGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCa                                                                         <  2:423152/151‑1 (MQ=255)
                                                                      gTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGc                                                                          >  2:82703/1‑149 (MQ=255)
                                                                             gCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGcc                                                                   >  2:348513/1‑149 (MQ=255)
                                                                               aCTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTg                                                               >  1:492642/1‑151 (MQ=255)
                                                                                   ccGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGc                                                             >  2:271388/1‑149 (MQ=255)
                                                                                      ggCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGata                                                        <  1:87660/151‑1 (MQ=255)
                                                                                        cGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATAtt                                                      <  1:628031/151‑1 (MQ=255)
                                                                                        cGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATAtt                                                      >  2:641531/1‑151 (MQ=255)
                                                                                             aaaaCGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGc                                                 >  1:600740/1‑151 (MQ=255)
                                                                                                aCGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCa                                              <  2:830833/151‑1 (MQ=255)
                                                                                                 cGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGccc                                               >  2:243954/1‑149 (MQ=255)
                                                                                                 cGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGccc                                               >  2:571564/1‑149 (MQ=255)
                                                                                                     gATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACg                                         <  2:504288/151‑1 (MQ=255)
                                                                                                           aaaCTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCAtt                                    >  2:685990/1‑150 (MQ=255)
                                                                                                            aaCTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTAt                                  >  1:296356/1‑151 (MQ=255)
                                                                                                            aaCTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTAt                                  >  1:540478/1‑151 (MQ=255)
                                                                                                            aaCTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTAt                                  <  1:753021/151‑1 (MQ=255)
                                                                                                             aCTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTATg                                 >  1:321940/1‑151 (MQ=255)
                                                                                                             aCTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTATg                                 <  2:632668/151‑1 (MQ=255)
                                                                                                                       ggTGCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTATGCATGACGGc                        <  1:754159/150‑1 (MQ=255)
                                                                                                                          gCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGGAGTTGCCAGAGGCAGAAATTAACGCGCACAACACCATTTAGCATGATGGCa                       >  1:456104/1‑148 (MQ=255)
                                                                                                                          gCGAATTATGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTATGCATGACGGCAAAg                    <  2:563344/151‑1 (MQ=255)
                                                                                                                                 aTGGCTTCCCGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTATGCATGACGGCAAAGCCTCAGa             <  1:220331/151‑1 (MQ=255)
                                                                                                                                          cGGTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTATGCATGACGGCAAAGCGTCAGAAGGCGatat    <  2:693352/151‑1 (MQ=255)
                                                                                                                                            gTCAGCGCCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTATGCATGACGGCAAAGCCTCAGAAGGCGATATTc  >  2:155234/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACACCTGGACGCCGTACTGGGTGAGTAACGAACTGAAGCCGGGCAAAGATGTCGTCTGGTTGCAGGTGCCGTTCTCCGCACTGCCGGGCGATAAAAACGCCGATACCAAACTGCCGAATGGTGCGAATTATGGCTTCCCGGTCAGCACCATGCATATCGTTGCCAACAAAGCCTGGGCCGAGAAAAACCCGGCAGCAGCGAAACTGTTTGCCATTATGCAGTTGCCAGTGGCAGATATTAACGCCCAGAACGCCATTATGCATGACGGCAAAGCCTCAGAAGGCGATATTC  >  NC_000913/2807751‑2808041

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