Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,041,467 A→G K31K (AAA→AAG zinT → metal‑binding protein ZinT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,041,4670AG97.5% 136.9 / NA 40K31K (AAA→AAGzinTmetal‑binding protein ZinT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (18/21);  minor base T (0/1);  total (18/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 = 9.99e-01

TAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGC  >  NC_000913/2041238‑2041713
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
acgggactgggtggagaaacactgtcgaagagtgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACaa                                                                                                                                                                                                                                   <  2:94025‑M2/110‑1 (MQ=255)
     actggggggagaaacactgtcgaagattgaatggcttgcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAgc                                                                                                                                                                                                                              <  2:451533‑M2/115‑1 (MQ=255)
                            agattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCCTGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGCACAAAAAGATGCTAATGGGGTTTTTGatgatg                                                                                                                                                                                                       >  2:468547‑M2/114‑251 (MQ=255)
                                         ggcgctaggcgtaattgtcttcgcggtgggantggtcgcgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGGCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTGTGATGATGCCAATGTACaaaa                                                                                                                                                                                          <  2:230436‑M2/151‑1 (MQ=255)
                                          gcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAAc                                                                                                                                                                                         <  1:154176‑M2/152‑1 (MQ=255)
                                            gctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCg                                                                                                                                                                                       <  1:132376‑M2/154‑1 (MQ=255)
                                                     aattgtcttcgcggtggtactggtcacgtttggcaaatatctcttccccgcaaaacccgtagtggcgccagttattcaggacgcatcaCGATTCGTCTTTACAATCTGGCTGTTGCTTTAGTTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGTCAATCCCTTAACAGAGGTCGATCAAAAAGCTTCTAATGGTGTTTGTGATGATGCCAATGTACAAAACCGAACGCtct                                                                                                                                                                               <  2:287305‑M2/162‑2 (MQ=255)
                                                                  gtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGGGACTGGGAtggg                                                                                                                                                                 >  2:254348‑M2/76‑250 (MQ=255)
                                                                     gtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCCCAGTGACTGGGATGGAGTc                                                                                                                                                              >  2:47653‑M2/73‑251 (MQ=255)
                                                                                    ggaaaatatctcgtcccggcaaaacccgtagttgcgacagttattgtggacgcatctCGATTCGTCTTTACAAGCTGGCTGTTGCTTTAGGTTTCTTTATTGTTAGCGCGCCTGCCTTTTCGCATGGTCATCACTCACGCGGCAAGTCCTTAACAGAGGTCGAACAAAATGCTGCTAAGGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCATTGTCTGGGATGGAGTCTGGCAATCCGTTTAt                                                                                                                                               <  2:476419‑M2/194‑1 (MQ=255)
                                                                                       aaataacgctgtccggcaaaacccgtagttgtgccagttattcaggacggatcaCGATTCGTCTGTACAAACTGGCTGTTGCTTTAGGTGTCTTTAGTGGTAGCGCTCCTGCCGTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCt                                                                                                                                            <  2:106212‑M2/197‑1 (MQ=255)
                                                                                        aatatctcttcccggcaaaacccgtcgttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCtt                                                                                                                                           >  1:348930‑M2/54‑251 (MQ=255)
                                                                                                      gcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGc                                                                                                                             >  1:336269‑M2/40‑251 (MQ=255)
                                                                                                             ccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTg                                                                                                                      <  1:193173‑M2/219‑1 (MQ=255)
                                                                                                             ccgtagttgcgccagttagtcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACGGCAAAGTGGCAAACTTg                                                                                                                      <  1:351442‑M2/219‑1 (MQ=255)
                                                                                                                          agttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTc                                                                                                          >  2:109984‑M2/20‑250 (MQ=255)
                                                                                                                             tattcaggacgactcaAGTTTCGTCTTTACAAACTGGCTTTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGACGCCCATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCaga                                                                                                       <  2:263723‑M2/234‑1 (MQ=255)
                                                                                                                                 caggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTTCTCATGGTGTTTTTGATGATGCCACTGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAgaagaa                                                                                                   >  2:185575‑M2/13‑250 (MQ=255)
                                                                                                                                  aggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAgaagaa                                                                                                   >  2:162973‑M2/12‑249 (MQ=255)
                                                                                                                                      cgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGa                                                                                             <  1:109984‑M2/244‑1 (MQ=255)
                                                                                                                                           caCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAg                                                                                        <  1:8159‑M2/249‑1 (MQ=255)
                                                                                                                                                  cGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAtaaaac                                                                                 <  2:393771/251‑1 (MQ=255)
                                                                                                                                                        tACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAtaaaactaaaa                                                                            <  1:360660/250‑1 (MQ=255)
                                                                                                                                                              cTGGCTGTTGCTTTATGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGACAGCGGATGCAGATAAAACTAAAACATTTGc                                                                     >  2:327775/1‑251 (MQ=255)
                                                                                                                                                               tGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTg                                                                                                                                                            >  1:400157/1‑163 (MQ=255)
                                                                                                                                                               tGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCACTGACTGGGATGGAGTCTg                                                                                                                                                            <  2:400157/163‑1 (MQ=255)
                                                                                                                                                                 gCTGTTGATTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGGCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCTGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGa                                                                  <  1:211058/251‑1 (MQ=255)
                                                                                                                                                                    gTTGCTTTAGGTGTCTTTATTGTTAGCGATCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCCGCGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGGGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAAt                                                               >  2:296639/1‑251 (MQ=255)
                                                                                                                                                                      tGCTTTAGGGGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAACCATTTGCTGAAAtt                                                              <  2:139581/250‑1 (MQ=255)
                                                                                                                                                                        cTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTAATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAATGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATTCAGATAAAACTAAAACATTTGCTGAAATTaa                                                            >  1:210283/1‑250 (MQ=255)
                                                                                                                                                                             ggTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGattt                                                      >  1:65648/1‑250 (MQ=255)
                                                                                                                                                                                    ttATTGGTAGCTCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATc                                                 <  2:279494/249‑1 (MQ=255)
                                                                                                                                                                                      aTTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATcaca                                              >  1:435429/1‑250 (MQ=255)
                                                                                                                                                                                         gTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATACAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAgg                                          <  1:264398/251‑1 (MQ=255)
                                                                                                                                                                                                            cGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCg                         >  2:283243/1‑249 (MQ=255)
                                                                                                                                                                                                             gCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGAc                                                                                                                    >  2:3536/1‑157 (MQ=255)
                                                                                                                                                                                                             gCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGAc                                                                                                                    <  1:3536/157‑1 (MQ=255)
                                                                                                                                                                                                               aTGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg                    <  1:355278/251‑1 (MQ=255)
                                                                                                                                                                                                                         ctcacACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTg          >  1:165771/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gCAAGCCCTTAACAGCGGTCGAACAAAAACCTGCTAATGGTGTTTTTGATGATGCAAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAACACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAAAAGATATCGAGATTATTGGCATTGAGGACGGc  >  2:273675/1‑251 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
TAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGACGGC  >  NC_000913/2041238‑2041713

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