Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,975,091 T→G 18.6% T308P (ACT→CCT)  aas ← fused 2‑acylglycerophospho‑ethanolamine acyl transferase/acyl‑acyl carrier protein synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,975,0910TG18.6% 108.3 / 2.5 70T308P (ACT→CCT) aasfused 2‑acylglycerophospho‑ethanolamine acyl transferase/acyl‑acyl carrier protein synthetase
Reads supporting (aligned to +/- strand):  ref base T (23/30);  new base G (11/6);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGGACGACGCCTTTCGGATGGCCTTCAGAACCGGAGGTAAAAA‑GGATCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAACACTGTATTTTTCAAGGATGCGTCCAACAAACAGCGTTTTCGTAAGCAATTTGCGATAGGAGTCTGGGG  >  NC_000913/2974857‑2975324
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
tATGGACGACGCCTTTCGGATGGCCTTCAGAACCGGAGGTAAAAA‑GGATCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAGGGGaaa                                                                                                                                                                                                                                          >  2:855794/1‑238 (MQ=255)
 attggccgggCCTTTCGGATGGCCTTTAGAACCGGAGGTAAAAA‑GGATCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTTGG‑CCAGACGCGGCATCAGCAAATGAGCGAAGAGCCATACTTTGGCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGa                                                                                                                                                                                                                                         <  2:433525/230‑1 (MQ=255)
 aTGGACGACGCCTTTCGGATGGCCTTCAGAACCGGAGGTAAAAA‑GGATCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCCGCTAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGCAA‑CTGGCGGGAAGGGAAGa                                                                                                                                                                                                                                         >  2:384861/1‑239 (MQ=255)
  tGGACGACGCCTTTCGGATGGCCTTCAGAACCGGAGGTAAAAA‑GGATCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGCGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGGGAAA‑CTGGCGGGGAGTGACGAt                                                                                                                                                                                                                                        >  2:131570/1‑239 (MQ=255)
  tGGACGACGCCTTTCGGATGGCCTTCAGAACCGGAGGTAAAAA‑GGATCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGTCGGGAAGTGAAGAt                                                                                                                                                                                                                                        >  2:877963/1‑239 (MQ=255)
         cgCCTTTCGGATGGCCTTCAGAACCGGAGGTAAAAA‑GGATCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGGGAAA‑CTGGCGGGAAGTGAAGATGGgttgg                                                                                                                                                                                                                                 >  2:668662/1‑237 (MQ=255)
                 gggttgCCTTCAGAACCGGGGGTAAAAG‑GTGTCCGCGCCTTTTCTTCCGGCTGCTTTTT‑AACGTGTG‑CCAGCCGGGGCACAAGCAAAGGAGGGAAGGTCAATGCTTTGTCGGCAGTGGT‑GAAATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCCCCTGAGTAAGTTGCTCCGGAAGATGCCAGATTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑gc                                                                                                                                                                                                                        <  1:851961/234‑1 (MQ=255)
                  gaggggcTTAAGACCCGGAGGTAAAAA‑GGATACGCGCCTTTTTTTCCGGCTGCTGTTT‑AACCTTGT‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGGGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCGGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑gc                                                                                                                                                                                                                        <  2:680120/232‑1 (MQ=255)
                  gATGGCCTTGAAAACCGGGGGAAAAAA‑GGACCAGCGCCTCTCTTTCCGGCTGGTGTTT‑ACCCTGTG‑CAGGACGCGGCAGCAGAAAAGGAGCGAAGATCCATACTTTGTCGGCAGGGGT‑GAGTGTTGCTTTTAAAT‑CTGCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCGGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGGGGGAAGTGAAGATGGTTTTGATTTCA‑gc                                                                                                                                                                                                                        <  1:131570/238‑1 (MQ=255)
                     ggctttccaaccgggaggaAAAAG‑GGACCAGGCCCTTTTTCTCCGGCT‑CTTGTT‑CACCCGTGCCCGGACGCGGCATCGCGAAAGGAGCAAGGATCAATACTTTGGCGGCAGTGGT‑AACATCTG‑TTTTAAATGCTCCCAGATAGCCCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCAGGAGTTTGCCTTTACCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑gccgc                                                                                                                                                                                                                     <  1:171322/225‑1 (MQ=255)
                      gCCTTCAGAACCGGAGGTAAAAA‑GGATCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTTGCGGGAAGTGAAGAGGGGTTTGATTTCA‑Gcagccg                                                                                                                                                                                                                   >  2:565276/1‑239 (MQ=255)
                         ctgaaaaCCGGGGGGAAAAAGGGATCAGCGCCTGTTCTGCCGGCTGCTGTTT‑AACCTTGG‑CCAGCCGGGTCATCAGCAAATGAGCGAAGCTCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCCCCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGt                                                                                                                                                                                                                  <  2:287898/233‑1 (MQ=255)
                              acgcGGAGGTAAAAG‑GGTTCTGCCCTTTTTCTCCGGGTTGTTGTTT‑ACCTGGTG‑CAGGGGGCGGCTAAAAAAATTGAGCGAAGATCAATACTTTGGGGGCGGTGGG‑TACATGTTCTTTAAAAT‑CTTCCAGATAGCCCCAGCGCCCTTGAGTAAGTTGCGCCGGCAGATGCCAGTGTTTGCCTTTATCGAGAAA‑CGGGCGGGAAGTGAAGTTGGTTTTGATTTCA‑GCCGCCGTAATAGCa                                                                                                                                                                                                           <  1:214776/236‑1 (MQ=255)
                              aaCCGGAGGAAAAAG‑GGTACAGCGCCTCTTCTTCCGGCTGCTGTTT‑ACCCTGTC‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGCCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGTGTAAGTTGCTCCGGCAGAGGCCAGGGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTGCA‑GCCGCCGTAATAGCa                                                                                                                                                                                                           <  2:22273/239‑1 (MQ=255)
                                   gAGGTAAAAA‑GGATCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGGGAAGATGGGTTTGGTTTCA‑GCCGCCGGAATAGCACTgg                                                                                                                                                                                                       >  2:292484/1‑238 (MQ=255)
                                      gatgtgtaaaagggacagctctTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGCCGCGGCATCAGCAAAGGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGGGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTc                                                                                                                                                                                                     <  2:200258/220‑1 (MQ=255)
                                       tAAAAG‑GGATCGGCGCCTCTTCTTCCGGCTGTTTTTA‑AACCTGGC‑CCAGACGCGGCATCAGAAAATGAGGAAAGATCCATTCTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAg                                                                                                                                                                                                   <  2:522162/238‑1 (MQ=255)
                                       tAAAAA‑GGATCAGCGCCTCTTCCTCCGGCTGCTGTTT‑AACCTGGG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGCCCCAGCGCCCCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAg                                                                                                                                                                                                   <  2:62650/238‑1 (MQ=255)
                                        aaaaaaggtgctcccccTCTTCTTCGGGTTTCTTTTT‑AACCGGTG‑CCGGAGGCGGCATCAGCAAATGAGCGAAGATCCATATTGTGCCGGCGGGGGT‑GACATCTGCTTTAAAAT‑CTCCAAAATAGACCCAGCGCCCCGGATTAAGTTGCTCCGGCAGATGCCAGAGTTTGCTTTTATCGAGAAC‑CTGGCGGAAAGTGAAGAGGGTTTTGATTCCA‑GCCGCCGTAATAGCACTGGTCAg                                                                                                                                                                                                   <  1:476575/223‑1 (MQ=255)
                                                 tCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGG‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGGAGGGAAGATGGCTTTTATTTTA‑ACCCCCGGAAAAGAACAGGGCAAGCCTTTTaa                                                                                                                                                                                          >  1:240760/1‑237 (MQ=255)
                                                     cgcCTCTTCTCCGGGCTGCTGTTT‑ACCCTGGG‑CCAGACGCGGCATCAGAAAATGAGAGAAGCTCCATACTTTGCCGGCAGGGGT‑GACATTTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGAGGCCAGATTTTGCCTTTATCGAGAAA‑CGGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCgg                                                                                                                                                                                     <  2:824741/239‑1 (MQ=255)
                                                       ccTCTTCCTGCGCCTGTTTTTA‑ATCCGGCG‑CCAGACGCCGCATCAGAAATTGCACGAAGCTCCATCCTTGCGCGGCAGGGGA‑GACATGTGCTTGTAAAT‑CTGCCAGATAGCCCCAGCGCGCTTGAGTAAGTTGCTCCGGCAGATGCCAGTGTTTGCTTTTATCGAAAAA‑CTGGAGGGAAGGGAAGATGGTTTTGTTTTAA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCa                                                                                                                                                                                   <  1:668662/239‑1 (MQ=255)
                                                         ttttcttcCGGCTGCTTTTT‑ACCCTGTG‑CCAGACGCGGCATCAGAAAATGAGCGAAGCTCCATACTGTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTCCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGTTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAgt                                                                                                                                                                                 <  2:168925/237‑1 (MQ=255)
                                                          tttcttcCGGTTG‑TGTTTCACCCTGTG‑CCAGAGGCGGCATCAACAAAGGAGCGAAAATCCATCTTTTGTCGGCGGTGTT‑GACTTCTGTTTTAAAAT‑CTTCAAAATAGACCCAGCCCCCCTAGGTAGGTGGCTCCGGAAGAGGCCAGGTTTTGCCTTTACCGAAAAA‑CGGGGGGGAAGGGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAgtg                                                                                                                                                                                <  1:557718/238‑1 (MQ=255)
                                                          cttctcCCGGCTGTTGTTT‑ACCCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGACCCATACTTGGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTGGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAgtg                                                                                                                                                                                <  2:584271/239‑1 (MQ=255)
                                                             cttcCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CCGCCGGGGAGGGAAGATGGTTTTGGTTTTT‑GCCCCCGCATTTGCACTGGTCACCCCCTTTACCCCCGCAgggtca                                                                                                                                                                             >  1:186752/1‑237 (MQ=255)
                                                                   gctgctGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGGGAAA‑CTGGCGGGAGGGGAAGATGGGTTTGATTTCA‑GCCGCCGGAAAAGAACCGGGCAGCCCTTTTa                                                                                                                                                                                           >  1:200258/1‑219 (MQ=255)
                                                                          ttt‑AACTGGGG‑CCAGACGGGGCATCAGCATAGGAGCGAAGTGCCAGACTTTGTGGGCAGTGGT‑GAAATCTGGTTTAAAAT‑CTGCCAGAGAGGCCCAGCACACGGGAGTAAGTTGCGCGGGCAAATGCAAGATTCTGTCTTTATGGAAAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTa                                                                                                                                                                              <  2:324767/225‑1 (MQ=255)
                                                                                   gtg‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGGTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGGGTAGTTCATCATTTGGGGCAttggg                                                                                                                                                        >  2:456252/1‑235 (MQ=255)
                                                                                     g‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGTAGTGCAGATGGTTTTGATTTTC‑ACCACAGTAAAAACCATCGGCAGCCCCTTTATCCCCGCAGGGTAGTTCTTCACTTCCGGCACgcggcc                                                                                                                                                      >  1:633509/1‑238 (MQ=255)
                                                                                     g‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGGTTTGATTTTA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTTCGGGCATgcggcg                                                                                                                                                      >  2:607649/1‑239 (MQ=255)
                                                                                               ggCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGGAATAGCCCTGGTCCGCCCTTTTTCCCCCGCCTGGTAGTTCATCTTTGTGGGGATGCGGGGA‑Cgg‑gggca                                                                                                                                           >  1:728630/1‑235 (MQ=255)
                                                                                                         aaggAGCGAATATCCATACTTTGCGGGCA‑TGGTGGACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCaaaa                                                                                                                                 <  2:231937/236‑1 (MQ=255)
                                                                                                         aaTGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGGAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCCCTGGTCAGCCCTTTTACCCCGGCCGGGTAGTTCATCATTCCGGGCATCCGGCGA‑ACG‑GCGCTGGGCCCaaaa                                                                                                                                 >  2:729505/1‑239 (MQ=255)
                                                                                                          aTGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTTCGGGCATGCGGGGC‑CGG‑GCGCTGGCCCCcaaaa                                                                                                                                >  2:821891/1‑239 (MQ=255)
                                                                                                           tGAGGGAACATCTATGTTTTGTGGGCAGGGGG‑GACATTTTCTTTAAAAT‑CTTCAAGATAGCCCCACCGCACGGGGGTAAGTGGCTCCGGCAGAGGCCATAGTTTCCTTTTATCGAGAAA‑CTGGGGGGAAGTGAAGATGGTTTTGTTTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAAt                                                                                                                               <  1:821891/239‑1 (MQ=255)
                                                                                                               cGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGTAGTGAAGATGGTTTTTATTTTA‑GCCGCCGTAAAAGCCCTGGTCAGCCCCTTTACCCCGGCATGGTAGTTCACCATTTCGGGACTGCGGGGAGCGG‑GCGA‑GGCCCCacaaaaaca                                                                                                                            >  1:580883/1‑234 (MQ=255)
                                                                                                               cGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGGGAAGATGGGTTTGATTTCC‑GCCCCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCCGGGTAGTTCATCATTTCGGGCCTGCGCCAC‑CGGAGCGTGGCCCCCAAAAACCa                                                                                                                             >  1:194905/1‑238 (MQ=255)
                                                                                                                 aaGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGGGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGGGTAGTTCATCATTCCGGGCCTTCGGCCA‑CGG‑GGGATGGCCCCCAAAATCACTg                                                                                                                          >  2:479329/1‑238 (MQ=255)
                                                                                                                    aTCCAGACTTTGGCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGAGAGACCCAGCGCACCTGAGTAAGTGGCTCCGGCAGAGGCCAGAGTTTGCCTTTACCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTc                                                                                                                                                                                                     <  2:60509/162‑1 (MQ=255)
                                                                                                                       ctttttttGTCGGGAGTGTT‑GACATTTGTTTTAAAAT‑CTCCCAAGTAGCCCCGGCGCACGGGAGTAAGTCGCTCCGGGAGATGCCAGTGTTTGCCTTTAGCGAGAAA‑CGGGCGGGAAGTGAAGGTGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACt                                                                                                                   <  1:260230/234‑1 (MQ=255)
                                                                                                                       catcatttGGCGGGAGGGGT‑GACATCTGTTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAATAAGTTTCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACt                                                                                                                   <  2:187006/234‑1 (MQ=255)
                                                                                                                       cATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGGGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGGCGGGTAGGTTCTCATTCCGGGCGTGCGGCCA‑CGG‑GCGATGGGCCCaaaaaaccacgccccccc                                                                                                                   >  1:413831/1‑226 (MQ=255)
                                                                                                                              tGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTTATTTCC‑GCCGCCGTAATAGCCCTGGTCAGCCCTTTTACCCCGGCAGGGTAGTTCATCATTGCCGGCATGCGGCGA‑CGG‑GCGATGGGCCCAAAAATCACTGCCGCCatgatg                                                                                                               >  2:568274/1‑236 (MQ=255)
                                                                                                                              tGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAAAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCC‑GCCCCCGTAATAGCACTGGTCAGCCCTTTTTACCCGGCAGGGTAGTTCCTCATTGCGGGCATGGCGGGC‑CGG‑GGGGTGGCCCCAAAAATCACCCCCGccccggtg                                                                                                               >  1:520754/1‑229 (MQ=255)
                                                                                                                                tCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATATACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGTGAAAGGAAGATGGCTTTGATTTCC‑GCCGCCCTAATAGCACCTGTCAGCCCTTTTACCCCGGCAGAGTAGTTCATCCTTGCCTGCACTCGGCCA‑CGG‑GCCATGGCCCCCAACATCACTGTCCCACTGAAGcacgc                                                                                                          >  2:655841/1‑235 (MQ=255)
                                                                                                                                    cAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCCCTGATGCCTGCTTTg                                                                                                      >  2:537998/1‑239 (MQ=255)
                                                                                                                                       tggt‑gCCTTCTGCTTTGAAAT‑CTTCCAGAGAGCCCCAGCGCACCTGAGTAAGTTGCTCCGGCAGAGGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTggg                                                                                                    <  2:588939/238‑1 (MQ=255)
                                                                                                                                       gggt‑aACATCTCCTTTTAAAT‑CTTCCAGATAGACCCGGCGCACCTAAGTACGTTCCTCCGGCAGTTCCCAGAGTTTGCCTTTATGGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTCCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTggg                                                                                                    <  2:768078/237‑1 (MQ=255)
                                                                                                                                         gt‑gACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGGTTTGATTTCA‑GACGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGGTCATAATTTCGGGAATGCGGCCA‑CGG‑GCGAGGGCCCCAAAAATCACTCCCGCCCTTAGTCCTGCATTGGgcag                                                                                                 >  2:115554/1‑239 (MQ=255)
                                                                                                                                                       ttAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGGTTTGATTTCA‑GCCGCCGTAATAGCACTGGCCAGCCCTTTTACCCCGGCAGGGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGTAGGCCCCAAAAATCACTGCCCACCTGATGCATGCATTGGGCGGCCTTAGACCGATg                                                                                    >  2:126093/1‑239 (MQ=255)
                                                                                                                                                             tcctCCAGATAGACCCAGCGTCCCTGAGTAAGTTGCTCCGGCAGATGCCAGATTTTGCCTTTAGCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTc                                                                               <  2:23676/236‑1 (MQ=255)
                                                                                                                                                             tcctCCAGAAAGGCCCAGGCCCCCTGAGTAAGTTGCTCCGGGAGATGCCAGGGTTTCCCTTTAGCGAGAAA‑CGGGCGGGAAGGGAAGTTGGTTTTGATTCCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTc                                                                               <  2:767852/236‑1 (MQ=255)
                                                                                                                                                                  cagttacgcccagcccagCTGAGTATGGTGCCCGGGGAGATGAATGTTTTTGCTTTGATCGAGACA‑CGGGGGGGAAGGGATGGTGGTTTTGTTTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCtt                                                                          <  1:877963/231‑1 (MQ=255)
                                                                                                                                                                       tgcccctcccccttggtttatgTTGCTCCGGCAGATGCCAGAGTTTGCCTTTACCGAGAAT‑CGGGCGGGAAGTGAAGGTTGTTTTGATTCCC‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAaca                                                                     <  1:104881/218‑1 (MQ=255)
                                                                                                                                                                       tacccctccccactgaaGTATGTTTGTCCGGGAGATCCTAGTGTTTTCCTTAACGGTGAAA‑GGGGGGGGAAGTAAAGTTGGTTTGGTTTCGA‑GCCGCCTTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAaca                                                                     <  1:860797/223‑1 (MQ=255)
                                                                                                                                                                        gACCCAGGCGACGTGAGTAGTGTGCACCGGCAGATGCCAGATGTTGCCTTTATCGAGAAA‑CTGGCGGGAGTGGAAGATGTTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAaca                                                                     <  2:441798/238‑1 (MQ=255)
                                                                                                                                                                         aCCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGCAGTGAAGATGGTTTTGATTTCA‑GCCCCCGTAATAGCCCTGGTCAGCCCTTTTAACCCGGCCGTGTAGTTCCTCATTTCGGGCATTCGGCGA‑CGG‑GGCATGGCCCCAAAAATCACTGCCGCACGTGTGCCCGCACTGGGCGGCTTTAAGCCGATGGGTTGGCCCTCacaac                                                                    >  1:767852/1‑234 (MQ=255)
                                                                                                                                                                         aCCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGACGCACTGATGCCTGCCTTGGGCAGCATTAAGCCGATGCGTTCGCCTTTAACACt                                                                   >  2:539197/1‑239 (MQ=255)
                                                                                                                                                                              gcgcACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCCGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCACAAATCACTGCCGCACTGATGCCTGCATTGGGCAGACTTAAGCCGAAGCGTTCGCCttataacccgtatt                                                              >  2:319035/1‑227 (MQ=255)
                                                                                                                                                                               cgcACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGGCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTTCCCTTCAACACTGTAttt                                                             >  2:342297/1‑239 (MQ=255)
                                                                                                                                                                               cgcACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAACACTGTAttt                                                             >  2:171060/1‑239 (MQ=255)
                                                                                                                                                                               cgcACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGGCAGCCCCTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCCTTGGGCAGCATTAAGCCGATGCGGTCGCCCTCAACACTGTAttt                                                             >  2:826387/1‑239 (MQ=255)
                                                                                                                                                                                          aatttGCCCGGGCAAGTGCCAGGTTTTGCCTTTACGGAAAAA‑CGGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAACACTGTATTTTTCAAGGAt                                                    <  2:883833/234‑1 (MQ=255)
                                                                                                                                                                                              tGCTCCGGCAGATGCCAGAGATGGACTTTATCGAGAAA‑CTGGCGGGAAGGGAAGGTGGGTTTTATTTCA‑GCCGCCCTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCTTTGCGGGCCTGCGGCGC‑CGG‑GCGATGGCCCCAAAAATCACTGCCCCACCGTTCCCTCGCTTGGGACGCATTAACCCGTTGCCTTCCCCTTCAACCCTCGATTTTTTAACGAAGCGGcc                                              >  1:724154/1‑239 (MQ=255)
                                                                                                                                                                                                    ggaggtggcaaaagTTTGCCTTTACCGAGAAC‑CTGGGGGAAAGGAAAGGTGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAACACTGTATTTTTCAAGGATGCGTCCaacaaa                                        <  2:654060/227‑1 (MQ=255)
                                                                                                                                                                                                          tGCCAAAGTTTGCCTTTA‑CCAGAAACCGGGGGGGAGGTGAAGTGGGTTTTGATTCCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAACACTGTATTTTTCAAGGATGCGTCCAACAAACAGCg                                   <  2:520754/238‑1 (MQ=255)
                                                                                                                                                                                                                     gCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGCTGCCTGCATTGGGCAGCATTAAACCGATGGGTTCGCCTTCCAAACTTTATTTTTTAAAGAAGGGTCCAACAAACAGCGGTTTCGTAAACa                       >  1:793086/1‑239 (MQ=255)
                                                                                                                                                                                                                                      tGGGGGGAAGTGAATGTGGTTTGGTTTTCA‑GCCGCCTTAATAGCCCTGGTCCGCCCCTTTACCCCGGCAGTTTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCTATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAACACTGTATTTTTCAAGGATGCGTCCAACAAACAGCGTTTTCGTAAGCAATTTGCGATAGGAGTc       <  1:739493/239‑1 (MQ=255)
                                                                                                                                                                                                                                            gAAGTGAAGGT‑GTTTTTATTTCACGCCGCGGTAATAGCCCGGGGCAGCCCTTTTACCCCGGCAGTTTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAACACTGTATTTTTCAAGGATGCGTCCAACAAACAGCGTTTTCGTAAGCAATTTGCGATAGGAGTCTgggg  <  1:710343/238‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
TATGGACGACGCCTTTCGGATGGCCTTCAGAACCGGAGGTAAAAA‑GGATCAGCGCCTCTTCTTCCGGCTGCTGTTT‑AACCTGTG‑CCAGACGCGGCATCAGCAAATGAGCGAAGATCCATACTTTGTCGGCAGTGGT‑GACATCTGCTTTTAAAT‑CTTCCAGATAGACCCAGCGCACCTGAGTAAGTTGCTCCGGCAGATGCCAGAGTTTGCCTTTATCGAGAAA‑CTGGCGGGAAGTGAAGATGGTTTTGATTTCA‑GCCGCCGTAATAGCACTGGTCAGCCCTTTTACCCCGGCAGTGTAGTTCATCATTGCGGGCATGCGGCGA‑CGG‑GCGATGGCCCCAAAAATCACTGCCGCACTGATGCCTGCATTGGGCAGCATTAAGCCGATGCGTTCGCCTTCAACACTGTATTTTTCAAGGATGCGTCCAACAAACAGCGTTTTCGTAAGCAATTTGCGATAGGAGTCTGGGG  >  NC_000913/2974857‑2975324

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