Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,051 C→G 6.2% 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 6.2% 235.3 / 5.8 78intergenic (+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 (36/36);  new base G (3/2);  total (40/38)
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.33e-01

CGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCC  >  NC_000913/3780905‑3781198
                                                                                                                                                  |                                                                                                                                                    
cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                                                    <  7:78411/149‑1 (MQ=255)
 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGTAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCTAACCCACACCCGCAgg                                                                                                                                                   >  8:149075/1‑149 (MQ=255)
 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAgg                                                                                                                                                   <  7:37466/149‑1 (MQ=255)
 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                                   <  8:204516/149‑1 (MQ=255)
 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGAGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                                   >  4:301904/1‑149 (MQ=255)
   tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAgg                                                                                                                                                   >  8:179719/1‑147 (MQ=255)
   tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                                   >  4:239333/1‑147 (MQ=255)
   tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                                   >  7:261621/1‑147 (MQ=255)
   tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCACTCCCTCAgg                                                                                                                                                   >  7:267657/1‑147 (MQ=255)
         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAggg                                                                                                                                            <  6:243047/148‑1 (MQ=255)
         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                                                            <  6:292581/148‑1 (MQ=255)
         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                                                            <  3:160171/148‑1 (MQ=255)
                aTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTc                                                                                                                                   <  4:75166/149‑1 (MQ=255)
                    aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACACTATACCTCAGGGAGAGGGGACCGGTCGgcg                                                                                                                               >  7:168993/1‑149 (MQ=255)
                       aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                                                                            >  1:19418/1‑149 (MQ=255)
                       aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                             <  2:77027/148‑1 (MQ=255)
                           aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatg                                                                                                                        <  5:170635/149‑1 (MQ=255)
                           aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatg                                                                                                                        <  5:290449/149‑1 (MQ=255)
                           aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCGGCGGCACTGGCTCCCATGGCGAAATGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCAATCCCACAAGGATAGGGTCACCGTCGGCGCtggaat                                                                                                                        >  6:240784/1‑145 (MQ=255)
                            aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgt                                                                                                                       <  2:19418/149‑1 (MQ=255)
                                 cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTc                                                                                                                   <  4:288058/148‑1 (MQ=255)
                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGTCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCTGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                                 <  2:351616/149‑1 (MQ=255)
                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGAc                                                                                                                                        >  2:274162/1‑126 (MQ=255)
                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGAc                                                                                                                                        <  1:274162/126‑1 (MQ=255)
                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGAGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                                 >  3:86812/1‑149 (MQ=255)
                                     aGATGCGCGGGTGCAGGGGCTGGGAACAGCGTCGCCTGCGGCACTGGCTCCCACGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTACCCCTCTCCCTCAGCGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                                              <  4:205509/149‑1 (MQ=255)
                                       aTTCGCTGGTGCAGGGTCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                                            <  6:102134/149‑1 (MQ=255)
                                          cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActct                                                                                                          <  8:261621/148‑1 (MQ=255)
                                          cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActct                                                                                                          <  8:267657/148‑1 (MQ=255)
                                           gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTAAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCACTg                                                                                                         >  6:193563/1‑148 (MQ=255)
                                                 gCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                                 <  5:155433/134‑1 (MQ=255)
                                                 gCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                                 >  6:155433/1‑134 (MQ=255)
                                                   aGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGTGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGacac                                                                                                <  5:240784/149‑1 (MQ=255)
                                                   aGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGacac                                                                                                <  7:76328/149‑1 (MQ=255)
                                                       gCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                                                    >  7:326204/1‑94 (MQ=255)
                                                       gCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                                                    <  8:326204/94‑1 (MQ=255)
                                                          gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCACTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCActc                                                                                         >  2:378295/1‑149 (MQ=255)
                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGGCCGTTCGGCGCTGTATGTACTCCCTCACGCTGAAACGCCACAGCACTCatt                                                                                      >  6:15863/1‑149 (MQ=255)
                                                             tAAAGAGTTGCCTGAGGCACTGTCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTATCCCTCAGGGAGATGGGACAGTTCGGCGGAGTATGTACTCCCTCACACAGAAACGACACAGccatcatt                                                                                      >  4:34997/1‑142 (MQ=255)
                                                               aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                                    <  7:149075/149‑1 (MQ=255)
                                                                   gTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                                                                                        >  5:134898/1‑125 (MQ=255)
                                                                   gTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                                                                                        <  6:134898/125‑1 (MQ=255)
                                                                    ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTtctc                                                                               >  4:49904/1‑149 (MQ=255)
                                                                       ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTccc                                                                             <  3:324328/148‑1 (MQ=255)
                                                                           cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                                                                        >  3:110584/1‑149 (MQ=255)
                                                                           cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTCCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                                                                          <  4:149460/147‑1 (MQ=255)
                                                                             gCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTAGTTTGCGCCCCATCCCCATCACCCTCGCCCACAGGGAGACGGGCCCGGTCGGAGCGGGATGTTATACCGCACTCGGACCCGACACCGCCCGCTTTTGTTCTCCCTCGcccc                                                                      >  6:254375/1‑149 (MQ=255)
                                                                              cACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCACCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACGCTTTTTTTCTCCCTCGCCCCt                                                                     >  6:85531/1‑149 (MQ=255)
                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTACCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCg                                                                  >  2:36688/1‑149 (MQ=255)
                                                                                   gCTCCCATGGCGACAGGGAATGCGGCATAGTCGGTGGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCGGTATGTACTCCCGCACTCTGAACCGACCCCGCACGCTTTTTTTCTCCCTCGCCCCTCCggg                                                                <  6:135/149‑1 (MQ=255)
                                                                                   gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                <  4:110584/149‑1 (MQ=255)
                                                                                        cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCt                                                                            >  6:178813/1‑132 (MQ=255)
                                                                                        cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCt                                                                            <  5:178813/132‑1 (MQ=255)
                                                                                              gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCt                                                                                                                >  2:96815/1‑90 (MQ=255)
                                                                                              gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCt                                                                                                                <  1:96815/90‑1 (MQ=255)
                                                                                               aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCGAACCCACACCCGCAGGGAGAGGGGACCGGTCGGCGCGGTATGTTAGCCCAGCCTCTGAAACGACCCCGCGCTCTTTTTTTATCCCTCCCCCCGCCGGGGGGAGGGCCcgg                                                    >  1:194173/1‑149 (MQ=255)
                                                                                                    gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACAGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGTAGAGGGCCGGGGTGAg                                               <  3:34997/149‑1 (MQ=255)
                                                                                                     aaggCGGCCTAGGCGGTTGCCCCCCCGCACCCTCACCCGCTCCCTCAGGGATAGGGGACAGTTCGCCGCTGTATGTCATCCCTCACTCTGAACCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                              <  5:254375/146‑1 (MQ=255)
                                                                                                     aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                              <  5:15863/149‑1 (MQ=255)
                                                                                                         cGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCACCTCCGGGGAGAGGGCCGGGGTGAGGGGa                                           >  6:332646/1‑148 (MQ=255)
                                                                                                           gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                        <  5:193563/149‑1 (MQ=255)
                                                                                                              tAGTCGTTTGCGCCCCCTCACCCAAACCCTCACCCTCAGGGAGAGGGGACCGGTCGGCGCTGTATGTACACCCTCACTCAGAAACGCCACCCCACAAGTTTTTTCTACCACCCCCCACCGGGGAGAGGGCCGGGGGGCGGGGAAAAGGc                                     >  8:179301/1‑149 (MQ=255)
                                                                                                                 tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCg                                   >  1:138653/1‑148 (MQ=255)
                                                                                                                 tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCGCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACACAGAAACGCCACCGCACTCATTTTTTCTCCCACCCCCCACCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGc                                  >  8:176426/1‑149 (MQ=255)
                                                                                                                            cccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                               <  7:62514/125‑1 (MQ=255)
                                                                                                                             ccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCaa                      >  2:92624/1‑149 (MQ=255)
                                                                                                                               tCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGGACTGCTCACCAAtt                    >  5:78281/1‑149 (MQ=255)
                                                                                                                                    cTACCCCTCTCCCTCTGGGAGTGGGGACCGTTAGGAGCTGTATGTAGTCCCTCACTCTCAACCGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGTGAAAAGGCCGCACTGCTCACCAATTTTCTg               <  4:86812/149‑1 (MQ=255)
                                                                                                                                     tAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc              >  3:290480/1‑149 (MQ=255)
                                                                                                                                      aaCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                             <  7:35987/117‑1 (MQ=255)
                                                                                                                                      aaCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                             >  8:35987/1‑117 (MQ=255)
                                                                                                                                        cccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTAt           >  1:67530/1‑149 (MQ=255)
                                                                                                                                         ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCACCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc          >  5:338943/1‑149 (MQ=255)
                                                                                                                                          ctctCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc         >  5:17021/1‑149 (MQ=255)
                                                                                                                                            ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg       >  7:64606/1‑149 (MQ=255)
                                                                                                                                            ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCACCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg       <  6:338943/149‑1 (MQ=255)
                                                                                                                                                cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc   >  8:108600/1‑149 (MQ=255)
                                                                                                                                                 tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTGTTTTTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGTGAGGGGAAAAGGCCGCTCTGCTCACCAATTTTCTGCTATCCTGccccc  >  6:203037/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCC  >  NC_000913/3780905‑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: