Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 944,551 C→T 34.9% P22L (CCT→CTT)  perR → transcriptional regulator (Fur family)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome944,5510CT34.9% 486.2 / inf 321P22L (CCT→CTT) perRtranscriptional regulator (Fur family)
Reads supporting (aligned to +/- strand):  ref base C (106/103);  new base T (56/56);  total (162/159)
Fisher's exact test for biased strand distribution p-value = 9.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

Only 100 of 331 total aligned reads displayed.

ATTGACAAATTATCGTAGAAAGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTCCGGGAATCAGGTTTGG  >  BS168_reference_genome/944403‑944700
                                                                                                                                                    |                                                                                                                                                     
aTTGACAAATTATCGTAGAAAGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTActcct                                                                                                                                                      >  1:3723704/1‑150 (MQ=255)
 ttGACAAATTATCGTAGAAAGAGATACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTActcctc                                                                                                                                                     <  1:4671701/150‑1 (MQ=255)
  tGACAAATTATCGTAGAAAGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCa                                                                                                                                                    >  2:3517430/1‑150 (MQ=255)
      aaaTTATCGTAGAAAGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGt                                                                                                                                                <  2:1900767/150‑1 (MQ=255)
       aaTTATCGTAGAAAGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTc                                                                                                                                               <  1:78917/150‑1 (MQ=255)
         ttATCGTAGAAAGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCAt                                                                                                                                             >  2:3924369/1‑150 (MQ=255)
          tATCGTAGAAAGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGGTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGGCATg                                                                                                                                            >  2:40978/1‑150 (MQ=255)
             cGTAGAAAGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGAAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGa                                                                                                                                         >  1:3876932/1‑150 (MQ=255)
                   aaGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTgg                                                                                                                                   >  1:2166207/1‑150 (MQ=255)
                   aaGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTgg                                                                                                                                   >  2:498054/1‑150 (MQ=255)
                      agTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAat                                                                                                                                <  1:2021421/150‑1 (MQ=255)
                         tACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATAtc                                                                                                                             <  2:4324061/150‑1 (MQ=255)
                         tACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATAtc                                                                                                                             >  2:418452/1‑150 (MQ=255)
                          acacTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATAtct                                                                                                                            >  2:2972744/1‑150 (MQ=255)
                           cacTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATAtctc                                                                                                                           <  2:2038407/150‑1 (MQ=255)
                           cacTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATAtctc                                                                                                                           <  1:1134163/150‑1 (MQ=255)
                              tAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGtt                                                                                                                        >  1:632264/1‑150 (MQ=255)
                               aaTTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTa                                                                                                                       >  2:4352381/1‑150 (MQ=255)
                                aTTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTaa                                                                                                                      >  1:906419/1‑150 (MQ=255)
                                  tataAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAAct                                                                                                                    <  2:3606991/150‑1 (MQ=255)
                                    taAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAActct                                                                                                                  <  1:2907016/150‑1 (MQ=255)
                                    taAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAActat                                                                                                                  <  1:2740773/150‑3 (MQ=255)
                                     aaaCATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTa                                                                                                                 <  2:106241/150‑1 (MQ=255)
                                      aaCATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTAt                                                                                                                >  1:7468/1‑150 (MQ=255)
                                      aaCATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTAt                                                                                                                >  2:4675660/1‑150 (MQ=255)
                                          ttACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCt                                                                                                            <  1:705444/150‑1 (MQ=255)
                                             cAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCAt                                                                                                         <  2:2460597/150‑1 (MQ=255)
                                             cAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCAt                                                                                                         <  1:4528470/150‑1 (MQ=255)
                                                tGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCa                                                                                                      >  1:3535055/1‑150 (MQ=255)
                                                  tAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAAc                                                                                                    >  2:3230104/1‑150 (MQ=255)
                                                    aGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCCACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAg                                                                                                  <  1:164833/150‑1 (MQ=255)
                                                     gAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAAGCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGc                                                                                                 <  1:3642793/150‑1 (MQ=255)
                                                     gAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGc                                                                                                 >  2:3713491/1‑150 (MQ=255)
                                                      aaTTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCg                                                                                                >  1:4819260/1‑150 (MQ=255)
                                                         tttttttAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGAc                                                                                             <  2:1368215/150‑1 (MQ=255)
                                                            ttttAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGat                                                                                          >  1:1752486/1‑150 (MQ=255)
                                                              ttAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGatat                                                                                        <  2:2255970/150‑1 (MQ=255)
                                                              ttAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGatat                                                                                        >  2:3736382/1‑150 (MQ=255)
                                                               tAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGatata                                                                                       <  1:4329831/150‑1 (MQ=255)
                                                                 gaagaGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGatatata                                                                                     >  2:876150/1‑150 (MQ=255)
                                                                 gaagaGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGatatata                                                                                     <  1:1002089/150‑1 (MQ=255)
                                                                 gaagaGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGatatata                                                                                     >  1:4910506/1‑150 (MQ=255)
                                                                  aagaGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGatatatat                                                                                    >  2:4903631/1‑150 (MQ=255)
                                                                   agagagGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGAtatatata                                                                                   >  2:343689/1‑150 (MQ=255)
                                                                     agagGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCACAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATTCAACAGCGGACGATATATATaaa                                                                                 <  1:539056/150‑1 (MQ=255)
                                                                     agagGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATaaa                                                                                 >  1:4013233/1‑150 (MQ=255)
                                                                       agGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGc                                                                               <  2:474985/150‑1 (MQ=255)
                                                                            cATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTgg                                                                          >  1:1722816/1‑150 (MQ=255)
                                                                             aTGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGa                                                                         >  2:2838533/1‑150 (MQ=255)
                                                                             aTGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGCACGATATATATAAAGCTCTGGa                                                                         >  1:1305522/1‑150 (MQ=255)
                                                                             aTGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGa                                                                         <  1:404889/150‑1 (MQ=255)
                                                                              tGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGaa                                                                        >  1:4400707/1‑150 (MQ=255)
                                                                                aCGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAgg                                                                      <  1:3020323/150‑1 (MQ=255)
                                                                                 cGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAggg                                                                     <  1:2460828/150‑1 (MQ=255)
                                                                                       gCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAAttt                                                               >  2:3901898/1‑150 (MQ=255)
                                                                                        cTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTc                                                              >  1:2242170/1‑150 (MQ=255)
                                                                                         tGCACATGAACTAAAAGAAGCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTcc                                                             <  1:3819598/150‑1 (MQ=255)
                                                                                         tGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTcc                                                             <  2:3314072/150‑1 (MQ=255)
                                                                                          gCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCt                                                            <  2:1536734/150‑1 (MQ=255)
                                                                                           cccaTGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTa                                                           <  2:1319483/148‑1 (MQ=255)
                                                                                           cacaTGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTa                                                           >  1:806239/1‑150 (MQ=255)
                                                                                               tGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACAt                                                       <  2:4505981/150‑1 (MQ=255)
                                                                                               tGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACAt                                                       <  1:398682/150‑1 (MQ=255)
                                                                                                gAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATg                                                      >  2:4620242/1‑150 (MQ=255)
                                                                                                 aaCTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGa                                                     >  2:801973/1‑150 (MQ=255)
                                                                                                       aagaagCCTTAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAg                                               >  2:1657683/1‑150 (MQ=255)
                                                                                                       aagaagCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTTAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAg                                               >  1:2990620/1‑150 (MQ=255)
                                                                                                        agaagCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGc                                              >  1:1268864/1‑150 (MQ=255)
                                                                                                           agCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGAc                                           <  2:4073780/150‑1 (MQ=255)
                                                                                                           agCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGAc                                           <  2:1584724/150‑1 (MQ=255)
                                                                                                           agCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGAc                                           <  2:524934/150‑1 (MQ=255)
                                                                                                           agCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGAc                                           <  2:251506/150‑1 (MQ=255)
                                                                                                             ccTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACgg                                         <  1:4283386/150‑1 (MQ=255)
                                                                                                                tAGAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGtat                                      >  2:4820410/1‑150 (MQ=255)
                                                                                                                 aGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGtata                                     <  2:4264196/150‑1 (MQ=255)
                                                                                                                 aGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGtata                                     >  2:1147678/1‑150 (MQ=255)
                                                                                                                  gAAACGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGtatat                                    <  1:4502216/150‑1 (MQ=255)
                                                                                                                   aaaCGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGtatata                                   <  2:3055256/150‑1 (MQ=255)
                                                                                                                    aaCGTTGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATaa                                  <  2:2076326/150‑1 (MQ=255)
                                                                                                                     accttGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTACTAACATGAGCGTAGCGACGGTATATAAc                                 <  1:3362051/147‑1 (MQ=255)
                                                                                                                     aCGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAAc                                 <  1:3116224/150‑1 (MQ=255)
                                                                                                                      cGTTGGAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACa                                <  2:4024547/150‑1 (MQ=255)
                                                                                                                         tGAAGGAAACCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAAtt                             <  2:2625391/150‑1 (MQ=255)
                                                                                                                            aGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGc                          >  2:2124407/1‑150 (MQ=255)
                                                                                                                                aaCCGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCgtgt                      >  1:4481253/1‑150 (MQ=255)
                                                                                                                                 aCCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCgtgtg                     >  2:1754495/1‑150 (MQ=255)
                                                                                                                                  ccGGAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCgtgtgt                    <  1:1590536/150‑1 (MQ=255)
                                                                                                                                     gAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTcc                 <  2:2378508/150‑1 (MQ=255)
                                                                                                                                     gAGTTCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTcc                 <  1:106675/150‑1 (MQ=255)
                                                                                                                                     gAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTcc                 <  2:676583/150‑1 (MQ=255)
                                                                                                                                     gAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACAGTATATAACAATTTGCGTGTGTTcc                 <  1:4262599/150‑1 (MQ=255)
                                                                                                                                      aGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTCCg                <  1:3169286/150‑1 (MQ=255)
                                                                                                                                        ttCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTCCggg              <  2:1080767/150‑1 (MQ=255)
                                                                                                                                        ttCGCATTACTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTCCggg              <  1:2596012/150‑1 (MQ=255)
                                                                                                                                             aTTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGGCGGTATATAACAATTTGCGTGTGTTCCGGGAATCa         >  2:2162345/1‑150 (MQ=255)
                                                                                                                                             aTTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTCCGGGAATCa         >  2:577737/1‑150 (MQ=255)
                                                                                                                                             aTTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACCCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTCCGGGAATCa         >  1:1361764/1‑150 (MQ=255)
                                                                                                                                                aCTCTTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTCCGGGAATCAGGt      <  2:925707/150‑1 (MQ=255)
                                                                                                                                                aCTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTCCGGGAATCAGGt      <  1:4105583/150‑1 (MQ=255)
                                                                                                                                                    ctcAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGAAGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTCCGGGAATCAGGTTTgg  >  1:2486116/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTGACAAATTATCGTAGAAAGAGTTACACTAATTATAAACATTACAATGTAAGAATTTTTTTTAGAAGAGAGGTGCATGACGGATGGCTGCACATGAACTAAAAGAAGCCTTAGAAACGTTGAAGGAAACCGGAGTTCGCATTACTCCTCAACGTCATGCGATTCTGGAATATCTCGTTAACTCTATGGCTCATCCAACAGCGGACGATATATATAAAGCTCTGGAAGGGAAATTTCCTAACATGAGCGTAGCGACGGTATATAACAATTTGCGTGTGTTCCGGGAATCAGGTTTGG  >  BS168_reference_genome/944403‑944700

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