Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 326,474 T→C 5.2% A258A (GCA→GCG betA ← choline dehydrogenase, a flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913326,4740TC 5.2% 240.9 / 2.3 77A258A (GCA→GCGbetAcholine dehydrogenase, a flavoprotein
Reads supporting (aligned to +/- strand):  ref base T (33/40);  new base C (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

ATATTGCAGATACATCTCCAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGG  >  NC_000913/326330‑326614
                                                                                                                                                |                                                                                                                                            
atatTGCAGATACATCTCCAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGcac                                                                                                                                          <  5:294233/149‑1 (MQ=255)
   ttGCAGATACATCTCCAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACAt                                                                                                                                        >  6:124098/1‑148 (MQ=255)
   ttGCAGATACATCTCCAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACAt                                                                                                                                        >  6:265241/1‑148 (MQ=255)
              tctcCAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATACCAGCActtct                                                                                                                            >  1:86336/1‑149 (MQ=255)
               ctcCAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCGGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCActtctt                                                                                                                           >  8:242436/1‑149 (MQ=255)
               ctcCAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCActtctt                                                                                                                           >  2:151481/1‑149 (MQ=255)
                tcCAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCttt                                                                                                                          >  4:196720/1‑149 (MQ=255)
                tcCAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCttt                                                                                                                          >  4:12411/1‑149 (MQ=255)
                  cAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTttg                                                                                                                         <  7:405216/148‑1 (MQ=255)
                    gatgatCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTttgttg                                                                                                                      <  2:100593/149‑1 (MQ=255)
                          ccTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGtt                                                                                                                >  1:320665/1‑149 (MQ=255)
                           cTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCTGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTATTTGTTGGCCGtt                                                                                                                <  3:170661/148‑1 (MQ=255)
                           cTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGtt                                                                                                                <  5:197883/148‑1 (MQ=255)
                           cTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTg                                                                                                               <  1:219411/149‑1 (MQ=255)
                            tGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTgc                                                                                                              >  3:260539/1‑149 (MQ=255)
                            tGAAGATTTTCGCCGAAGCCGGGTAATTCATGCGCCAGCGGAATATCAAACTCCGCCAGCACTTCCGCGTTGCCGACGCCGGAGCGTTGCAGGAACTGCGGTGAGGCAATCGCGACTGCACATAACAGCACTTCTTTGTTGGCCGTTgc                                                                                                              >  8:312652/1‑149 (MQ=255)
                             gAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCCGCACATAACAGCACTTCTTTGTTGGCCGTTgcg                                                                                                             >  3:428548/1‑149 (MQ=255)
                                 aTTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCggg                                                                                                         <  8:60254/149‑1 (MQ=255)
                                    ttCGCCGACGCCGGGTAATTCATGCACCAGGGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGtt                                                                                                       <  5:124098/148‑1 (MQ=255)
                                    ttCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCCCTTCTTTGTTGGCCGTTGCGCGGGTTg                                                                                                      >  7:209788/1‑149 (MQ=255)
                                    ttCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTg                                                                                                      >  2:144230/1‑149 (MQ=255)
                                     tCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTg                                                                                                      <  7:256350/148‑1 (MQ=255)
                                     tCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTg                                                                                                      <  1:151481/148‑1 (MQ=255)
                                       gCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTggg                                                                                                    <  6:435031/148‑1 (MQ=255)
                                       gCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTggg                                                                                                    <  1:144230/148‑1 (MQ=255)
                                        ccGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGAt                                                                                                  <  6:178606/149‑1 (MQ=255)
                                            cGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTg                                                                                              <  3:12411/149‑1 (MQ=255)
                                               cGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTg                                                                                           >  3:386616/1‑149 (MQ=255)
                                                 ggTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCGGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTc                                                                                         <  7:242436/149‑1 (MQ=255)
                                                      ttCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCtt                                                                                    >  7:177154/1‑149 (MQ=255)
                                                         aTGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCa                                                                                 <  5:265241/149‑1 (MQ=255)
                                                         aTGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCa                                                                                 <  5:197776/149‑1 (MQ=255)
                                                         aTGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCa                                                                                 <  7:108553/149‑1 (MQ=255)
                                                          tGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAg                                                                                <  2:146387/149‑1 (MQ=255)
                                                           gCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGc                                                                               >  3:16284/1‑149 (MQ=255)
                                                                  cGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGa                                                                        >  3:376578/1‑149 (MQ=255)
                                                                       tatCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGcc                                                                    >  1:187542/1‑148 (MQ=255)
                                                                          cAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCTTTGGAGGATCTGCGTTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTTTCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTcc                                                                                  <  1:128794/131‑1 (MQ=255)
                                                                          cAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCTTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTTCGCGGGTTGGGATGGTGCTGTCGCCTTcc                                                                                  >  2:128794/1‑131 (MQ=255)
                                                                            aaCTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCCGTCGCCTTCCAGCCATTCGACGCCCACcgc                                                              >  4:32050/1‑149 (MQ=255)
                                                                             aCTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCCCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCATTCGCCTTCCAGCCATTCGACGCCCCCcgcg                                                             >  4:242717/1‑149 (MQ=255)
                                                                               tCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGCCGCCCACcgcgcg                                                           >  2:447274/1‑149 (MQ=255)
                                                                               tCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACcgcgcg                                                           >  8:322692/1‑149 (MQ=255)
                                                                               tCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATAGCGCCTTCACATAACATCACATCATTGTTGGCCGTTGCGCGGGTTGGGATGCTGCTGTCGCCTTCCAGCCATTCGACTCCCACcgcgcg                                                           >  2:103720/1‑149 (MQ=255)
                                                                                 cGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACcgcgcg                                                           <  3:196720/147‑1 (MQ=255)
                                                                                  gCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGAGAGGCAATCGCGCCTGCACTTAACAGCACTTCTTTGATGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGTCGCCCACCGCGCGttt                                                        <  7:312652/149‑1 (MQ=255)
                                                                                    cagcagTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGc                                                      >  1:176378/1‑149 (MQ=255)
                                                                                     agcagTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGcc                                                     <  4:376578/149‑1 (MQ=255)
                                                                                           tCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGcc                                                     <  1:213662/143‑1 (MQ=255)
                                                                                           tCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGcc                                                     >  2:213662/1‑143 (MQ=255)
                                                                                            cAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCCGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCaaa                                               >  7:265980/1‑148 (MQ=255)
                                                                                               cGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAatg                                           <  2:176378/149‑1 (MQ=255)
                                                                                                gTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAatga                                          <  2:187542/149‑1 (MQ=255)
                                                                                                 ttGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAatgat                                         <  2:320665/149‑1 (MQ=255)
                                                                                                 ttGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAatgat                                         <  4:260539/149‑1 (MQ=255)
                                                                                                   gCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAatgatg                                        <  7:322692/148‑1 (MQ=255)
                                                                                                           cGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGt                                >  8:347712/1‑148 (MQ=255)
                                                                                                                  ttGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGt                        <  8:177154/149‑1 (MQ=255)
                                                                                                                  ttGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGt                        >  7:79265/1‑149 (MQ=255)
                                                                                                                     cAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGa                      <  4:386616/148‑1 (MQ=255)
                                                                                                                     cAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGa                      <  7:247631/148‑1 (MQ=255)
                                                                                                                     cAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAg                     >  2:221326/1‑149 (MQ=255)
                                                                                                                      aGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGt                    <  4:16284/149‑1 (MQ=255)
                                                                                                                       ggATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTa                   <  8:79265/149‑1 (MQ=255)
                                                                                                                         aTCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAat                                            <  4:226928/122‑1 (MQ=255)
                                                                                                                         aTCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAat                                            >  3:226928/1‑122 (MQ=255)
                                                                                                                              cGGTGAGGCAATCGCGCCCGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGaa               <  4:428548/146‑1 (MQ=255)
                                                                                                                               ggTGAGGCATTCGCGCCTGCACATAAGAGCTCTTCTTTGTTGGCCGTTGCGCGGGTTGGTATTGTGCTGTCGCCTTCCAGGCATTCGACTCCCACCGCGAGTTTGACGTCAAAATTGATGTGATCGGTCATAGCGTGTGTACGAATGGt           <  1:103720/149‑1 (MQ=255)
                                                                                                                                  gAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAg        <  1:285719/149‑1 (MQ=255)
                                                                                                                                    ggCAATCGCGCCCGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTTATCGGTCATAGCGTGAGTACGAATGGTCAGGt      <  8:265980/149‑1 (MQ=255)
                                                                                                                                     gCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGt      >  5:46184/1‑148 (MQ=255)
                                                                                                                                     gCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGt      >  8:160494/1‑148 (MQ=255)
                                                                                                                                       attCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAg   <  3:242717/147‑1 (MQ=255)
                                                                                                                                       aaTCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAg   >  7:41876/1‑149 (MQ=255)
                                                                                                                                        aTCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAgg  <  2:86336/149‑1 (MQ=255)
                                                                                                                                               cTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTc                                                  <  4:424092/94‑1 (MQ=255)
                                                                                                                                               cTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTc                                                  >  3:424092/1‑94 (MQ=255)
                                                                                                                                                |                                                                                                                                            
ATATTGCAGATACATCTCCAGATGATCCTGAAGATTTTCGCCGACGCCGGGTAATTCATGCACCAGCGGAATATCAAACTCCGCCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGG  >  NC_000913/326330‑326614

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