Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,722,946 G→A A74T (GCT→ACT)  pssA → phosphatidylserine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,722,9460GA95.8% 129.4 / ‑4.7 48A74T (GCT→ACT) pssAphosphatidylserine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/23);  minor base C (2/0);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

AATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCC  >  NC_000913/2722701‑2723185
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
aaTGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTaa                                                                                                                                                                                                                                             >  1:21196/1‑250 (MQ=255)
    aaCAGCGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGc                                                                                                                                                                                                                                        >  1:95012/1‑251 (MQ=255)
                   cACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTgg                                                                                                                                                                                                                         <  1:458721/251‑1 (MQ=255)
                             ttgtCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCgggg                                                                                                                                                                                                               >  1:466650/1‑250 (MQ=255)
                                         aaGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAACAGGCAGCGTCCGGAACTGGATGTGAGGGTGCTGGTCGACTg                                                                                                                                                                                                   >  2:247442/1‑251 (MQ=255)
                                          aGCGAAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCATTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTgg                                                                                                                                                                                                  <  2:212886/251‑1 (MQ=255)
                                             gTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATAAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAAGCAGCGTCCGGCACTGGCTGTGCAGGTGCTGGTCGACTGGCAt                                                                                                                                                                                               >  2:416303/1‑251 (MQ=255)
                                              taataaACATCAACAAAACCTTACCAAACTACACAAGCTTTCTCAATCAGTTGATGATGTCGCTTTCTTTTACGCTCCCGCCGACTTCCAGGCGACGCTGCTGGAAAAAATAGACAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATATCGAACAGGATGACGGTGGAAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGAGCGGGTGCTGGTCGATTGGCATc                                                                                                                                                                                              >  1:149487/1‑251 (MQ=255)
                                                      aTCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTATGCACGCGTTGTATGAGCCTAAAAGGCAGCGTCCGGAACTGGCTGTGAGGGTGCTGGTCGACTGGCATCGTGccaa                                                                                                                                                                                       >  2:308267/1‑247 (MQ=255)
                                                         aacaacACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGCTTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGCCTGGCATCGTGCACAACgg                                                                                                                                                                                   >  2:185057/1‑250 (MQ=255)
                                                          acaacaCCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTg                                                                                                                                                                                  <  1:276664/251‑1 (MQ=255)
                                                                ccTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGc                                                                                                                                                                             >  2:113852/1‑250 (MQ=255)
                                                                   tGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCAtt                                                                                                                                                                          <  1:113852/250‑1 (MQ=255)
                                                                    gCCCAACTATCCAAGTTTACTGAATCAGGTGATGTTTTCGATTTCTTTTCCGCTTCCGCCTACTTCCGGGAGCAGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGTATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCCTTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGTATGTGCGGGTGTTGGTCGACTGGCATCGTGCACAATGTGGACGCATTg                                                                                                                                                                         <  2:284286/250‑1 (MQ=255)
                                                                          cTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGc                                                                                                                                                                  >  1:145985/1‑251 (MQ=255)
                                                                                          aaTCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAAc                                                                                                                                                  <  1:391847/251‑1 (MQ=255)
                                                                                           agcAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTg                                                                                                                                                                   <  2:413793/231‑1 (MQ=255)
                                                                                           aTCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTg                                                                                                                                                                   >  1:413793/1‑233 (MQ=255)
                                                                                               gtggaggatgtcgatgtcTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGa                                                                                                                                             <  2:466650/245‑1 (MQ=255)
                                                                                                       tgTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCAGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGtcct                                                                                                                                     >  2:110553/1‑248 (MQ=255)
                                                                                                                tttttACGCCCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATg                                                                                                                             <  2:268815/249‑1 (MQ=255)
                                                                                                                   ttACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGgcgc                                                                                                                         >  1:138584/1‑251 (MQ=255)
                                                                                                                   ttACGCTACCGCCGACTTCCGGGAGACGCTGCTGGACAAAATAAACAGCTCGCAGAAGCGCCTTTGCATTGTCGCACTGTATCTAGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAACGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGCAGCATGGGCGCTGCGGCATCTAACACTAAAGCTGACTGGTACTGCCGCATGgcgc                                                                                                                         >  2:477554/1‑251 (MQ=255)
                                                                                                                           ccgccgACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAAt                                                                                                                 <  1:108861/251‑1 (MQ=255)
                                                                                                                            cgccgACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATc                                                                                                                <  1:357891/251‑1 (MQ=255)
                                                                                                                             gccgACTTCCGGTGGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCCTTTGCATTTTCGCCCTGTATCTCGAACATGATGACGGTGGCAAAGGCATTGTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACATCGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGATAATcc                                                                                                               <  2:422013/251‑1 (MQ=255)
                                                                                                                               cgACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTACCACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCgg                                                                                                             >  2:51077/1‑251 (MQ=255)
                                                                                                                                                                   cgcgAAGCAGCGCGTTTGAATCGTCGCCCTGTATCACGAACAGGCAGACGGTGGCAACGGCATTCTGAACGCGTTGTATGAGCCTAAAAGGAAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGACATCGTGCACAACGTGGACGCATTGGCGCGGCGGCATCTAACACTAACGCGGAATGGTAATGCCGCCGGGCACAGGAACATCCGGGCATAGATGTTCCGGTTTAGGGCGTTCcaataaat                                                                          >  2:448511/1‑250 (MQ=255)
                                                                                                                                                                     cgAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGGCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACt                                                                       <  1:479065/251‑1 (MQ=255)
                                                                                                                                                                     cgAAGCAGCGCATTTGCATTATCGCCCTGTATCTCGCAAAGGATGACGGTGGCAATGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGAAACTGGCTGTGCGGGTGATGGTCGACTGGCATCATGCACAACGTGGACGCATTGGCGCTGCGGCATCTACAACTAACGCTGACTGGTACTGCCGCATGGCGCCGAAAAATCCGGGCGGAGAGGTTCCGGTTTATGGCGTTCCAATCAATaat                                                                       >  2:257070/1‑249 (MQ=255)
                                                                                                                                                                        agcagcGCATTTGCAGTGTTGCCCTGCATCTCCAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGGCCGGAACTGGATGTGCTGGTGCGGGTCGACTGGAATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACTCTATCGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGATCCGGTTTATGGCGTTCCAATCAATACTCGt                                                                    <  2:440623/251‑1 (MQ=255)
                                                                                                                                                                         gcatcgcATTTTCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACTCTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTCATGGCGTGCCAAGCAATACTCGTg                                                                   <  2:241959/247‑1 (MQ=255)
                                                                                                                                                                          cagcGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTCACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGATTATGGCGTTCCAATCAATACTCGTg                                                                   >  1:287869/1‑250 (MQ=255)
                                                                                                                                                                                      aTTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGgtgt                                                      <  2:265793/251‑1 (MQ=255)
                                                                                                                                                                                         gTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCt                                                   <  1:51077/251‑1 (MQ=255)
                                                                                                                                                                                              ccTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACt                                              >  1:226040/1‑251 (MQ=255)
                                                                                                                                                                                                              gATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTatcat                              >  1:168126/1‑251 (MQ=255)
                                                                                                                                                                                                               aTGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCTACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAAc                                                                                                                                                  <  2:62504/134‑1 (MQ=255)
                                                                                                                                                                                                               aTGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAAc                                                                                                                                                  >  1:62504/1‑134 (MQ=255)
                                                                                                                                                                                                                  aCGGTGGCAAAGGCATTCTGAACGCGTGGTATGAGACGAAAAGGCAGCGTCCGGGACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCGCTTGGTGTTCTGCACTTTAAAGGCTTTATCATcgac                          <  2:260725/251‑1 (MQ=255)
                                                                                                                                                                                                                       ggCAAAGGCATTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGAt                       >  1:411200/1‑249 (MQ=255)
                                                                                                                                                                                                                        gCAATGGCATTCTGAACTCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCTCAACGTGGACGCAGTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCTCTTTAAAGGCTTTATCATCGACGATAGc                    <  2:408945/251‑1 (MQ=255)
                                                                                                                                                                                                                            aGGCATTCTGAACGCGTTGTATGAGACGAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGGCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTAc                <  2:168126/251‑1 (MQ=255)
                                                                                                                                                                                                                                aTTCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACGGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTta            <  2:21196/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ttCTGAACGCGTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTACACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGACCTTTAAAGGCTTTATCATCGACGATAGAGTACTTtat           >  2:289858/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGCAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAACGGCTTTat                                 >  2:137230/1‑219 (MQ=255)
                                                                                                                                                                                                                                           gTTGTATGAGACTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTat                                 <  1:137230/219‑1 (MQ=255)
                                                                                                                                                                                                                                            ttGTATGAGACTAAAAGGCAGCGTCAGGAACTGGATGTGCGGGTGCTGGTAGACTGGCATCGCGCACAACGTGGCCGCATTGGCGCTGAGGCATCTAACACTAACGCTGACTGGTACTGCAGCATGGCGCAGGAAAATCCGGGCGTAGGTGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTCCTTTATAGCGGTGcc  >  2:193967/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
AATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACATCAACAACACCTTGCCCAACTACCCAAGATTTCTCAATCAGTTGATGATGTCGATTTCTTTTACGCTCCCGCCGACTTCCGGGAGACGCTGCTGGAAAAAATAGCCAGCGCGAAGCAGCGCATTTGCATTGTCGCCCTGTATCTCGAACAGGATGACGGTGGCAAAGGCATTCTGAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCC  >  NC_000913/2722701‑2723185

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