Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,051 C→G 5.8% intergenic (+34/‑164) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0510CG 5.8% 209.9 / 6.0 69intergenic (+34/‑164)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base G (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

ACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCC  >  NC_000913/3780904‑3781198
                                                                                                                                                   |                                                                                                                                                   
aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTc                                                                                                                                                     <  5:50943/148‑1 (MQ=255)
  gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCTCTCCCTCAgg                                                                                                                                                  >  6:290235/1‑148 (MQ=255)
       cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTGGCCGGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTGACCCTAACCCTCTCCCTCAGGgaga                                                                                                                                              <  6:397226/148‑1 (MQ=255)
         ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                                                                            >  3:192173/1‑148 (MQ=255)
          gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGagag                                                                                                                                             <  5:185087/146‑1 (MQ=255)
          gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                                                                            <  8:56532/147‑1 (MQ=255)
            gcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                                                         <  4:192173/148‑1 (MQ=255)
                     aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                                                >  6:322652/1‑148 (MQ=255)
                          cGAAATTACGCAAGATTCGCTGGTGCCGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                                                                          <  5:290235/148‑1 (MQ=255)
                           gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggtt                                                                                                                                         <  6:234985/133‑3 (MQ=255)
                           gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggtt                                                                                                                                         >  5:234985/1‑131 (MQ=255)
                           gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                                                                          <  5:1385/148‑1 (MQ=255)
                            aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCATTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtat                                                                                                                         <  8:7228/148‑1 (MQ=255)
                              aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCTAAACCACGCCCTCAGGGAGAGGGGACCCGTCGGCGCagtatg                                                                                                                        >  8:277194/1‑147 (MQ=255)
                                tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgt                                                                                                                           >  3:91451/1‑142 (MQ=255)
                                tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgt                                                                                                                           <  4:91451/142‑1 (MQ=255)
                                   gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTcc                                                                                                                  <  2:313943/148‑1 (MQ=255)
                                   gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTcc                                                                                                                  <  8:366574/148‑1 (MQ=255)
                                    cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTc                                                                                                                                                     >  6:256930/1‑112 (MQ=255)
                                    cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTc                                                                                                                                                     <  5:256930/112‑1 (MQ=255)
                                    cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                                 >  3:115578/1‑148 (MQ=255)
                                      aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTc                                                                                                               <  7:109113/148‑1 (MQ=255)
                                      aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCCTTCGGCGCTGTATGTACTCCCt                                                                                                                >  4:23224/1‑147 (MQ=255)
                                      aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCTAACCCACTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTacca                                                                                                                >  2:444981/1‑146 (MQ=255)
                                                 tGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACg                                                                                                    <  7:277194/148‑1 (MQ=255)
                                                 tGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCACATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACg                                                                                                    <  7:66016/148‑1 (MQ=255)
                                                           gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                                              <  5:55474/96‑1 (MQ=255)
                                                           gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                                              >  6:55474/1‑96 (MQ=255)
                                                               aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCACATGGAGAGGGGACCGTTTGGCGCTGTATGTACTCCCTCACTCTGAAACGCAACCGCCCACttt                                                                                      >  1:411391/1‑148 (MQ=255)
                                                                aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttt                                                                                     <  5:322652/148‑1 (MQ=255)
                                                                     ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCATTTTTtct                                                                                >  2:171201/1‑148 (MQ=255)
                                                                       gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                               <  4:115578/147‑1 (MQ=255)
                                                                          tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCACAGGGAGAGGGGACCGTTCGGCGCTGTATGTTCTCCCGCACTCTGAAACGCCCCCGCACTCTTTTTTTCTCCCTc                                                                           >  7:137620/1‑148 (MQ=255)
                                                                                   ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCGTTTTTTCTCCCTCGCCCCTCCg                                                                  >  3:73885/1‑148 (MQ=255)
                                                                                        ccATGGCGAAAGGTAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGATGGGACCGTTCGGCGCTTTATTTACTCCTGCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                              <  6:306794/147‑1 (MQ=255)
                                                                                         cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCACCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTAATCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                              >  4:115660/1‑146 (MQ=255)
                                                                                         cATGGCGAAAGGGAAGGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                            <  1:361345/148‑1 (MQ=255)
                                                                                                 aaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                    >  1:113356/1‑148 (MQ=255)
                                                                                                    ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                 <  2:113356/148‑1 (MQ=255)
                                                                                                     gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAACCACTACCTCAGGGAGAGGGGACCGTTAGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACAATTTTTTTATCCCAAGCCCCTCCGGGGAGAGGGCCggggg                                                  >  6:273366/1‑145 (MQ=255)
                                                                                                            gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                 <  4:55910/124‑1 (MQ=255)
                                                                                                            gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                 >  3:55910/1‑124 (MQ=255)
                                                                                                            gCATAGGGGTGTGCTCCCCCTCACCCTACCCCTCTCCCTCATGGAGAGGGGACCGTGCGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTTAGGGGaaa                                         <  2:411391/148‑1 (MQ=255)
                                                                                                               tAGTCGTTTGGGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaa                                                                                                       <  3:45983/83‑1 (MQ=255)
                                                                                                               tAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                        <  3:235460/130‑1 (MQ=255)
                                                                                                               tAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                 <  5:357432/137‑1 (MQ=255)
                                                                                                                aGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCGGTATGTACTACCTCACTCTGAACCGACACCGCCCTCTTTTTGTCTACCGCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                     >  2:289323/1‑148 (MQ=255)
                                                                                                                aGTCGTTTGCGCCCCCTCACCCTAACCCACTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTACCTCACTCTGAAACGACACCGCACTCTTTTTTTATCCCACGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                     >  8:348317/1‑148 (MQ=255)
                                                                                                                aGTCGTTTGCGCCCCCTCACCCTAACCCACACCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTAATCCCTCACTCTGAAACGCCACCGCACGCTTTTTTTCACCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGa                                                >  6:357432/1‑137 (MQ=255)
                                                                                                                 ggcGTTTGCGCCCCCTCACCCAAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGAACCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                     <  7:394973/145‑1 (MQ=255)
                                                                                                                 gTCGTTGGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTGCGGCGCTGTATGTCCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTGCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGcc                                    <  3:23224/148‑1 (MQ=255)
                                                                                                                  tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                      >  3:328139/1‑129 (MQ=255)
                                                                                                                  tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                      <  4:328139/129‑1 (MQ=255)
                                                                                                                  tCGTTGGCGCCCCCACACCCAAACCCTCTCCCTCCGGGCGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGTGAGGGCCGGGGTGAGGGGAAAAGGcc                                    <  1:444981/147‑1 (MQ=255)
                                                                                                                           cccccTCACCCTAACCCTCTCCCTCAGGGAGCGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACCACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCa                          <  3:212082/148‑1 (MQ=255)
                                                                                                                              ccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAcc                        >  6:323840/1‑147 (MQ=255)
                                                                                                                                tCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                     >  3:307198/1‑148 (MQ=255)
                                                                                                                                   acggtccccTCTCCCTGAGGGAGAGGGTACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                  >  7:77903/7‑148 (MQ=255)
                                                                                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                             <  7:130117/104‑1 (MQ=255)
                                                                                                                                    ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                             >  8:130117/1‑104 (MQ=255)
                                                                                                                                     cTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                >  7:177310/1‑148 (MQ=255)
                                                                                                                                      aaaCCCTCTCCCGCGGGGAGAGGGGACCGTTCGGCTCTGTATGTAGCCACTCACTCTGAAACCACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGTGAGGGCGGGGGTTAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                <  2:246999/146‑1 (MQ=255)
                                                                                                                                       aaCCCTCTCCCTCGGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg               >  1:246999/1‑147 (MQ=255)
                                                                                                                                       aaCCCTCTCCCTCCGGGAGTGGGGACCGTTCGGCGCTTTATGTACTCCCTCACTCTTAACCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc              <  4:307198/148‑1 (MQ=255)
                                                                                                                                       aaCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                      <  2:42785/124‑1 (MQ=255)
                                                                                                                                       aaCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                      >  1:42785/1‑124 (MQ=255)
                                                                                                                                        ccccTCACCCTCCGCGCGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTCAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCGGGGGAGTGGGCGGGGGTGAGGGTAAAAGGCCGCACTGCTCACCAATTTTCTg               <  7:410054/145‑1 (MQ=255)
                                                                                                                                           ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc          >  6:237954/1‑148 (MQ=255)
                                                                                                                                            tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                     <  7:453421/104‑1 (MQ=255)
                                                                                                                                            tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                     >  8:453421/1‑104 (MQ=255)
                                                                                                                                            tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc         >  3:354406/1‑148 (MQ=255)
                                                                                                                                               cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc      >  3:84489/1‑148 (MQ=255)
                                                                                                                                                   cAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc  >  1:382619/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCC  >  NC_000913/3780904‑3781198

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