Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,196,220 C→T H366H (CAC→CAT)  icd → isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,2200CT95.8% 132.6 / NA 48H366H (CAC→CAT) icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (13/33);  minor base A (2/0);  total (15/33)
Fisher's exact test for biased strand distribution p-value = 9.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCA  >  NC_000913/1195922‑1196448
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
aagaCCCGAACACTGGGAAAGAGTTCTGAAGTGAAGAGGTGATTGCTGATGCATTCCTGGTACAGGTCCTGCTGCGTCCGGCTGAAGAGGAGGTTATCGCCTGGATGAACCTGAACGGTGACTACATTTCTGACGCCCGGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCatat                                                                                                                                                                                                                                    <  1:17208/298‑1 (MQ=255)
    cccGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCCCCCACGGTACTCCGCCGAAATATGCCGGTCAGGACAACGTCAATCCTGGCCCTATTATTCTCTCCGTTTGAG‑TGGTGCTGCGCCAACTGGGt                                                                                                                                                                                                                                >  1:762962/1‑301 (MQ=255)
     ccGAACACTAGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGtt                                                                                                                                                                                                                               <  2:34774/301‑1 (MQ=255)
         acacTTGCAAAGAGTGCGTCATTAAAGACGTGATTGCTGCTGCATTCCTGCAACAGATCCTGCGGCGTCCGGCTGAATATGAGGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTTTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGAc                                                                                                                                                                                                                           <  2:808308/301‑1 (MQ=255)
                    gagaTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGCAGCCACCCACGGTACTGCGCCGAAATATGCCGGTTAGGACAAAGTAAATCCTGGCCCTATTCTTCTCTCCG‑CTGAGATGATTCTGCGCCATATGGGTTGGACTGAAgggg                                                                                                                                                                                                                   >  2:837436/1‑298 (MQ=255)
                                     cctgagtgcgcaTGCATTCCTGCAACTGATCCTGCTGCGTCCGGCCTAATATCATGGTATGTCCTGTATGAACCTGAACGGTGACTACATCTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGt                                                                                                                                                                                               <  2:84058/290‑1 (MQ=255)
                                     cGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGt                                                                                                                                                                                               <  2:778840/301‑1 (MQ=255)
                                     cGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGACGTTATCGCCTGTATTAACCTGAACGGTGACTACATGTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCTCTGGTGCAAACATCGGTGACGAATTCGCCCTGTTTGATGCCACCCTCGGCACTGCGCCGAAATATCCCGGTCAGGTCAAAGTAAATCCTGGCTCTTTTGTTCTCTCCG‑CTGTGAGTTTGCTTCCCCAAATGGGTCGGACTGAAGCGGCTGAccgtacttggaaaggg                                                                                                                                                                                               >  1:456994/1‑286 (MQ=255)
                                           tGCTGATTCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCGGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGaa                                                                                                                                                                                         <  1:137669/301‑1 (MQ=255)
                                           tGCTGATGCATTCCTGCAACAGCTCCTGCTGCGTCCGGCTGAATATTATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGaa                                                                                                                                                                                         <  2:457007/301‑1 (MQ=255)
                                            gCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGCCCAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCCCCATATGGGGTGGACTTAAGCGGCTGCCCTGATTGGTAAAGGTATgggaa                                                                                                                                                                                        >  1:349742/1‑298 (MQ=255)
                                            gCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCTCCATATGGGGTGGACTGAAGCGGCTGACCTGATTGGTAAAGGTATGGaaa                                                                                                                                                                                        >  1:305948/1‑300 (MQ=255)
                                                aTGCATTCCTGCAACAGGTCCTGGTTCGTCCGGCTGAATATGATGTGATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcgc                                                                                                                                                                                    <  2:433962/301‑1 (MQ=255)
                                                aTGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGAATGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                                                                                                                                                                     >  2:874912/1‑300 (MQ=255)
                                                        cTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaatg                                                                                                                                                                            >  2:396414/1‑299 (MQ=255)
                                                               tgagccgggggtgtcacgCTGAATATGAGGGTATGGCCTGTATGAACCTGCACGGTGAATACATTTCTGAAAGCC‑GGCACGGGCAGGTTGGGGGGAACGGTTTCGCCCCTGGTGCAAAAATCGGTGACGAATGCGCCCTGATTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAGATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAGAGGTATGGAAGGCGCAATCAATGCCAAGAc                                                                                                                                                                     <  1:270277/284‑1 (MQ=25)
                                                                 aTCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTTCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCg                                                                                                                                                                   >  2:687194/1‑301 (MQ=255)
                                                                      gctgcGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAc                                                                                                                                                               <  2:57109/300‑1 (MQ=255)
                                                                       ttgcGTCCGGCTGAATATGATGTTTACGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCG‑GGGCAGGTTGGCGGTATCGGTATCGCCCCTGGTTCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAActt                                                                                                                                                             <  1:687168/300‑3 (MQ=38)
                                                                                          atgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGGCTgat                                                                                                                                          >  2:103098/1‑301 (MQ=255)
                                                                                           tgtgctCGCCAGTATGACGCTGAACGGTTGTCACATGTCGGACGACCTGGCA‑GCGCAGGTTGGCGGGATCGGTATCGCCCCTGGTGCAAAAATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatg                                                                                                                                         <  1:291961/296‑1 (MQ=25)
                                                                                           tgTTTTAGCCCGTATGAACCTGAACGGTGACTACATTGCTTACGCCCGGGCA‑GCGCAGGGTGGCGGGATCGGGAGCGCCCCTGGTTCAAACATCGGTGACGAAGTCGCCCTGTTTGAAGACACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTAGTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatg                                                                                                                                         <  1:730673/301‑1 (MQ=25)
                                                                                            gTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGCTGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAACGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACGTATGACGTcagacttttgatgg                                                                                                                                        >  1:845976/1‑288 (MQ=38)
                                                                                             ttATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatgg                                                                                                                                        <  2:297141/300‑1 (MQ=255)
                                                                                             ttATCGCCTGTATGAACCTGAACGGTGACTACATTTCGGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatgg                                                                                                                                        <  2:196608/300‑1 (MQ=255)
                                                                                                  gCCTGTATGAACCTGAACGGTAACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGCGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGCTGTTGCTGCGCCATATGGGTTGGTCTGCAGCGGCTGACCTGCTTGTTAAAGGTAGGGAACGCGCAATCAATGACACGCCCGTAACTTATGACTTCGCACGCCTGAGGGAAGgcg                                                                                                                                  >  2:99320/1‑301 (MQ=255)
                                                                                                      gTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTaa                                                                                                                              >  2:361161/1‑301 (MQ=255)
                                                                                                        tcaacactaaccggtggtcacgtttCGGACGCCCTGGCA‑GCGGAGGTTGGCGGTACTGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCGCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCtgagc                                                                                                                            <  1:78165/279‑5 (MQ=21)
                                                                                                        aTGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTaagc                                                                                                                            <  2:765321/301‑3 (MQ=37)
                                                                                                        aTGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTaagc                                                                                                                            <  2:445294/301‑3 (MQ=37)
                                                                                                             ccTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCATGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGATATGGAAGGCGCAATCAATGCCAAGAACGTAACTTATGAATTCGCAAGGCTGCTGGAAGGCGCTAAGccgcta                                                                                                                       >  1:743553/1‑300 (MQ=37)
                                                                                                                         ttccATGTCTGACGCGCTGGCA‑GAGCAGGTTGGCGGTATCGGGATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCTAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCagag                                                                                                           <  1:655987/298‑1 (MQ=21)
                                                                                                                                    aCGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGaag                                                                                                 <  2:477317/300‑3 (MQ=25)
                                                                                                                                                           cggtatcggtatCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtc                                                                           <  2:762990/300‑2 (MQ=21)
                                                                                                                                                                   gtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTTTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtctctccatg                                                                   <  2:826949/300‑10 (MQ=18)
                                                                                                                                                                    tatcgCCCCTGGTGCAAACATCGGTGACGAATGCGACCTGTTTGAAGCCACCCACGGTACTGCGCCGAAAGATGCCGGGCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGAGGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtctctccatgtg                                                                 <  1:396403/301‑12 (MQ=17)
                                                                                                                                                                          cccTGGTGCAAACAGCCGGGACGAATCAGCCCGTGTTGAAGCCACCCCCGTTACTGCGCCGAAATATGGCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCGGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtctctccatgtgttaaat                                                           <  1:103093/301‑18 (MQ=14)
                                                                                                                                                                          cccTGGTGAAAACATAGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAATTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtctctccatgtgttaaat                                                           <  2:806189/301‑18 (MQ=17)
                                                                                                                                                                          cccTGGGGCAAAAATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtctctccatgtgttaaat                                                           <  2:846008/301‑18 (MQ=17)
                                                                                                                                                                                                               aaaGCCACCCACGGTAATGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCGCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttattt                      <  1:361151/300‑5 (MQ=11)
                                                                                                                                                                                                                          cGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatatt           >  2:265037/1‑286 (MQ=14)
                                                                                                                                                                                                                            ttACTTCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattat         <  2:197130/300‑18 (MQ=14)
                                                                                                                                                                                                                             tACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTTTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattat         <  2:349752/300‑18 (MQ=14)
                                                                                                                                                                                                                              tcTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGCCTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattatta       <  1:265027/300‑20 (MQ=14)
                                                                                                                                                                                                                                tGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattattaat     <  2:305958/301‑22 (MQ=14)
                                                                                                                                                                                                                                 gcgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtc                                                                           >  2:630980/1‑229 (MQ=14)
                                                                                                                                                                                                                                 gcgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtc                                                                           <  1:630960/230‑2 (MQ=14)
                                                                                                                                                                                                                                 gcgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattattaat     <  2:467450/300‑22 (MQ=14)
                                                                                                                                                                                                                                    ccGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattattaatagc  <  2:901291/300‑25 (MQ=17)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
AAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCG‑CTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCA  >  NC_000913/1195922‑1196448

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