Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,509 A→G A391A (GCA→GCG)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,5090AG98.2% 168.6 / ‑4.5 56A391A (GCA→GCG) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (26/28);  total (26/30)
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.62e-01

TAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGAT  >  NC_000913/1275‑1748
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
taataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGTTTACGCAATCATCTTCCGCATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGaacgaaccgccgg                                                                                                                                                                                                                                   >  1:209783/5‑239 (MQ=255)
   taTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTggac                                                                                                                                                                                                                                 >  1:400077/2‑248 (MQ=255)
    aTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGTATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTGCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAt                                                                                                                                                                                                                               <  1:305385/251‑1 (MQ=255)
         aaTGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTAtct                                                                                                                                                                                                                          >  1:178941/1‑251 (MQ=255)
          atggtCAGCGTGTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCAGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTAGCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCGGGCGGTGACGGTTCTCCTGGCCATTAtctc                                                                                                                                                                                                                         <  1:197310/247‑1 (MQ=255)
          aggTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTAtctc                                                                                                                                                                                                                         <  1:234003/249‑1 (MQ=255)
            gTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACCCGCCCGTCTTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGAGACTGTG‑TGCGAGCTGAACGGGCAAGGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGTCGGAAAGGCGGGCCACTTTCTCgg                                                                                                                                                                                                                       >  2:499872/1‑251 (MQ=255)
                      tccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTATCTCggggggaggggg                                                                                                                                                                                                             >  1:211457/7‑241 (MQ=255)
                       ccggcccGGGGATGAAAGGGATGGTTGGCATGGGGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGa                                                                                                                                                                                                            <  2:178941/245‑1 (MQ=255)
                           cctggggATGAAAGGGATGATTGGCACGGGGGCTCGCGTGTTTGCAGGGATTCCACGCGCCTGTATTTCCGTTGGT‑CTGATTTCGCAATCATCCTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGGG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGATCTGAAAGAAGTCTTACTGGAGCCGCTGGCGGGGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTa                                                                                                                                                                                                       <  2:363392/248‑1 (MQ=255)
                                ggATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCCTTATCCCGGGGG‑GAGGTGATGGGATgggg                                                                                                                                                                                                  >  1:114108/1‑248 (MQ=255)
                                       aGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGc                                                                                                                                                                                          <  1:263026/251‑1 (MQ=255)
                                            tGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGCTGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACCCAATCCTCTTCCGAATACCGTATCAGTTTCTGCGTTCCGCAAAGCGCCTGTG‑GGCGAGCTGAACGGGCAATGCTGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAAGGGAAGGCCGGCCATTATCTCGGGGG‑GAGGGGAGGGGTGGC‑ACACCTTGCGTggg                                                                                                                                                                                     >  2:155773/1‑251 (MQ=255)
                                                 ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTGCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAGTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGAtctc                                                                                                                                                                                <  1:23288/251‑1 (MQ=255)
                                                               ggcttttttaGCGATGTCACGCGTCCGTATTTCCGG‑GGGGCTGATTACGCAATCCTCGTTCTACTTCAGTATGAGTTTCTGCGTTCCGCAAAGCGACTGGG‑TGCGAGCTGAACGGGCAATGCGGTAAGAGTGCTACCTGTAACTGAAAGAAGGCTTATTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCCCGGCGAAACTCTTTg                                                                                                                                                                  <  2:107071/242‑1 (MQ=255)
                                                                      gCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACt                                                                                                                                                           <  1:183381/251‑1 (MQ=255)
                                                                          cGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCCTTATCTCGGGGG‑TAGGTGATGGGATGC‑GCACCTTGCGGGGGATCTCGGCGGAATTCTTTTCCGcgcctgcc                                                                                                                                                       >  2:371728/1‑244 (MQ=255)
                                                                           gATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAAGGAACGGCTGGCCCTTATCTCGGTGG‑TAGGTGCTGGTATTC‑GCACCATGCTTGGGATTTCGGCGGAATTTTTTTCCGcggctgccc                                                                                                                                                      >  2:265921/1‑243 (MQ=255)
                                                                           gATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCGGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGGTCTACCTGGAACTGAAAGAAGGATTGCTGGAGCCGCTGGAGGTGTCGGGACGGCTGGCCATTATATCGGTGGGTGGGGGAGGGGTTGG‑ACACCT‑GCGTGGGATCTCGGGGAAATTATTTGCCGcggatgccc                                                                                                                                                      >  2:242115/1‑243 (MQ=255)
                                                                                cACGCGCCCGTATTTCCGT‑GGGGCTGGTTACGCATTCTTCTTCCGATGTCAGGATCAGTTTCTGCGGTCCGCAAATCGACTGTG‑TGCGAGCTGAACGGGGAATGCAGGAAGAGTTCTACCGGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGGCATTATCGCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCTCGCTGGCCCGCGcc                                                                                                                                                 <  2:319268/251‑1 (MQ=255)
                                                                                     gcCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGGGATGGGATGG‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAatat                                                                                                                                            >  1:283040/1‑251 (MQ=255)
                                                                                      cccGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGGGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCtt                                                                                                                                                                    >  1:39199/1‑226 (MQ=255)
                                                                                      cccGTATTCCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTCCCTGGAACTGAAAGAAGGCGTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCtt                                                                                                                                                                    <  2:39199/226‑1 (MQ=255)
                                                                                        cGTAGTTCCGA‑GGTGCTGATTTCTCAATCATCTCCCGAATACAGTAACAGTTTCTGCGGTCCGCAAAGCTACCGTG‑TGCGAGCTGAGCGGGCAAGGCAGGGAGAGGTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCGCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCaa                                                                                                                                         <  2:5130/251‑1 (MQ=255)
                                                                                                       gCTGATTACGCATTCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGTGATCTCGGCGAAATTTTTTGCCGCGCTGGGCCGCGTCAATATCAACATTGTCGCCAtt                                                                                                                            <  1:189230/250‑1 (MQ=255)
                                                                                                       gCTGATTACGCAATCAACTTCCGAAGGCAGTATCAGTTTCTGCGTTCCGCAAGGCGAC‑GTGGTGCGAGCGGAACGGGCAATGCAGGAAGAGTTCTCCCTGGATCTGAAAGAAGGCTCACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCGTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTg                                                                                                                           <  2:50197/251‑1 (MQ=255)
                                                                                                          gATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAACGCGACTGTG‑TGCGAGCTGAACGGGCACTGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAAGGAACGGCTGGCCATTATCTCGGGGG‑GAGGGGATGGTTAAC‑GCACCATGCGGGGGAACACGGCGAAATTCTTTGCCGGGCTGGCCCGCGCCAATATCCACAATGTACCCATTGCTc                                                                                                                        >  2:224212/1‑251 (MQ=255)
                                                                                                           aTTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCACCCTTGTCGCCATTtctcc                                                                                                                       >  1:216865/1‑250 (MQ=255)
                                                                                                                      atcTTCCGCATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGAGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGGATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTTCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAtcttatg                                                                                                            >  1:232817/1‑248 (MQ=255)
                                                                                                                                 aCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGTAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTg                                                                                                                                                                                           <  2:490530/161‑1 (MQ=255)
                                                                                                                                 aCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTg                                                                                                                                                                                           >  1:490530/1‑161 (MQ=255)
                                                                                                                                        caggttCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCg                                                                                          <  1:129368/247‑1 (MQ=255)
                                                                                                                                          gTTTCTGCGTTCCGCAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCCGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCAATATCTCGGTGG‑TGGGGGGTGGTATGC‑GAAACTTGCGTGGGTTCTCGGCGAAATTCTTTGCACCGCATGCCCGCGCCACAATCAACATAGTTGCCCTTTCTACGGGGTCTTTTGCACGCCCAATCTCTGTtgtg                                                                                        >  2:145162/1‑251 (MQ=255)
                                                                                                                                                    cccgcaAAGCGGCTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTAGCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCACTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACg                                                                              <  1:169415/247‑1 (MQ=255)
                                                                                                                                                               aCTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGAcc                                                                   <  1:481842/251‑1 (MQ=255)
                                                                                                                                                                      tgCGAGCTGAACGGGCAATGCAGGACGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAtct                                                                                                                >  1:136147/1‑200 (MQ=255)
                                                                                                                                                                      tgCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAtct                                                                                                                <  2:136147/200‑1 (MQ=255)
                                                                                                                                                                         gAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTTAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCCACCACTGGCGTg                                                          >  1:303873/1‑251 (MQ=255)
                                                                                                                                                                           gCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTgcg                                                        <  1:170842/251‑1 (MQ=255)
                                                                                                                                                                              gAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTg                                                                                                            <  2:92941/196‑1 (MQ=255)
                                                                                                                                                                              gAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTg                                                                                                            >  1:92941/1‑196 (MQ=255)
                                                                                                                                                                               aaCGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGtt                                                    <  1:70311/251‑1 (MQ=255)
                                                                                                                                                                               aaCGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACGGGGGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGtt                                                    <  2:211457/251‑1 (MQ=255)
                                                                                                                                                                                aCGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTa                                                   <  1:258914/251‑1 (MQ=255)
                                                                                                                                                                                     cAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACCGAACGGCCGGCCCTTATCTCGGTGG‑TAGGGGATGGTATGCTACCCCTT‑CGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCCCCAATATCAAAATTGTCCCCATTGCTCAGGGGCCTTCTGGACGCTCAATCTCCTTCGCGGGAAATAACGATGATGCACCCACTGGCGTGCGGGGTACtcac                                              >  2:104476/1‑250 (MQ=255)
                                                                                                                                                                                                gagTTCTACCTGGAACTGAAAGAAGTCTTACTGTAGCCGCTGGCTGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGTGATCTCGGCGAAATTCTTTGCCGCGCTGGTCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCGGGCGTGGTAAATAACGATGATGCGTCCACTGGCGTGCGCGTTACTCATCAGATGCTGtt                                   <  1:163130/251‑1 (MQ=255)
                                                                                                                                                                                                          tGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACCGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCgtt                         >  1:429785/1‑249 (MQ=255)
                                                                                                                                                                                                           ggAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGAt                         <  1:443914/250‑1 (MQ=255)
                                                                                                                                                                                                                  aaaGAAGGCTTACTGGAGCGTCTGGTGGTTACGTGACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGTGAAAGTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGACGATGCGACCGCTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTAt                 <  2:327912/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaaGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGtt                                   <  2:404845/233‑1 (MQ=255)
                                                                                                                                                                                                                  aaaGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGGGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGtt                                   >  1:404845/1‑233 (MQ=255)
                                                                                                                                                                                                                               tGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGGTGGTATGC‑GCACCTTGCGGGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATCACGATCAGGTTATCGCAGTGTTtgtg    >  2:6229/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTTCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGGGGGCGCGGTTCTCCTCAGGTGCTTTTTAATACCGATTAGGTTATCGAAGGGTTtgtgg   >  2:144133/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gggcctcTGGCGGTGCCGGACCGGCGGGCCAGTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTTGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGTCGTGCGAGTTACTCATCAGATGCTGTGCAATACCGAGCAGGTTATCGAAGTGTTTGTGAt  <  2:176575/245‑1 (MQ=255)
                                                                                                                                                                                                                                   gCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGaa                                                                                                          >  1:258610/1‑145 (MQ=255)
                                                                                                                                                                                                                                   gCCGCGGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCGGGCCCGCGCCAATATCAACATTGTCGCCAGTGCTCAGGGATCTGCTGaa                                                                                                          <  2:258610/145‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
TAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGT‑GGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTG‑TGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGG‑TAGGTGATGGTATGC‑GCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGAT  >  NC_000913/1275‑1748

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