Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,657,901 A→T 7.0% intergenic (+31/‑168) ynfD → / → ynfE DUF1161 domain‑containing protein YnfD/putative selenate reductase YnfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,657,9010AT 7.0% 311.7 / 8.6 99intergenic (+31/‑168)ynfD/ynfEDUF1161 domain‑containing protein YnfD/putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base A (41/51);  new base T (4/3);  total (45/54)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGA  >  NC_000913/1657752‑1658048
                                                                                                                                                     |                                                                                                                                                   
aGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACaa                                                                                                                                                     <  1:1181125/150‑1 (MQ=255)
  gtcgtcGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAAt                                                                                                                                                   >  1:825929/1‑150 (MQ=255)
   tcgtcgGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATc                                                                                                                                                  <  1:547176/150‑1 (MQ=255)
    cgtcgGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                                                                                                                 <  2:1390626/150‑1 (MQ=255)
     gtcgGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCg                                                                                                                                                >  1:419323/1‑150 (MQ=255)
     gtcgGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg                                                                                                                                                >  1:862283/1‑150 (MQ=255)
         gCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                                                                            >  2:749124/1‑149 (MQ=255)
         gCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                                                                            <  1:1488161/150‑2 (MQ=255)
          ccATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAt                                                                                                                                           >  2:1074932/1‑150 (MQ=255)
            aTTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCTACTGAACGCCGGATATGACAAATCCGGCGAttt                                                                                                                                         <  2:1210948/150‑1 (MQ=255)
             ttGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTg                                                                                                                                        <  1:1223949/150‑1 (MQ=255)
              tGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcagt                                                                                                                                       <  2:419323/150‑7 (MQ=255)
              tGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGa                                                                                                                                       <  2:1190290/150‑1 (MQ=255)
              tGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGa                                                                                                                                       >  1:502676/1‑150 (MQ=255)
                cgcTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc                                                                                                                                     >  2:1602056/1‑150 (MQ=255)
                cgcTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc                                                                                                                                     <  2:590337/150‑1 (MQ=255)
                cgcTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc                                                                                                                                     <  1:706499/150‑1 (MQ=255)
                 gcTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                                                                                                    <  1:461768/150‑1 (MQ=255)
                 gcTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                                                                                                    >  2:308364/1‑150 (MQ=255)
                    aaTGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                                                                                                                 <  2:1592421/150‑1 (MQ=255)
                    aaTGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                                                                                                                 >  2:570797/1‑150 (MQ=255)
                    aaTGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                                                                                                                 <  1:1135118/150‑1 (MQ=255)
                     aTGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCACGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca                                                                                                                                <  2:829766/150‑1 (MQ=255)
                       gATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaac                                                                                                                              >  1:640172/1‑150 (MQ=255)
                        aTACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                                                                                                             >  1:585969/1‑150 (MQ=255)
                         tACGCATTAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATGATCGACTGAACGCCGGAGATGCCCACTCCGGCGATTTGAACATACAACAt                                                                                                                            >  2:763225/1‑150 (MQ=255)
                         tACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACAt                                                                                                                            >  2:307471/1‑150 (MQ=255)
                              aTAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATcc                                                                                                                       <  1:1126228/150‑1 (MQ=255)
                               tAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATccc                                                                                                                      >  1:690555/1‑150 (MQ=255)
                                aaaaTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCa                                                                                                                     >  2:1207806/1‑150 (MQ=255)
                                 aaaTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGa                                                                                                                                       >  1:441539/1‑131 (MQ=255)
                                 aaaTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGa                                                                                                                                       <  2:441539/131‑1 (MQ=255)
                                   aTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttc                                                                                                                                          <  2:1096164/126‑4 (MQ=255)
                                   aTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttc                                                                                                                                          >  1:1096164/1‑123 (MQ=255)
                                     tctcTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCtta                                                                                                                <  1:570797/150‑1 (MQ=255)
                                       tctATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                                                                                                              >  1:122552/1‑150 (MQ=255)
                                          ataCCCGCACCACCAGCGGTAACGTCTCAGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt                                                                                                           >  2:1105475/1‑150 (MQ=255)
                                           taCCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTc                                                                                                          >  1:176185/1‑150 (MQ=255)
                                            aCCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCa                                                                                                         <  2:122552/150‑1 (MQ=255)
                                             cccGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCAt                                                                                                        >  2:682023/1‑150 (MQ=255)
                                               cGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAc                                                                                                      >  1:1112063/1‑150 (MQ=255)
                                                        aGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt                                                                                             >  1:1200882/1‑150 (MQ=255)
                                                         gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTg                                                                                            <  1:612223/150‑1 (MQ=255)
                                                           ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat                                                                                          >  1:1397550/1‑150 (MQ=255)
                                                           ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat                                                                                          >  1:1572883/1‑150 (MQ=255)
                                                            gTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACACATACCCTTCTATTGata                                                                                         >  2:1334924/1‑150 (MQ=255)
                                                                cGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGa                                                                                     >  1:119070/1‑150 (MQ=255)
                                                                 gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                                    <  1:1074932/150‑1 (MQ=255)
                                                                 gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                                    <  1:1190304/150‑1 (MQ=255)
                                                                  tctcTGCTAACGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATAAGACAAATCCGGCGATATGAACCAACAACATAAACCCACCTTATTACTCATGCACATCTATTGATATGGAtt                                                                                   >  2:1170413/1‑150 (MQ=255)
                                                                   ctctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATta                                                                                  >  2:575535/1‑150 (MQ=255)
                                                                     ctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaat                                                                                <  1:1151050/150‑1 (MQ=255)
                                                                       gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGCACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATGCCACCTTATTACTCATACGCTTCTATTGATATGGATtaataa                                                                              >  1:173805/1‑150 (MQ=255)
                                                                          cccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGCACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAAtta                                                                           >  2:1060004/1‑149 (MQ=255)
                                                                           ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCt                                                                          <  1:682023/150‑1 (MQ=255)
                                                                            cgcgCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAtt                                                                                   <  1:311193/140‑1 (MQ=255)
                                                                            cgcgCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAtt                                                                                   >  2:311193/1‑140 (MQ=255)
                                                                             gcgcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTa                                                                        >  2:811220/1‑150 (MQ=255)
                                                                              cgcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                                                                       <  2:640172/150‑1 (MQ=255)
                                                                                cAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAcc                                                                     <  2:579440/150‑1 (MQ=255)
                                                                                 aGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAccc                                                                    <  2:464906/150‑1 (MQ=255)
                                                                                 aGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAccc                                                                    <  1:1464511/150‑1 (MQ=255)
                                                                                    cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaa                                                                 <  1:871792/150‑1 (MQ=255)
                                                                                         cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATggg                                                            >  2:1335455/1‑150 (MQ=255)
                                                                                         cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATggg                                                            <  2:1338737/150‑1 (MQ=255)
                                                                                          aGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGt                                                           <  1:1263298/150‑1 (MQ=255)
                                                                                           ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTa                                                          >  1:1521733/1‑150 (MQ=255)
                                                                                            gATGGCTCGCCTGCCGAACCTCATTAACACATTATCTACTGAACGCCGTATATGACAAATCCGGCGATTTTAACATACAACATAATCCCACATTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATgtttag                                                         <  1:334881/150‑6 (MQ=255)
                                                                                            gATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAg                                                         <  1:334879/150‑1 (MQ=255)
                                                                                             aTGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGa                                                        <  2:624197/150‑1 (MQ=255)
                                                                                                  gcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTccc                                                   <  2:15767/150‑1 (MQ=255)
                                                                                                    gcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCctc                                                 <  1:307471/150‑1 (MQ=255)
                                                                                                     ccTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCctct                                                <  2:825929/150‑1 (MQ=255)
                                                                                                        tccGAAACGCAGTAAAACATTATCGTCTGAACGCCGGATATTACAAATCCGTCGATTTGAACATACAACATAATCCCAACTTATTAATCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATTGGTAGACTCCCTCTAtt                                             <  2:167213/149‑1 (MQ=255)
                                                                                                        gCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTAtt                                             <  2:532784/150‑1 (MQ=255)
                                                                                                         ccGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTg                                            <  1:1602056/150‑1 (MQ=255)
                                                                                                           gAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGtt                                          <  2:1004122/150‑1 (MQ=255)
                                                                                                           gAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACACATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGtt                                          <  1:1334924/150‑1 (MQ=255)
                                                                                                             aCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAg                                        <  2:821419/150‑1 (MQ=255)
                                                                                                               cGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAgcg                                      >  2:1523021/1‑150 (MQ=255)
                                                                                                               cGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAgcg                                      <  2:1200882/150‑1 (MQ=255)
                                                                                                                gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAATTGAACATACAACATAATCCCACCTTATTACTAATACACTTCTATTGATATGGATTAATAATTCTAAACACAAAATGGGTAGACTCCCTCTATTGTTAgccc                                     >  2:1593197/1‑148 (MQ=255)
                                                                                                                   gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCt                                  <  2:862283/150‑1 (MQ=255)
                                                                                                                   gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCt                                  >  1:1169878/1‑150 (MQ=255)
                                                                                                                     aaCACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTaa                                >  2:52477/1‑150 (MQ=255)
                                                                                                                        acaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAata                             <  1:579620/150‑1 (MQ=255)
                                                                                                                        acaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATACCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGTTAGACTCCCTATATTCTTAGCGCGCTAAata                             >  1:721407/1‑150 (MQ=255)
                                                                                                                           ttATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTc                          <  1:1105475/150‑1 (MQ=255)
                                                                                                                             aTCGACTGAATGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCaa                        <  1:1007109/150‑1 (MQ=255)
                                                                                                                              tCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAat                       <  2:1113043/150‑1 (MQ=255)
                                                                                                                              tCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAat                       <  2:268439/150‑1 (MQ=255)
                                                                                                                                gACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAatat                     >  1:1052622/1‑150 (MQ=255)
                                                                                                                                    gAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATaa                 >  1:853794/1‑150 (MQ=255)
                                                                                                                                     atcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATaaa                >  1:1102236/3‑150 (MQ=255)
                                                                                                                                       cGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACt              >  1:1191994/1‑150 (MQ=255)
                                                                                                                                          cGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACATTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACtttt           >  2:786081/1‑150 (MQ=255)
                                                                                                                                           ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTta          <  2:1397550/150‑1 (MQ=255)
                                                                                                                                             atatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTtata        <  2:1521733/150‑1 (MQ=255)
                                                                                                                                                 gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAAc    <  1:52477/150‑1 (MQ=255)
                                                                                                                                                   cAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGa  <  2:1627372/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGA  >  NC_000913/1657752‑1658048

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