Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,151,954 T→C V265A (GTG→GCG)  unknown → Molybdenum transport system protein ModD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,151,9540TC97.9% 129.8 / ‑6.3 48V265A (GTG→GCG) unknownMolybdenum transport system protein ModD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/28);  minor base G (0/1);  total (19/29)
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.96e-01

AGACTGGAAGCTCCGGAAAATAAAATCACCTTGGAAGCTGATAATTTTTCGCAATTTGAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGTGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCACTTTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTCCCCTCCACAATTCAACAAAATAGAATTTATACCCTTAAAT  >  1/1151710‑1152178
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                
agaCGGGAAGCTCCGGAAAATAAAATCACCTTGGAAGCTGATAATTTTTCGCAATTTGAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCtt                                                                                                                                                                                                                             <  1:28469/250‑1 (MQ=255)
       aaGCTCCGGAAAATAAAATCACCTTGGAAGCTGATAATTTTTCGCAATTTGAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGc                                                                                                                                                                                                                      <  1:500883/250‑1 (MQ=255)
         gCTCCGGAAAATAAAATCACCTTGGAAGCTGATAATTTTTCGCAATTTGAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGGAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGGGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTAc                                                                                                                                                                                                                           >  1:737065/1‑243 (MQ=255)
                 aaataaaatCACCTTGGAAGCTGATAATTTTTCGCAATTTGAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAg                                                                                                                                                                                                           <  1:109085/251‑1 (MQ=255)
                 aaataaaatCACCTTGGAAGCTGATAATTTTTCGCAATTTGAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAg                                                                                                                                                                                                           >  1:874729/1‑251 (MQ=255)
                        atCACCTTGGAAGCTGATAATTTTTCGCAATTTGAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGGCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTaa                                                                                                                                                                                                    <  1:437478/251‑1 (MQ=255)
                          cACCTTGGAAGCTGATAATTTTTCGCAATTTGAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAg                                                                                                                                                                                                  <  1:398428/251‑1 (MQ=255)
                                     cTGATAATTTTTCGCAATTTGAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGGATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAg                                                                                                                                                                                       <  1:759424/251‑1 (MQ=255)
                                                         gAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGGGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAg                                                                                                                                                                                                  >  1:768620/1‑220 (MQ=255)
                                                              aaTGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCAccgc                                                                                                                                                                                                                 <  1:583865/200‑1 (MQ=255)
                                                                 gCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTaa                                                                                                                                                                                                    >  1:13373/1‑210 (MQ=255)
                                                                  cTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTaaaaaaa                                                                                                                                                          <  1:761078/251‑1 (MQ=255)
                                                                        gCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTaaaaaaa                                                                                                                                                          >  1:160742/1‑245 (MQ=255)
                                                                            aaCCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCACCTTACTATGCACCGCCGGAAGATATTAAAgtgg                                                                                                                                                                                               >  1:431550/1‑204 (MQ=255)
                                                                                ggATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAg                                                                                                                                                                                       <  1:870808/208‑1 (MQ=255)
                                                                                  atatTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCACttttttt                                                                                                                                          >  1:73125/1‑251 (MQ=255)
                                                                                  atatTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCACttttttt                                                                                                                                          >  1:417628/1‑251 (MQ=255)
                                                                                   tattatTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCACtttttttt                                                                                                                                         <  1:387208/251‑1 (MQ=255)
                                                                                     ttattCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAgt                                                                                                                                                                                                 <  1:243480/193‑1 (MQ=255)
                                                                                       attCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCACttttttttaaat                                                                                                                                     >  1:546478/1‑247 (MQ=255)
                                                                                            gTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTa                                                                                                                               >  1:709077/1‑251 (MQ=255)
                                                                                                    aaTGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAacac                                                                                                                       <  1:230957/251‑1 (MQ=255)
                                                                                                          cGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAac                                                                                                                                                                     >  1:500541/1‑200 (MQ=255)
                                                                                                                aaCAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAAt                                                                                                            <  1:281025/250‑1 (MQ=255)
                                                                                                                   aGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCtt                                                                                                        >  1:575854/1‑251 (MQ=255)
                                                                                                                    ggTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGc                                                                                                                                                                                                                      <  1:679231/141‑1 (MQ=255)
                                                                                                                    ggTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAg                                                                                                                                                                                                  >  1:205324/1‑161 (MQ=255)
                                                                                                                    ggTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGAGAATAAAAAAAACATAGCAGATTATGTCAAGTTGGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGGACCGCCGGAAg                                                                                                                                                                                                           >  1:399013/1‑152 (MQ=255)
                                                                                                                        cAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTaaaaa                                                                                                   >  1:678145/1‑251 (MQ=255)
                                                                                                                              gCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCAc                                                                                                                                                                                                                    <  1:413200/133‑1 (MQ=255)
                                                                                                                                   ggCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGa                                                                                                                                                                                                          >  1:425739/1‑138 (MQ=255)
                                                                                                                                                                                   cggGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTcc                                         <  1:501845/250‑1 (MQ=255)
                                                                                                                                                                                   cggGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTcc                                         <  1:174574/250‑1 (MQ=255)
                                                                                                                                                                                            tAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTCCCCTCCACAAt                               <  1:844050/251‑1 (MQ=255)
                                                                                                                                                                                                 acaacaTAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTCCCCTCCACAATTcaac                          <  1:828359/251‑1 (MQ=255)
                                                                                                                                                                                                    acaTAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTAt                                                                                                                                                                                                                        <  1:876896/59‑1 (MQ=255)
                                                                                                                                                                                                    acaTAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAAc                                                                                                                                                                           >  1:130566/1‑104 (MQ=255)
                                                                                                                                                                                                    acaTAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTc                                                                                                          <  1:428885/168‑1 (MQ=255)
                                                                                                                                                                                                    acaTAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTc                                          <  1:174160/232‑1 (MQ=255)
                                                                                                                                                                                                       tAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAac                                                                                                                                                                     <  1:628469/107‑1 (MQ=255)
                                                                                                                                                                                                         gCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTcccc                                       <  1:422546/230‑1 (MQ=255)
                                                                                                                                                                                                         gCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTCCCCTCCACAATTCAACAAAATAGa                  >  1:278401/1‑251 (MQ=255)
                                                                                                                                                                                                             aTTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTCCCCTCCACAATTCAACAAAATAGAAttt              <  1:214345/251‑1 (MQ=255)
                                                                                                                                                                                                                   ccAAGTTAGGGATAAATTTATTTATCATCTCCGCGCCTTACAATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTCCCCTCCACAATTCAACAAAATAGAATTTATAcc         <  1:644628/250‑1 (MQ=255)
                                                                                                                                                                                                                   ccAAGTTAGGGATAAATTTATTTATCATCTCCCGGCCTAACAAGGCCCCGCCGAAGGATATTAAATTGTTGA‑TGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTCCCCTCCACAATTCAACAAAATAGAATTTATAccc        <  1:699015/250‑1 (MQ=255)
                                                                                                                                                                                                                         tAGGGATAAATTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTCCCCTCCACAATTCAACAAAATAGAATTTATACCCTtaaat  <  1:57777/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aaTTTATTTATCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTa                                                                             >  1:784427/1‑168 (MQ=255)
                                                                                                                                                                                                                                         tatCATCTCCGCGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCAC‑TTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTCCCCTCCACaa                                <  1:885191/205‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                
AGACTGGAAGCTCCGGAAAATAAAATCACCTTGGAAGCTGATAATTTTTCGCAATTTGAACAAATGCTTGCAGCCGAACCGGATATTATTCAGTTGGATAAATGGTCGCCGGAACAGGTTCAACAGGCATTGGCTTTGTTAAAAAAGCAAAGTAAAAAAACGCTCCTTTCTGTTGCCGGCGGGGTGAATAAAAACAACATAGCAGATTATGCCAAGTTAGGGATAAATTTATTTATCATCTCCGTGCCTTACTATGCACCGCCGGAAGATATTAAAGTGGTGATTGAGCGGATGCAGAACAAAAACACTTAAAAAAATGACCGCACTTTTTTTTTAAATAGTTAAAAAACACAACGCAGAAATTCTTAAAAAAAAGTTATTAAAAGTGCGGTTAAAAAACAAATTAAAAATTCACCGTACTTTTTATTCCCCTCCACAATTCAACAAAATAGAATTTATACCCTTAAAT  >  1/1151710‑1152178

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