Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 555,049 C→T P147S (CCG→TCG)  cysS → cysteine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913555,0490CT96.4% 176.7 / ‑3.4 55P147S (CCG→TCG) cysScysteine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (35/18);  total (37/18)
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.35e-01

CATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCCCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGG  >  NC_000913/554808‑555284
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           
cATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATc                                                                                                                                                                                                                                    <  1:202136/251‑1 (MQ=255)
      cGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACtt                                                                                                                                                                                                                              <  2:209489/251‑1 (MQ=255)
       gATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACACCGGCGACGTGATGTTCGACGTCTCGACCGATCCAACtt                                                                                                                                                                                                                              >  1:184173/1‑250 (MQ=255)
             gacgacAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGt                                                                                                                                                                                                                       >  2:208718/1‑251 (MQ=255)
                                gcgcCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGAtctc                                                                                                                                                                                                    <  1:179478/251‑1 (MQ=255)
                                   ccAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGAc                                                                                                                                                                                                 >  1:200821/1‑251 (MQ=255)
                                     aaTGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGAcc                                                                                                                                                                                                >  2:279383/1‑250 (MQ=255)
                                      aTGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCac                                                                                                                                                                                              >  2:124222/1‑250 (MQ=255)
                                        gaaaaaggcgaaatagttGTGGAGATTGTGGATCGCATGATCGACGAAATGCGGAAAGATTTTGATGCTTTGAACATGCTGCGCACGGATATGGAGCCGCGCGCGACGGACGATATCGCAGAAATTAGTGAAATCACTGAACAACTGATCGCCAAAGGTCACGCATAGGTGGCGGATAACGGCGACGTGATGTTCGACGTCTCAGCCGGTCCAACTTATGGCGTGCTGTCGCGTCTGGAGCGCCACCAGCt                                                                                                                                                                                            <  2:177093/235‑1 (MQ=255)
                                           aaTGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTgc                                                                                                                                                                                          >  1:228568/1‑250 (MQ=255)
                                                cGAAAGCTTTGTGGAGATGGTGACTCGCATGATCGCCGAAAAGAACAAAGGTTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTACTGTCGCGTCAGGATCTCGACCAGCTgcagng                                                                                                                                                                                      >  2:260444/1‑247 (MQ=255)
                                                 gAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGAGTGCTGTCGCATCAGGATCTCGACCAGCTGcaggcagg                                                                                                                                                                                   >  1:67/1‑251 (MQ=255)
                                                                 tggtggATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCGGGCGGGCGCGTTggcggg                                                                                                                                                                   >  1:235707/1‑246 (MQ=255)
                                                                  ggtggATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTg                                                                                                                                                                   <  1:45614/250‑1 (MQ=255)
                                                                                  gCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCa                                                                                                                                                   >  1:92980/1‑250 (MQ=255)
                                                                                     gAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACCAACGCAAccc                                                                                                                                               >  2:107251/1‑251 (MQ=255)
                                                                                       aaTGCACCAAGATTTTGATGCTTTGAACATTCTGCGCACGGATATGGAGCCGCGCGCGACGCACCATATCGCAGACATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGGCGACGACAAACGCAAccaa                                                                                                                                              >  2:42621/1‑248 (MQ=255)
                                                                                           cacaAAGATGTTGATGCTTTGAACATTCTTCGCCCGGATATGGAGCCGCGCGCGACGCAACATAGCGCAGAAATTATTGAACGCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTTGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGCCCAGCGTCAGGCAGGCGCGCTCGTTGACGTGGTCGACGACAAACGCAACCCAATGGa                                                                                                                                         <  2:160001/251‑1 (MQ=255)
                                                                                                  aTTTTGATGCTTTGAACATTCTGCGCCCGGATATGGACCAGGGCGCGACGCACCATATAGCAGAAATTTTTAAACTAACTGAACAACTGATCGCCAAAGGCCACGTTTATGTGGCGGACAAAGGCGACGCGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGTTGTCGCGTCAGGTTTTCGACCTGCTGAAGGCCGGCGCGGGGGTTACCGGGGGCGACGACAAACGCACCCCCACGGCCTTCGtt                                                                                                                                  >  2:152371/1‑251 (MQ=255)
                                                                                                       gATGCTTTGAACATTCTGCGCCCGGATGTGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGGCGACGACAAACGCCACCCAATGGACTTCGGTCtgtg                                                                                                                             >  2:44763/1‑251 (MQ=255)
                                                                                                        aTGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCCCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGCCAAACGCAACCCAATGGACTTCGTTCTGTgg                                                                                                                            >  1:168820/1‑251 (MQ=255)
                                                                                                             ttGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAg                                                                                                                         <  2:244269/249‑1 (MQ=255)
                                                                                                             ttGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGGGGTCGACGACAAACGCAACCCAATGGACCTCGTTCAGTGGa                                                                                                                           >  2:160341/1‑247 (MQ=255)
                                                                                                              tGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCGGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGAt                                                                                                                       <  1:107548/250‑1 (MQ=255)
                                                                                                              tGAACATTCTGAGCCCGGATATGGAGCACCGCGCGAAGCACCATATCGCAGAAATAATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGACGACGTGAAGTTCCACGTGTCGACCGATCTAACTCATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCATGCAGGCGGGGGCGGTGACGGGGGATACGGCAAACGCAACCCAGAGGACTCCCTTCTGTCGatgatg                                                                                                                      >  2:95401/1‑251 (MQ=255)
                                                                                                                     tCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTATGTGGAAGATGTCGaca                                                                                                                >  2:135171/1‑248 (MQ=255)
                                                                                                                     tCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGGTGACGGGGTCGACGACAAACGCAACCCCATGGACTTCGTTCTGGGGAAGATGTCGaa                                                                                                                 >  2:42370/1‑249 (MQ=255)
                                                                                                                              ggATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGc                                                                                                                                                                                                                     <  1:57470/140‑1 (MQ=255)
                                                                                                                              ggATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGc                                                                                                                                                                                                                     >  2:57470/1‑140 (MQ=255)
                                                                                                                              ggATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGCCGTGGTCGCCGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAAc                                                                                                      >  2:271402/1‑251 (MQ=255)
                                                                                                                                          gcgcgcgACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGCTCGCCAAAGGTCACGCTTATGTGGAGGACAACGGCGACGTGATGTTCGACGTCTAGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGt                                                                                          >  1:41004/1‑251 (MQ=255)
                                                                                                                                              gcgACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGCTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGCAAGAGGGCGAACCGAGCTGGCCGTCTcc                                                                                      >  2:212433/1‑251 (MQ=255)
                                                                                                                                                  cGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCCCGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCCATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTgg                                                                                  >  1:36371/1‑251 (MQ=255)
                                                                                                                                                     aCCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGgc                                                                               >  1:97057/1‑251 (MQ=255)
                                                                                                                                                        atatCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGcgcg                                                                            <  2:200821/251‑1 (MQ=255)
                                                                                                                                                         tatCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCgg                                                                           >  1:61398/1‑251 (MQ=255)
                                                                                                                                                         tatCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCgg                                                                           >  1:171347/1‑251 (MQ=255)
                                                                                                                                                          atCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCggg                                                                          <  1:44763/251‑1 (MQ=255)
                                                                                                                                                              cAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGgtcgt                                                                      <  1:134110/251‑1 (MQ=255)
                                                                                                                                                                     tattGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGCAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGG‑CGCGGGTCGTCctggctg                                                              >  1:69864/1‑251 (MQ=255)
                                                                                                                                                                               aCTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTAGTCCTGGCTGGCACATTg                                                      >  2:169310/1‑250 (MQ=255)
                                                                                                                                                                                  gAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGaa                                                    <  1:40526/249‑1 (MQ=255)
                                                                                                                                                                                   aacaacTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCCGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGCCAAACGCAACCCAATGGACTTCGTTCTGTGGACGATGTCGAAAGCGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACGTTGAATGt                                                 >  1:161432/1‑251 (MQ=255)
                                                                                                                                                                                           aTCGCCAAAGCTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGACAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGGTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGCTGTCGCAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACCTTGAATGTTCGGCAAt                                         >  2:268380/1‑251 (MQ=255)
                                                                                                                                                                                            tctccAAAGGTCACGCTTATGTGGCGGCCAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAAGGGACCTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGATGGCCGTCTCCGCGGGGCGCGGGTGGTCCTGGCTGGCACATTGAATGTTCGGCAATg                                        >  2:15512/4‑251 (MQ=255)
                                                                                                                                                                                                cAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACt                                    <  1:124222/251‑1 (MQ=255)
                                                                                                                                                                                                   aGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCa                                 <  2:219808/251‑1 (MQ=255)
                                                                                                                                                                                                           cTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTg                         >  1:189374/1‑251 (MQ=255)
                                                                                                                                                                                                               tgtgGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGGCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTTGGCAATGAACTGCAAGCAGCTggga                      >  2:210067/1‑249 (MQ=255)
                                                                                                                                                                                                                gtgGCGGACAACGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGC‑ATCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAAc                   >  1:228853/1‑251 (MQ=255)
                                                                                                                                                                                                                          aCGGCGCCGTGATGTTCGACGTCCCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGat          >  1:79548/1‑251 (MQ=255)
                                                                                                                                                                                                                          aCGGCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGat          <  2:97057/251‑1 (MQ=255)
                                                                                                                                                                                                                            ggCGACGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAAAGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGAAAGCAGCTGGGGAACCACTTTGatat        >  2:95100/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cGTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCGGGCGCGCGCGTTGACGTGGTCGACGACAACCGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCAcg   <  2:235707/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gTGATGTTCGACGTCTCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCAcgg  <  1:208718/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           
CATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACCATATCGCAGAAATTATTGAACTCACTGAACAACTGATCGCCAAAGGTCACGCTTATGTGGCGGACAACGGCGACGTGATGTTCGACGTCCCGACCGATCCAACTTATGGCGTGCTGTCGCGTCAGGATCTCGACCAGCTGCAGGCAGGCGCGCGCGTTGACGTGGTCGACGACAAACGCAACCCAATGGACTTCGTTCTGTGGAAGATGTCGAAAGAGGGCGAACCGAGCTGGCCGTCTCCGTGGGGCGCGGGTCGTCCTGGCTGGCACATTGAATGTTCGGCAATGAACTGCAAGCAGCTGGGTAACCACTTTGATATCCACGG  >  NC_000913/554808‑555284

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