Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 81,044 C→T L305L (CTG→CTA)  leuB ← 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091381,0440CT97.9% 125.2 / ‑6.2 47L305L (CTG→CTA) leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/25);  minor base A (0/1);  total (21/26)
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.00e+00

GGGTTTCGTTTTCGGCTTCGTACACAACGTGAGCGTCGAACAATTTTTCGTATAACGTCTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCGAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGA  >  NC_000913/80798‑81282
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
gggTTTCGTTTTCGGCTTCGTACACAACGTGAGCGTCGAACAATTTTTCGTATAACGTCTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCg                                                                                                                                                                                                                                             <  1:269985/251‑1 (MQ=255)
      cGTTTTCGGCTTCGTACACATCGTGAGCGTCGAACAATTTTTCGTATAACTACTGAGCCATGATTATACCCCTGCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGTGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAgg                                                                                                                                                                                                                                        <  2:13096/250‑1 (MQ=255)
              gCTTCGTACACAACGTGAGCGTCGAACAATTTTTCGTATAACGTCTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCAATTTCATCGGTACTAACGGCGGAAGCGCCACGGGCTAAATCCCCGGTTCG‑AATGCCTTCTTCTAATGCGCGGATAATGGAGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCATGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGc                                                                                                                                                                                                                                >  1:131677/1‑250 (MQ=255)
                   gTACACAACGTGAGCGTCGAACAATTTTTCGTATAACGTCTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGGGCG‑AATGCCTTCTTCTAATGCGCGGATAATGGCGCGATCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCgg                                                                                                                                                                                                                          >  2:75431/1‑251 (MQ=255)
                   gTACACAACGGTAGGGTCGAACAATGTTTCGTATAACGTCTTAGCCATGATTACACCCCTTCGGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCTACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGACGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTTCAATCgg                                                                                                                                                                                                                          <  2:241491/251‑1 (MQ=255)
                           cGTGAGCGTCGAACAATTTTTCGTATAACGTCTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGAGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGa                                                                                                                                                                                                                                       >  1:268139/1‑230 (MQ=255)
                           cGTGAGCGGCGAACAAGTTTTCGTATAACGTCTTAGCCATGATTACACCCCTTCTGCGACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGa                                                                                                                                                                                                                                       <  2:268139/230‑1 (MQ=255)
                               aGCGTCGAACAATTTTTCGTATAACGTCTTAGCCAGGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGtt                                                                                                                                                                                                              <  2:143136/251‑1 (MQ=255)
                                    cGAACAATTTTTCGTATAACGTCTTAGCCATGATTACCCCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTg                                                                                                                                                                                                          <  2:256614/250‑1 (MQ=255)
                                    cGAACAATTTTTCGTATAACGTCTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGc                                                                                                                                                                                                         <  2:82509/251‑1 (MQ=255)
                                            ttttCGTATAACGTCTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTGGCCTGCGat                                                                                                                                                                                                  <  1:73180/250‑1 (MQ=255)
                                                  tataACGTCTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGt                                                                                                                                                                                           >  1:133124/1‑251 (MQ=255)
                                                        gTCTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGccaag                                                                                                                                                                                     <  1:87902/251‑4 (MQ=255)
                                                            tAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAgccg                                                                                                                                                                                  <  2:56769/250‑1 (MQ=255)
                                                                                      cATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAAGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCt                                                                                                                                                      >  2:238293/1‑251 (MQ=255)
                                                                                                     tatCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGGGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAACCGCCC‑GCCGGTTCATTCAGTCCAAAACCTTGCTCGTTCAGGCTg                                                                                                                                       >  2:10704/1‑251 (MQ=255)
                                                                                                         gCCCAATTCATAGGTCCTAACGGCGGCAGCGCCACGGGCTAAATCCCCCGTTCG‑AATGCCTTCTTCTAATGCGCGGATAATGGAGCGGTCAATGGCGCAATCCGCATCAACGGCATCCAGGCTGTCA‑CGAAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCACGCCGCCC‑GCCGGTTCATACAGTCCACAACCTTGCTCGTTCAGGCTGGCgg                                                                                                                                   >  2:323194/1‑251 (MQ=255)
                                                                                                              tgtcgtCGGTGCTAAGGGCGGCAGCGCCAGGGGTTAACTCCCCGTTGCA‑AATGCCTTCTTCTAAGTCGGGGGTAATGGCGCGTGCAATGGCGCAAGCCGCATACTCGGCATCCAGGTTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGTCTATTTTTTTGCCTGCGATGTCTGGTGCCAGGGCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTTGCTGAAgg                                                                                                                               <  2:313998/245‑1 (MQ=255)
                                                                                                                 cATCGGTACTAACGGCGGCTGCGCCACGGGCTAAATCCCGGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAAGGGCGCGTTCAAGGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCGGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAAc                                                                                                                           <  1:284382/251‑1 (MQ=255)
                                                                                                                 cATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAAc                                                                                                                           <  1:203707/251‑1 (MQ=255)
                                                                                                                    cGGTACTACCGGCGGAAGCGACAAGGGCGAACTCCCCGGTGCG‑AATGACTTCTTCTAATGCGCCGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCACTCAGGCTTAAA‑CGCAGCAGCAGTGCTAGCGCAAGGATTTGTGCAATCGGGTTGGCGAAGGTTTTGCCTGCGAAATATGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTACAAAACCTTTCTCGTTCCAGCTGGCGGGAggacacatc                                                                                                                        >  2:274925/1‑244 (MQ=255)
                                                                                                                    cGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGCGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCgg                                                                                                                                                                         >  1:60338/1‑202 (MQ=255)
                                                                                                                    cGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGCGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTGTGCCGGCGATATCGGGTGCCAAGCCGCCC‑GCCgg                                                                                                                                                                         <  2:60338/202‑1 (MQ=255)
                                                                                                                         cTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAACAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGGTTTTTCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTACAAAACCTTGCTCGTTCAGGCTGGCGGCAGGCAACATCCCCAt                                                                                                                   >  1:299966/1‑251 (MQ=255)
                                                                                                                           aaCGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGGTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGAAGTGCTAGCGAAAGGATTTGTGCAATCGGGGTGGCGATGGTTTTGCCTGCGATATCTGGGGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGGTCAGGCTGGCGGAAGGCAACATCCCCATCg                                                                                                                 >  1:43220/1‑251 (MQ=255)
                                                                                                                                cggcAGCGCCACGGGTTAAACCCCCGGTGCG‑AATGCTGTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCTAAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGAAGCAGTGCTAGCGAAAGGTTTTGTGCAATCGGGTTGGCGATGTTTTTGCCGGCGATATCGGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCa                                                                                                            <  1:275994/251‑1 (MQ=255)
                                                                                                                                     gcgccaggggtaaaaTCCCCGG‑GCGGAATGCCTTCTTCAAATGCGCGGTTAATGGCGCGTTCAATGGCGCAACCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGTTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGAt                                                                                                       <  2:103866/239‑1 (MQ=255)
                                                                                                                                       cccACGCGCTAAATCCCCGGTGCG‑AATGCTTTCTTCTAATGCGCGGTTTATGCCGCGTTCAGTGGTGCAAGCCGCTTTATCAGCAGCCAGGCTGTAA‑CGCAGCAGCAGGTCTAGCTAAAGTATTTTGGCAATCGGGTTGGCTATTTTTTTGCCTGCGATATCTGTTCCCAAGCCGCCC‑GCCGGTTCATACAGTCCAACACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCa                                                                                                     <  1:175167/250‑1 (MQ=255)
                                                                                                                                            gggCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTgcg                                                                                                >  1:61545/1‑251 (MQ=255)
                                                                                                                                                 aactCCCCGGGGCG‑AGTTCCGTCTTCTAATGCGCGGTTAAGGGCGCGTTCAGTGGCGCAAGCGGCATCATCGTCATCCAGGCTGTAA‑CGCAGCTGCAGTTCTAGCGAAGGGATGTGTGCAATCGGGTTGCTGTGGTTTTTGCCGGCGATAGCTGGTGCCAAGCCGCCC‑GCCGGTTCAGACATGCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTc                                                                                           <  1:187343/248‑1 (MQ=255)
                                                                                                                                                 aaacccccGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGACTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTc                                                                                           <  1:242340/247‑1 (MQ=255)
                                                                                                                                                 aaacccccGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAAGGGGGGGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAACCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTc                                                                                           <  1:75431/247‑1 (MQ=255)
                                                                                                                                                       ccGGTGCG‑ATTGCCTTCTTCTAATCCGCGGTTAATGGCGCGTTCAAGGGCGCAGGCCGCACCATCGGGAACCAGGCT‑TAATCGCAGCAGCAGTGCTACCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCGGCGATATCGGGTTCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGACAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTcaga                                                                                     <  1:312695/251‑1 (MQ=255)
                                                                                                                                                          gTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTcagacag                                                                                  <  1:16685/251‑1 (MQ=255)
                                                                                                                                                                 aTGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGaa                                                                                                                                 >  1:270177/1‑198 (MQ=255)
                                                                                                                                                                   gCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGc                                                                          >  1:172428/1‑251 (MQ=255)
                                                                                                                                                                     cttcttCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGcc                                                                         >  1:111661/1‑250 (MQ=255)
                                                                                                                                                                               tGCGAGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGAAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGCCGCCAAACAGGTTgg                                                              >  1:117124/1‑251 (MQ=255)
                                                                                                                                                                                        aagggCGCGTGCAAGGGCACAAGCCGTATAACCGCCAGCCAGGCTGAAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCGGGGGCCAAGCCCCCC‑GCCGGTTCATACAGTCCAAAACCGTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCAcagc                                                      <  1:315378/247‑1 (MQ=255)
                                                                                                                                                                                             cgcgTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACg                                                >  1:110578/1‑251 (MQ=255)
                                                                                                                                                                                                  acaaggccgcACGCCGCATCATCGGCATCCAGGTTGTAG‑CGCGGCAGCAGTGCTAGCGAAAGGTTTTTTGCAATCGGGTGGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCaaa                                           <  2:126484/244‑1 (MQ=255)
                                                                                                                                                                                                            aaGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATgg                                  >  1:186457/1‑250 (MQ=255)
                                                                                                                                                                                                            aaGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCAAGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCACAACCTTGCTCGTTCAGGCTGACGGAAGGCAACATCCCCATCGCGCAAGTGATCATTGCGCCCTCGTCAGCCAGAATGTCGCCAAACAGGTTGGAGCACACCATAACGTCAAACTGTGATggc                                 >  1:186747/1‑250 (MQ=255)
                                                                                                                                                                                                                 gCATCATCGGCATCCAGGCTGTAA‑AGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGCTATCTGGTGCCAAGCCGCAC‑GCCGGTTAATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCCGTGATCATTGCGACCTCGTCAGACAGAATGTCGCCAAACAGGATGGCGCACAGCAGAACGTCAGACTGAGATGGATCttt                            >  2:274216/1‑251 (MQ=255)
                                                                                                                                                                                                                            cacccaggcgttaccGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAAGCGGGTTGGCGATGTTTTTGTCTGCGATATCGGGTCCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGc                  <  2:259188/237‑1 (MQ=255)
                                                                                                                                                                                                                                   cTGTAA‑CGCAGCAGCAGTGCTAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGAATGGTGGc          >  1:325047/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cagcagcagTGCAGGCGAAAGGAGTTGTTCAATCGGGTGGGCGATGTTTTTGCTGGTGATATCGGGTGCC‑AGCCCCCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGAGCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGa  <  1:323194/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
GGGTTTCGTTTTCGGCTTCGTACACAACGTGAGCGTCGAACAATTTTTCGTATAACGTCTTAGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCG‑AATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAA‑CGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCGAGCCGCCC‑GCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGA  >  NC_000913/80798‑81282

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