Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,714,488 A→C L636W (TTG→TGG)  bisC ← biotin sulfoxide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,714,4880AC97.9% 154.6 / ‑3.8 49L636W (TTG→TGG) bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (25/22);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

AATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGA‑CCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCAACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGTCCGGGTTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAAC  >  NC_000913/3714259‑3714718
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       
aaggCAAGTCCCGCCAGGTTAACTTCCCT‑GCCAAAGACCCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACACCCGACCAGTATCGCCATCTTGTATGCCGCGCCCCTGGGAATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCCCCGCGTCCAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTCCCTGCCACTGtt                                                                                                                                                                                                                                   <  2:197854/233‑1 (MQ=255)
aaTGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGA‑CCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGtt                                                                                                                                                                                                                                   <  1:518190/236‑1 (MQ=255)
aaTGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGA‑CCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGACGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCCGACA‑GTACCTGCGAc                                                                                                                                                                                                                                       >  2:449589/1‑232 (MQ=255)
                            tCGCTAATGA‑CCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGACGGGCAGAAC‑GGACCTGCCCCTGGTCTGGGTCGGCATTTGCCTGGCAATCGt                                                                                                                                                                                                        >  2:264143/1‑235 (MQ=255)
                                               gCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCGGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTgg                                                                                                                                                                                     <  2:181943/236‑1 (MQ=255)
                                                                          aaCCGAAGTGTATCGTCATCTTGTATCCCGCGCTCCTGGGCATCGTCAGGATGAAGGGTGACAGGCTCACGTTTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTa                                                                                                                                                          <  1:86235/236‑1 (MQ=255)
                                                                            ccGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGGCCCAGCCACATTGGATGCCCACGGCAATCCCGGTGACCGTaa                                                                                                                                                         >  2:86235/1‑235 (MQ=255)
                                                                               aaCAGTATCGCCATCTTGTATGCCGCGCTCCTGGTCATCGTAAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCATGGCAATCCGGGTAACCGTACTCgg                                                                                                                                                     >  1:8715/1‑236 (MQ=255)
                                                                                            tCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGtt                                                                                                                                                                                                                                   >  2:660967/1‑145 (MQ=255)
                                                                                            tCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGtt                                                                                                                                                                                                                                   <  1:660967/145‑1 (MQ=255)
                                                                                            tCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGACCGGTGCCAGCCACATTGGATGCCCAGGGCAATCCGGGGAAACGTAATCGGCAATACGCtgtgg                                                                                                                                        >  2:484804/1‑235 (MQ=255)
                                                                                                      cgcgCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATAc                                                                                                                                               <  1:82568/219‑1 (MQ=255)
                                                                                                      cgcgCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATAc                                                                                                                                               >  2:82568/1‑219 (MQ=255)
                                                                                                        cgcTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGGATTGCCCTGCCATTCGTCTGGTTCCAGCCAAATTGGATGTCCAGGGCAATCTGGGTAAACGTAATCGGGAATATGATGTGATAAGATTTTAAt                                                                                                                            >  2:155477/1‑236 (MQ=255)
                                                                                                        cgcTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCg                                                                                                                                                                                                     <  1:140740/163‑1 (MQ=255)
                                                                                                        cgcTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCg                                                                                                                                                                                                     >  2:140740/1‑163 (MQ=255)
                                                                                                        cgcTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAAt                                                                                                                            >  1:514624/1‑236 (MQ=255)
                                                                                                                     cGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCCATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCg                                                                                                               >  1:172832/1‑236 (MQ=255)
                                                                                                                     cGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACAT‑‑‑‑TGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGgtt                                                                                                           >  1:614996/1‑236 (MQ=255)
                                                                                                                               aaTGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCcac                                                                                                                                                                                          <  1:560309/151‑1 (MQ=255)
                                                                                                                               aaTGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCcac                                                                                                                                                                                          >  2:560309/1‑151 (MQ=255)
                                                                                                                               aaTGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCCATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGGTTTTAACg                                                                                                     >  1:129417/1‑236 (MQ=255)
                                                                                                                                 tGGTGACAGGCTCACGATTGGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGa                                                                                                   <  1:449589/236‑1 (MQ=255)
                                                                                                                                 gggTGGCAGGCTAACGTTTTGCCACCGCGTACAATTCGCGAAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGa                                                                                                   <  1:477976/235‑1 (MQ=255)
                                                                                                                                     gACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGGTTTTAACGGATtcg                                                                                               >  1:218507/1‑236 (MQ=255)
                                                                                                                                     gACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGGTTTTAACGGATgcg                                                                                               >  1:178569/1‑236 (MQ=255)
                                                                                                                                     gACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCAACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGGAACCGTAATCGGCAATACGCTGTGAGACAAGTTTAATCCTGCCCCTGGCGGGTTTTAAAgggtatg                                                                                               >  2:507736/1‑231 (MQ=255)
                                                                                                                                      aCAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGc                                                                                                                                                                          >  1:29795/1‑160 (MQ=255)
                                                                                                                                              aCGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGGTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAAGCTTGCCGCTGGCGGATTTTAACGGGTGCGCCAGCCGAt                                                                                      >  1:312352/1‑236 (MQ=255)
                                                                                                                                                    tGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGGTTTTTAAGGATGCGCCAGCGGATCGCGGc                                                                                >  1:428660/1‑236 (MQ=255)
                                                                                                                                                            cgTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAAt                                                                                                                                                                       >  1:674025/1‑141 (MQ=255)
                                                                                                                                                            cgTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAAt                                                                                                                                                                       <  2:674025/141‑1 (MQ=255)
                                                                                                                                                                aaaTCCGCGCAAAGAACTGTAATCCAGCTGGCTGTGCAGGCGGTCCGCCGGATGGGCAGAAA‑GTACCTGCCTCTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCagcga                                                                    <  1:125619/235‑1 (MQ=255)
                                                                                                                                                                            agaACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGGTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATAAGCGAAGC‑GAATaaaa                                                        >  1:678323/1‑234 (MQ=255)
                                                                                                                                                                                                tgtgCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGGTTTTAACGGATACGCCACCGGATCACGGCGCAAATCAGCGAAGCGGAAAAAACCTCCCCCGGCCCGGGtttt                                     >  2:193106/1‑235 (MQ=255)
                                                                                                                                                                                                  tgCAGGGGGTGCCCCGGATGGGAAGAAA‑GCACCTGCCACTGTGCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACAT‑‑‑‑TGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGTCCGGGTTTTCCGGCAt                              <  2:614996/235‑1 (MQ=255)
                                                                                                                                                                                                    ggggcggcGCCCCGAATGGACAAAAT‑GCCCCTGCCACTTTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGTCCGGGTTTTCCGGCa                               <  2:428660/234‑1 (MQ=255)
                                                                                                                                                                                                     cggCGGGGCCCCGGGAGGGCAAAAAGGTCCCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGTCCGGGTTTTCCGGCa                               <  2:678323/235‑1 (MQ=255)
                                                                                                                                                                                                      gggggTGCGCCGGGTGGGCGGAAG‑GTCCCTGCCTCGGTGTTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGTCCGGGTTTTCCGGCAt                              <  1:264143/232‑1 (MQ=255)
                                                                                                                                                                                                       gCGGTGCGCCGGATGCGCAAAAT‑GCACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGTCCGGGTTTTCCGGCAtc                             <  2:218507/235‑1 (MQ=255)
                                                                                                                                                                                                           tccgcggaaggggAAGAAA‑TTCCCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGTCCGGGTTTTCCGGCAtc                             <  1:645778/221‑1 (MQ=255)
                                                                                                                                                                                                            gcgcCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCg                                                                                                                     >  2:251021/1‑143 (MQ=255)
                                                                                                                                                                                                            gcgcCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCg                                                                                                                     <  1:251021/143‑1 (MQ=255)
                                                                                                                                                                                                            gcgcCGGATGGGCAGAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGGTTTTAACGGATGCGCCAGCGGATCGCGGGAAAAAACAGCGAAGC‑GAATAAACCGCTCGCTGGTCGGGGGTTCCGGCAtctcca                         >  1:544154/1‑232 (MQ=255)
                                                                                                                                                                                                                   ttGGGCAGAAA‑GTCCCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGTCCGGGTTTTCCGGCATCTCAATTAACTGGt                <  1:606517/235‑1 (MQ=255)
                                                                                                                                                                                                                          gAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCa                                                                               >  2:279354/1‑167 (MQ=255)
                                                                                                                                                                                                                          gAAA‑GTACCTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCa                                                                               <  1:279354/167‑1 (MQ=255)
                                                                                                                                                                                                                          acactgtcccTGCCATTGTTCTGGTTCGGAATGGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCACAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGTCCGGGTTTTCCGGCATCTCAATTAACTGGTTGGCt           <  2:129417/227‑1 (MQ=255)
                                                                                                                                                                                                                                   cTGCCACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCCAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGGCCGGGGTTTCCGGCAACTCCATTAACTGGTTGGCTTGCCCGAAc  >  2:187412/1‑235 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       
AATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGA‑CCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAA‑GTACCTGCAACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCGGCTCCAGCCACATTGGATGCCCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGC‑GAATAAACCGCTCGCTGTCCGGGTTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAAC  >  NC_000913/3714259‑3714718

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