Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,152 A→C L142R (CTC→CGC)  thiF ← adenylyltransferase, modifies ThiS C‑terminus

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,1520AC76.9% 87.2 / 29.6 52L142R (CTC→CGC) thiFadenylyltransferase, modifies ThiS C‑terminus
Reads supporting (aligned to +/- strand):  ref base A (6/6);  new base C (18/22);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

TCCCGCAGGTGTCTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCT  >  NC_000913/4192926‑4193381
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      
tcCCGCAGGTTTCTCTATCCCGCTTACTAACTTAATGGCTACCAGTGCCTGCAAAGTGCCCAACACCCCGCCCCCCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAgg                                                                                                                                                                                                                            <  2:88509/239‑1 (MQ=255)
                  aCCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCCGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCTCCCAATCCGACCGCGCTGGCGGGGATAATCGGCGGGTGGCGTGCCCCGCCGCCGGCACTCATCTCCTGGc                                                                                                                                                                                                          >  2:668809/1‑239 (MQ=255)
                    cgttatataaCTTATTGGCTCCCGGTGCCGGCAAGGTGCCCATACCCCCAACCCCCGGGCCAACCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGa                                                                                                                                                                                                        <  1:1169016/233‑1 (MQ=255)
                     gtttAAAAACTAAATGGTTCCCAGTGCCGGAAAAGTCCCCATACCCCCGCCCCCCGGCCCACCCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAg                                                                                                                                                                                                       <  2:261939/237‑1 (MQ=255)
                     gCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATACCCCCGACCACCGGGCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAg                                                                                                                                                                                                       <  1:808703/239‑1 (MQ=255)
                      cTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATACCCTGGCgcg                                                                                                                                                                                                       >  1:761500/1‑236 (MQ=255)
                        taataaCTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGGGTTGCGTGCCACGCAGGCGGCCTTAATATCCTGGGGAGTc                                                                                                                                                                                                     >  2:175655/1‑238 (MQ=255)
                         aataaCTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGGGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCg                                                                                                                                                                                                    >  1:737173/1‑238 (MQ=255)
                               ttAATGGCTCCCAGGCCCGACAAGGGCCCCATACCCCCGACCCCCGGCCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATAtt                                                                                                                                                                                             <  2:370056/239‑1 (MQ=255)
                                 aatgttctCCAGTGCCGGCAAGGTGCCCATCCCCCCGCCCACGGGGCCAACCCCCCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGt                                                                                                                                                                                           <  2:896977/232‑1 (MQ=255)
                                 aaTGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTTCGTGCCCAGCAGGCGGCCTTAATCTCCTGGCGAGTCGccctgtttt                                                                                                                                                                                           >  2:69185/1‑232 (MQ=255)
                                 aaTGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCat                                                                                                                                                                                                >  1:155779/1‑234 (MQ=255)
                                              cccTGCAACGCGCCCTGACCCCCGACCACCGGCCCACCCGCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTc                                                                                                                                                                                 <  2:761500/235‑1 (MQ=255)
                                                          cccATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCCGGCGGCATTAATCTCCTGGCGCGTCGCCCTATTTTCGGTACAGTCGAGCACCACATCgg                                                                                                                                                                   >  2:450610/1‑238 (MQ=255)
                                                                     gACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCCTTAATCTCCTGGCGAGGCGCCATATTGTCGGGACAGTCGAGGACCACATCGGGCCGGGCaaac                                                                                                                                                       >  2:266458/1‑236 (MQ=255)
                                                                                      cgcccgccGGGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACCCCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCa                                                                                                                                     <  1:450610/239‑1 (MQ=255)
                                                                                        cccgccGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGGGTCGCCCTATTGTCGGTAACGTTCAGGACCACAACGGCCCCGTCCACCCCCACTTTTTAACCCTCAc                                                                                                                                    >  2:769384/1‑238 (MQ=255)
                                                                                                gcggccagtGCGTTCGGGCTCCCGGTTAGCTGGTCCAAGGCCGGGGTTACCCCCCTGGCCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAAc                                                                                                                            <  1:130316/230‑1 (MQ=255)
                                                                                                             tCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGc                                                                                                              <  2:155779/239‑1 (MQ=255)
                                                                                                                    ccTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGCCGCCATATTGTCGGTACAGTCGAGCACCACCTCGGCCCGGGCAACCGCATCTTTTAACCCCTCACCCGTTAAACCttgttggtat                                                                                                                >  2:732494/1‑226 (MQ=255)
                                                                                                                            tCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTc                                                                                               >  1:119773/1‑239 (MQ=255)
                                                                                                                                  cacaGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCCCATCGGCCCGTGCAACCGCATCCTTTAACGCCTCCACCGTTAACCCTTGTTGTAATGCTTTCAGTTTAATTTCGGGAtt                                                                                         >  2:588692/1‑239 (MQ=255)
                                                                                                                                   acaGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAg                                                                                                                                                                                   >  2:533125/1‑149 (MQ=255)
                                                                                                                                   acaGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAg                                                                                                                                                                                   <  1:533125/149‑1 (MQ=255)
                                                                                                                                          gCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGCTAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCCTCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCa                                                                                       >  1:197198/1‑233 (MQ=255)
                                                                                                                                           cGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGGGTTGAGTGCCACGCAGGCGGCATTAATTTCCTGGCGAGTCGCCATATTGTCGGGACAGTCGAGCACCACATCGGCCCCGGCAACCGCATCTTTTTAAGCCTCCCCCGTTAACCGTTGGTGGAATGCTTTCAGGTTGATGTTGGGGTTCacac                                                                                    >  2:182597/1‑232 (MQ=255)
                                                                                                                                                            ccAGGGCTGCGGCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTg                                                                <  1:486185/238‑1 (MQ=255)
                                                                                                                                                              aGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAgg                                                                                                                                                                                                                            <  2:924205/81‑1 (MQ=255)
                                                                                                                                                              aGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAgg                                                                                                                                                                                                                            >  1:924205/1‑81 (MQ=255)
                                                                                                                                                                     cgTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCg                                                                                                                                                                                >  1:213498/1‑118 (MQ=255)
                                                                                                                                                                     cgTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCg                                                                                                                                                                                <  2:213498/118‑1 (MQ=255)
                                                                                                                                                                             ccATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGACTTCGGGcga                                              <  1:175655/239‑1 (MQ=255)
                                                                                                                                                                             ccAGCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcga                                              <  1:909816/239‑1 (MQ=255)
                                                                                                                                                                                tCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatcg                                           >  1:210297/1‑239 (MQ=255)
                                                                                                                                                                                  aaCTGGCCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatcgat                                         <  1:588692/239‑1 (MQ=255)
                                                                                                                                                                                       aCCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGGGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTTACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCTGGCGCTCGCTATCt                                     >  2:176624/1‑238 (MQ=255)
                                                                                                                                                                                          gCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAg                                 <  2:119773/239‑1 (MQ=255)
                                                                                                                                                                                          gCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGc                                                                                                                                                                 >  1:719938/1‑111 (MQ=255)
                                                                                                                                                                                          gCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGc                                                                                                                                                                 <  2:719938/111‑1 (MQ=255)
                                                                                                                                                                                            cAAAGCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCtt                                    <  1:176624/234‑1 (MQ=255)
                                                                                                                                                                                                tgcgcgcgtgcTGGCGGTGATGAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCTAGTCGCCATATTGTCGGTACCGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCGCACCCGGTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTTGCTGGCCTGCGATTTGGGGCGATCGATATCTTCAGTGGTaa                           <  2:144258/230‑1 (MQ=255)
                                                                                                                                                                                                    aCCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACAt                                                                                                                                                                      >  1:704719/1‑97 (MQ=255)
                                                                                                                                                                                                    aCCGCGCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACAt                                                                                                                                                                      <  2:704719/97‑1 (MQ=255)
                                                                                                                                                                                                      cgcgCTGGCGGTGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAgg                      <  1:81024/238‑1 (MQ=255)
                                                                                                                                                                                                                    tAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGa       <  2:210297/239‑1 (MQ=255)
                                                                                                                                                                                                                       gcggcgTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTa                                                                                                                  >  2:273151/1‑129 (MQ=255)
                                                                                                                                                                                                                       gcggcgTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTa                                                                                                                  <  1:273151/129‑1 (MQ=255)
                                                                                                                                                                                                                       gcggcgTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACAt                                                                                                                                                                      <  1:770256/78‑1 (MQ=255)
                                                                                                                                                                                                                       gcggcgTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACAt                                                                                                                                                                      >  2:770256/1‑78 (MQ=255)
                                                                                                                                                                                                                          gcttgtTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCt  <  1:69185/235‑1 (MQ=255)
                                                                                                                                                                                                                           cgTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcga                                              >  2:635280/1‑193 (MQ=255)
                                                                                                                                                                                                                           cgTGTTGCGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcga                                              <  1:635280/193‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      
TCCCGCAGGTGTCTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCT  >  NC_000913/4192926‑4193381

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