Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,474,846 A→G G40G (GGT→GGC tusB ← sulfurtransferase complex subunit TusB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,474,8460AG98.4% 195.9 / NA 63G40G (GGT→GGCtusBsulfurtransferase complex subunit TusB
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (26/36);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAG  >  NC_000913/3474602‑3475086
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
cccTCACACAATATGAGGACGCCGAATTTTAGGGCGATTCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTAGAGTCAATGAGAATGATGTCGTGCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATATACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCa                                                                                                                                                                                                                                            <  1:393251/250‑1 (MQ=255)
         aaTATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGTGGGCATTACGCAGACTTTCAAGGGAGCGGTTGCCGTCAACTGCGGCa                                                                                                                                                                                                                                   <  2:369541/251‑1 (MQ=255)
          atatGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGc                                                                                                                                                                                                                                    >  1:81236/1‑249 (MQ=255)
              gAGGACGCCTAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCGTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTAc                                                                                                                                                                                                                              <  2:336154/251‑1 (MQ=255)
              gAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTa                                                                                                                                                                                                                               <  1:102300/250‑1 (MQ=255)
              gAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTAc                                                                                                                                                                                                                              <  1:199252/251‑1 (MQ=255)
                   cGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCAt                                                                                                                                                                                                                         <  1:132444/251‑1 (MQ=255)
                      cGAATTTTAGCGCGATGCCGCCAAGGCGTCACGACATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCACTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTTATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCtt                                                                                                                                                                                                                      >  2:192583/1‑251 (MQ=255)
                       gAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTg                                                                                                                                                                                                                     >  2:220210/1‑251 (MQ=255)
                         aTTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACAAGGCCATCTGGCTGGGGTGCTTCACCGTAAGTCTGCCGAAATCAGTATAGACAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGCATTTCAAGGTAGCGGTTGCCGTCAACTGCGGCCGTTACGCCATCTTGc                                                                                                                                                                                                                    >  1:200537/1‑250 (MQ=255)
                             tAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCaata                                                                                                                                                                                                                <  1:332870/250‑1 (MQ=255)
                                       ccGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCg                                                                                                                                                                                                     <  1:321260/251‑1 (MQ=255)
                                        agAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCATGGCATAGACCTTAATGGGGGCATTACGCAGCCTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTTCGCCATCTTGCAATAATAGCAGTTCGt                                                                                                                                                                                                    <  2:352908/250‑1 (MQ=255)
                                             atgtgtCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCt                                                                                                                                                                                               >  1:180586/3‑251 (MQ=255)
                                             aGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCt                                                                                                                                                                                               <  1:157427/251‑1 (MQ=255)
                                                  tCAAGAAATATACAACGATCACGACATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCGTAGACCTTAATGGGGGCATTAAGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACt                                                                                                                                                                                          >  2:354175/1‑251 (MQ=255)
                                                   cAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACt                                                                                                                                                                                          <  1:105853/250‑1 (MQ=255)
                                                        aaTATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGCAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGc                                                                                                                                                                                      >  2:287741/1‑249 (MQ=255)
                                                             aCAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACgc                                                                                                                                                                                >  2:137505/1‑250 (MQ=255)
                                                                              aCCAGGCCATCTGGCTGGGGTGCTTAACCGTACGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGCCCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCACTAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATcc                                                                                                                                                              >  1:108999/1‑251 (MQ=255)
                                                                                cAGGTCATCTGGCTGGGGTGCTTAACCTGAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACGGCGGCAGTTACGCCATCTTGCAATAAGAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAAGCCGt                                                                                                                                                            <  2:348165/251‑1 (MQ=255)
                                                                                   gCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAg                                                                                                                                                         <  1:190647/251‑1 (MQ=255)
                                                                                   gCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCGCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAg                                                                                                                                                         <  2:348742/251‑1 (MQ=255)
                                                                                    ccATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGc                                                                                                                                                        >  1:414698/1‑251 (MQ=255)
                                                                                              ggggTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGGCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCAGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGa                                                                                                                                                <  2:317623/249‑1 (MQ=255)
                                                                                              ggggTGCTTAACCGTAAGTCGGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCGCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGa                                                                                                                                                <  1:391039/249‑1 (MQ=255)
                                                                                                   gCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATg                                                                                                                                          <  1:197727/250‑1 (MQ=255)
                                                                                                     ttAACCGTACGTCTGAAGACATCAGTATAGTCAATGAGACTGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTTAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTa                                                                                                                                        >  2:298581/1‑250 (MQ=255)
                                                                                                        aCCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAAtgt                                                                                                                                    >  2:48427/1‑251 (MQ=255)
                                                                                                         ccGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAAtgtg                                                                                                                                   <  2:81236/251‑1 (MQ=255)
                                                                                                           gTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAAtgtgtg                                                                                                                                 >  1:54737/1‑251 (MQ=255)
                                                                                                                           cAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGc                                                                                                                                              >  1:176288/1‑222 (MQ=255)
                                                                                                                           cAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAAGCCGTCAGCCAGTGTGAGc                                                                                                                                              <  2:176288/222‑1 (MQ=255)
                                                                                                                            aGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGGCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGcctcaa                                                                                                                <  1:256019/251‑1 (MQ=255)
                                                                                                                                 aggcAATGAGAATGAGGTCGTTCGAAATTTGACCAGGCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTGGCAATAATAGCAGTTCGTCGCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGcctcaaaacc                                                                                                            <  2:108999/247‑1 (MQ=255)
                                                                                                                                  gTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGcctcaaaacctc                                                                                                          <  1:250239/251‑1 (MQ=255)
                                                                                                                                      aTGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCg                                                                                                                                             >  1:320029/1‑212 (MQ=255)
                                                                                                                                      aTGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCg                                                                                                                                             <  2:320029/212‑1 (MQ=255)
                                                                                                                                       tGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTCATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCAGCCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAAt                                                                                                     >  1:26336/1‑251 (MQ=255)
                                                                                                                                         agaATGATGTCGTTCGAACTTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATa                                                                                                    >  2:273615/1‑250 (MQ=255)
                                                                                                                                             agatgTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATc                                                                                               <  1:192679/250‑1 (MQ=255)
                                                                                                                                                   cGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGtt                                                                                         >  2:393626/1‑251 (MQ=255)
                                                                                                                                                    gTTCGAATTTTGACCAGACAAACCGCGGGCAATAAGGTCTTCGTTCAGAGCATAGACCTGAATGGGGGCATTACGCAGACTTTCAAGGCAGCGGTTGCCGTCAACTGCGGCAGTTACACCATCTTGCAATAATAGCAGTTCGCCTCCTTCACTGAGCAGACTCATCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGAGGTAATGTGCGCAGCATAAACGCCTGAACACCTCAAAATAATATCGTAGTTg                                                                                        >  1:302587/1‑251 (MQ=255)
                                                                                                                                                           aTTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGAAGACTTTCCAGGAAGAGGTTACAGTCAACTGAGGCAGTTACGCCATCTTGCAATAATAGCAGCTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGCGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAACCCTCAAAATAACATCGTAGTTGGCGAGtt                                                                                 >  1:212179/1‑251 (MQ=255)
                                                                                                                                                           aTTTGACCAGTCAAACCGCGGGCAATAAGGTCGTCGTTCAGGGCCGAGACCTTAATGGGGTCATTACTCAGACTTTCAAGGTAGTGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCTGCGCAGCAAAATCCGTCAGCCAGGTTGAGCGATGTAATGTCTGCAGTATAAACGTCTGAAAACTTCAAGATATCATCGTAGTTGGCTAGtt                                                                                 <  2:371239/251‑1 (MQ=255)
                                                                                                                                                                  cAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCACAACCTCAAACTAACATCGTAGTTGGCGAGTTCgcggcg                                                                          >  2:86313/1‑251 (MQ=255)
                                                                                                                                                                     tCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTa                                                                        >  1:151821/1‑250 (MQ=255)
                                                                                                                                                                        aaCCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGc                                                                    <  1:286739/251‑1 (MQ=255)
                                                                                                                                                                        aaCCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGc                                                                    <  1:220210/251‑1 (MQ=255)
                                                                                                                                                                               ggCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTGCGCGGCGTAAGGCATCTGCt                                                             <  1:124759/251‑1 (MQ=255)
                                                                                                                                                                                 caagaaGGTCTTTGTTCAGGGCATAGACCTGATTGGTGGCATTACGCAGACTTTCAAGGTAGCTGTTGCCGTCAACTGCGGCAGTTACACCATCTTGCAATAATAGCAGTTCGTCTCCGTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCTTCAGCCAGGGTGATCGATGTAATGGGTGCAGCATAAACGCCTCAAAACCTCAAAAGATCATCGTAGTGGTCGAGTTCGCGTCGTAAGGCATCTGCTTc                                                           <  2:200537/247‑1 (MQ=255)
                                                                                                                                                                                 cAATAAGGTCTTCGTGCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTAAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCCGTTCGACTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGTCGAGTTCGCGGCGTAAGGCATCTGCTTc                                                           >  2:236701/1‑251 (MQ=255)
                                                                                                                                                                                               ttCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCtt                                             >  1:323826/1‑251 (MQ=255)
                                                                                                                                                                                                     gCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACTCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCTAGCGGCGTGGCTTCGacaa                                       <  1:140311/251‑1 (MQ=255)
                                                                                                                                                                                                     gCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGaca                                        >  1:205456/1‑250 (MQ=255)
                                                                                                                                                                                                            ccTTAATGGGGGCAGTACGCAGACTGTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCATCGCAGCAAAATCCGGCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCTGCGTGGCTTCGACATCAAAggg                                <  2:37074/251‑1 (MQ=255)
                                                                                                                                                                                                                      ggCATTACGCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGa                      <  1:209787/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAAccc              >  2:34201/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCAGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAAccc              <  1:230548/251‑1 (MQ=255)
                                                                                                                                                                                                                                aGACTTTCAAGGTAGCGGTTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATCGCCGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGGCTTCGGCTCTAACCcgc            >  2:190108/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aCTTTCAAGGTAGCGGGTGCCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGATCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGt          <  2:164218/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ttCAAGGTAGCGGTTGCCGTCAACTGCGGCAGGTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGtt         <  2:43715/249‑1 (MQ=255)
                                                                                                                                                                                                                                          ggTAGCGGTTGCCGTCAACAGCGGCAGTTAAGCCATCTCGCAATAATAGCAGTTCGTCTCCTTCACTGATCACACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGGGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTATTTGGCGAGTTCGCGGCGTAATGCAGCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCAGATCTGCCCCGCGTTcgcgtag  <  2:184721/251‑4 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
TCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATACAACGATCCCGCCATCACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAG  >  NC_000913/3474602‑3475086

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