Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,985,112 C→G 16.9% intergenic (‑43/+27) araG ← / ← araF L‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,985,1120CG16.9% 231.0 / 39.9 83intergenic (‑43/+27)araG/araFL‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (40/29);  new base G (10/4);  total (50/33)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

ACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAA  >  NC_000913/1984971‑1985260
                                                                                                                                             |                                                                                                                                                     
aCCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga                                                                                                                                                <  3:456854/149‑1 (MQ=255)
  cGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑gg                                                                                                                                             <  2:491803/149‑1 (MQ=255)
  cGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑gg                                                                                                                                             <  3:684106/149‑1 (MQ=255)
   ggCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑ggg                                                                                                                                            <  3:401845/149‑1 (MQ=255)
        aGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAAtt                                                                                                                                       >  1:610663/1‑149 (MQ=255)
            gTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTc                                                                                                                                   <  2:231126/149‑1 (MQ=255)
             tCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga                                                                                                                                                <  2:357327/136‑1 (MQ=255)
             tCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga                                                                                                                                                >  1:357327/1‑136 (MQ=255)
               aaaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCCTCATGCAGA‑GGGGAATTTTTCCgg                                                                                                                                >  1:304134/1‑149 (MQ=255)
               aaaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑gg                                                                                                                                             >  3:307067/1‑136 (MQ=255)
               aaaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑gg                                                                                                                                             <  4:307067/136‑1 (MQ=255)
               aaaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCgg                                                                                                                                >  1:435181/1‑149 (MQ=255)
               aaaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCgg                                                                                                                                >  4:520476/1‑149 (MQ=255)
                       tatCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAAttac                                                                                                                        <  3:655239/149‑1 (MQ=255)
                       tatCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAAttac                                                                                                                        >  2:364915/1‑149 (MQ=255)
                           cGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑gggg                                                                                                                                           >  4:85757/1‑126 (MQ=255)
                           cGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑gggg                                                                                                                                           <  3:85757/126‑1 (MQ=255)
                            gTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTAcc                                                                                                                   >  2:286099/1‑149 (MQ=255)
                            gTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTAcc                                                                                                                   >  1:633656/1‑149 (MQ=255)
                                    cTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTAACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAAc                                                                                                           <  1:27816/149‑1 (MQ=255)
                                     ttAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAAcc                                                                                                          >  3:588105/1‑149 (MQ=255)
                                         cGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCtttt                                                                                                      <  3:520476/149‑1 (MQ=255)
                                          gCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttttt                                                                                                     >  3:309287/1‑149 (MQ=255)
                                                aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttt                                                                                                       >  4:344268/1‑141 (MQ=255)
                                                aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttt                                                                                                       <  3:344268/141‑1 (MQ=255)
                                                    gTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttc                                                                                           >  1:604581/1‑149 (MQ=255)
                                                       ttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc                                                                                        <  4:309287/149‑1 (MQ=255)
                                                                ccGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGtt                                                                               >  4:610691/1‑149 (MQ=255)
                                                                        tGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                                                                       >  3:151475/1‑149 (MQ=255)
                                                                        tGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                                                                       >  1:453488/1‑149 (MQ=255)
                                                                            agaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGAt                                                                   <  2:669208/149‑1 (MQ=255)
                                                                            agaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGAt                                                                   >  1:454717/1‑149 (MQ=255)
                                                                            agaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGAt                                                                   <  1:286099/149‑1 (MQ=255)
                                                                                aCGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt                                                               >  2:383915/1‑149 (MQ=255)
                                                                                       aGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGt                                                        >  3:396845/1‑149 (MQ=255)
                                                                                        gACTGTTGCATAACGTCTCGGTGAATCACTGGTCGTTCACACCCCCTCTGCACCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTc                                                       >  4:585962/1‑149 (MQ=255)
                                                                                        gACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTc                                                       >  1:376536/1‑149 (MQ=255)
                                                                                          cTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg                                                     <  2:453488/149‑1 (MQ=255)
                                                                                            gTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtctc                                                                                                 >  2:377148/1‑103 (MQ=255)
                                                                                            gTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtctc                                                                                                 <  1:377148/103‑1 (MQ=255)
                                                                                            gTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTa                                                   >  3:109888/1‑149 (MQ=255)
                                                                                               gCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACt                                                >  3:332088/1‑149 (MQ=255)
                                                                                                  tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCg                                             <  1:383915/149‑1 (MQ=255)
                                                                                                   aaCGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCgg                                            <  4:396845/149‑1 (MQ=255)
                                                                                                   aaCGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCgg                                            >  1:363349/1‑149 (MQ=255)
                                                                                                              gAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTgg                                  >  4:527845/1‑148 (MQ=255)
                                                                                                              gAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTgg                                  <  3:527845/148‑1 (MQ=255)
                                                                                                               aaTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCg                                >  1:678182/1‑149 (MQ=255)
                                                                                                                 tCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaaa                                                                                  <  2:82828/97‑1 (MQ=255)
                                                                                                                 tCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaaa                                                                                  >  1:82828/1‑97 (MQ=255)
                                                                                                                    cTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCa                           >  2:301878/1‑149 (MQ=255)
                                                                                                                     tGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCaa                          <  4:109888/149‑1 (MQ=255)
                                                                                                                     tGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCaa                          <  1:301878/149‑1 (MQ=255)
                                                                                                                      ggTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAAc                         <  1:265528/149‑1 (MQ=255)
                                                                                                                       gTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAACACGTCGGTAACATCGGTAAATTTTGGCGGTTCaaag                        >  3:551972/1‑147 (MQ=255)
                                                                                                                        tCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttc                                                                                           <  2:135056/81‑1 (MQ=255)
                                                                                                                        tCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttc                                                                                           >  1:135056/1‑81 (MQ=255)
                                                                                                                        tCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGt                       <  2:610663/149‑1 (MQ=255)
                                                                                                                            tCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCttt                   >  4:480517/1‑149 (MQ=255)
                                                                                                                             cacaCCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtttt                  >  3:374033/1‑149 (MQ=255)
                                                                                                                              acacCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTgg                                  <  3:459803/132‑1 (MQ=255)
                                                                                                                              acacCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTgg                                  >  4:459803/1‑132 (MQ=255)
                                                                                                                                ttccccTCTGCATGATGCAGAGGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGc                >  1:662215/3‑149 (MQ=255)
                                                                                                                                acCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCt               >  1:292470/1‑149 (MQ=255)
                                                                                                                                 tccccTCTGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTa              >  2:622091/2‑149 (MQ=255)
                                                                                                                                 cccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGt                              >  3:510220/1‑133 (MQ=255)
                                                                                                                                 cccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTa              <  4:332088/149‑1 (MQ=255)
                                                                                                                                 cccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCCGTAAATTTTGGCGGt                              <  4:510220/133‑1 (MQ=255)
                                                                                                                                  ccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGt                                                                      >  3:473189/1‑92 (MQ=255)
                                                                                                                                  ccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGt                                                                      <  4:473189/92‑1 (MQ=255)
                                                                                                                                   cccTCTGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAcc            >  3:279176/1‑149 (MQ=255)
                                                                                                                                   cccTCTGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAcc            >  2:524211/1‑149 (MQ=255)
                                                                                                                                   cccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAAtttt                                    <  2:383068/125‑1 (MQ=255)
                                                                                                                                   cccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAAtttt                                    >  1:383068/1‑125 (MQ=255)
                                                                                                                                    ccTCTGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc           >  4:350113/1‑149 (MQ=255)
                                                                                                                                    ccTCTGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc           >  1:317183/1‑149 (MQ=255)
                                                                                                                                    ccTCTGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc           >  4:103486/1‑149 (MQ=255)
                                                                                                                                     ctctGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCa          >  3:212338/1‑149 (MQ=255)
                                                                                                                                     ctctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTc                      >  2:544484/1‑137 (MQ=255)
                                                                                                                                     ctctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTc                      <  1:544484/137‑1 (MQ=255)
                                                                                                                                     ctctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCa          <  1:364915/149‑1 (MQ=255)
                                                                                                                                     ctctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCa          >  1:510553/1‑149 (MQ=255)
                                                                                                                                      tctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAg         >  1:685187/1‑149 (MQ=255)
                                                                                                                                             catGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTaaa  <  4:151475/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
ACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAA  >  NC_000913/1984971‑1985260

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